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digital  receivers  that  use  sophisticated  signal  processing  are  required.  This  thesis  de¬ 
scribes  a  novel  signal  processing  architecture,  and  shows  simulation  results  for  a  number 
of  LPI  waveforms.  The  LPI  signal  detection  receiver  is  based  on  Quadrature  Mirror  Filter 
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EXECUTIVE  SUMMARY 


Low  Probability  of  Intercept  (LPI)  radars  are  a  class  of  radar  systems  that  have 
certain  performance  characteristics  that  make  them  nearly  undetectable  with  most  com¬ 
mon  digital  intercept  receivers.  The  tenn  low  probability  of  intercept  (LPI)  is  that  prop¬ 
erty  of  an  emitter  that  because  of  its  low  power,  wide  bandwidth,  frequency  variability,  or 
other  design  attributes,  makes  it  difficult  to  be  detected  or  identified  by  means  of  passive 
intercept  devices  such  as  radar  warning,  electronic  support  and  electronic  intelligence  re¬ 
ceivers.  New  signal  processing  techniques  are  required  in  order  to  detect  LPI  radar  wave¬ 
forms.  Consequently,  LPI  radar  systems  can  operate  undetected  until  the  intercept 
receiver  is  much  closer  than  the  radar’s  target  detector.  As  a  result,  the  surprise  factor 
remains  one  of  the  most  important  elements  of  military  tactics,  which  poses  a  significant 
problem  for  the  enemy. 

To  detect  LPI  radars,  new  direct  digital  receivers  using  sophisticated  signal  proc¬ 
essing  are  required.  This  thesis  describes  a  type  of  spread  spectrum  intercept  receiver  and 
implements  a  simulation  of  it.  The  LPI  signal  detection  receiver  is  based  on  Quadrature 
Mirror  Filter  Bank  (QMFB)  Tree  processing.  The  receiver  uses  orthogonal  wavelet  tech¬ 
niques  and  a  QMFB  tree  structure  to  decompose  the  input  waveform  into  components 
representing  the  energy  in  rectangular  “tiles”  in  the  time-frequency  plane.  By  analyzing 
the  outputs  at  different  layers  of  the  tree  it  is  possible  to  do  feature  extraction,  identify 
and  classify  the  LPI  waveform  parameters,  and  distinguish  between  the  various  LPI  sig¬ 
nal  modulations.  Input  continuous  waveforms  used  as  input  signals  to  the  detection  algo¬ 
rithm  include  Frequency  Modulated  Continuous  Wave,  Polyphase  Codes,  Costas  Codes 
and  Frequency  Shift  Keying/Phase  Shift  Keying  waveforms.  The  output  matrices  result¬ 
ing  from  the  most  relevant  layers  of  the  QMFB  tree  processing  are  examined  and  the  LPI 
modulation  parameters  are  extracted  under  various  signal-to-noise  ratios. 

The  overall  architecture  of  the  QMFB  presented  in  this  thesis  is  a  good  alternative 
to  an  intercept  receiver,  particularly  for  detecting  hopped  LPI  spread  spectrum  radar  sig- 
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nals,  as  the  FSK/PSK  Costas.  The  feature  extraction  gives  a  time-frequency  plane  with 
good  estimation  of  position,  center  frequency  and  time  duration  of  the  signals  processed. 


I. 


INTRODUCTION 


A.  LPI  RADAR  SIGNALS  AND  THEIR  DETECTION 


In  this  chapter  the  LPI  radar  signals  and  their  detection  are  explored  giving  a 
background  of  the  problem  and  defining  the  research  objective. 

1.  To  See  and  Not  Be  Seen 

New  advanced  components  and  signal  processing  algorithms  have  been  devel¬ 
oped  by  modem  electronic  technology.  Some  of  these  new  electronic  components  and 
signal  processing  techniques  are  applied  in  radar  systems  with  complex  waveforms.  As  a 
result,  the  performance  of  these  radars  improves  significantly.  The  appearance  of  modern 
Electronic  Support  (ES),  Radar  Warning  Receivers  (RWR),  and  Anti  Radiation  Missiles 
(ARM)  constitute  the  most  important  threats  to  radar  operation  on  the  battlefield.  Conse¬ 
quently,  modem  radars  should  have  a  Low  Probability  of  Intercept  (LPI)  capability  in 
order  to  hide  from  and  survive  an  enemy  attack. 

LPI  is  considered  an  important  tactical  requirement  and  is  being  specified  by 
many  military  personnel  using  radar  today.  An  LPI  radar  has  specific  design  features  that 
render  it  difficult  to  detect.  Low  probability  of  intercept  is  that  property  of  an  emitter  that 
because  of  its  low  power,  wide  bandwidth,  frequency  variability,  or  other  design  attrib¬ 
utes,  makes  it  difficult  to  be  detected  or  identified  by  means  of  passive  intercept  receiver 
devices.  LPI  radars  attempt  to  detect  targets  at  longer  ranges  than  intercept  receivers  at 
the  target  can  detect  the  radar.  Thus,  the  objective  of  an  LPI  radar  is  “To  See  and  Not  Be 
Seen,”  or,  “To  Detect  and  Not  Be  Detected”  [1]. 

2.  The  Principle  of  LPI  Radar 

In  order  to  accomplish  interception  by  an  intercept  receiver,  the  detection  range  of 
the  radar  Rr  should  be  longer  than  the  range  Rf  at  which  the  intercept  receiver  can  detect 
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the  radar  emission  as  shown  in  Figure  1.  From  Figure  1,  a  range  factor  a  can  be  defined 
as: 


Rr 


(1.1) 


If  a  >  1 ,  the  radar  may  be  detected  at  a  range  greater  than  the  range  the  radar  can  detect 
targets.  If  a  <  1 ,  the  radar  can  detect  targets  at  a  range  greater  than  the  range  of  intercept. 
This  latter  type  of  radar  is  an  LPI  radar. 


Figure  1  LPI  radar,  target  and  intercept  receiver  configuration.  After  [2]. 


It  should  be  noted  that  Rr  and  R,  only  occur  when  certain  detection  probability 

and  false  alarm  probability  conditions  are  satisfied.  In  fact,  the  so-called  LPI  performance 
is  a  probabilistic  event.  LPI  radars  may  not  be  detected  by  certain  low  sensitivity  inter¬ 
cept  receivers,  but  high  sensitivity  intercept  receivers  may  have  no  problem.  Although, 
the  LPI  radars  take  advantage  in  the  signal-to-noise  ratio  (SNR)  conditions,  hiding  and 
spreading  the  transmitted  waveform  under  the  noise  level  in  a  wide  bandwidth. 

In  the  case  of  an  electronic  attack,  just  detecting  the  energy  will  not  help  program 
the  waveform  generator  of  a  jamming  system  to  perform  the  attack.  For  this  the  modula¬ 
tion  parameters  are  necessary  and  sophisticated  signal  processing  required  to  defining  the 
transmission  waveform  parameters  of  the  jamming  signal.  [2]. 
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3.  Characteristics  of  LPI  Radar  Signals  and  Their  Detection 

LPI  radars  have  many  combined  features  that  help  prevent  their  parameter  detec¬ 
tion  by  modem  intercept  receivers.  These  features  are  centered  on  the  antenna  and  trans¬ 
mitter.  The  first  antenna  characteristic  is  a  low  side-lobe  transmit  pattern.  The  use  of  low 
side-lobe  pattern  reduces  the  possibility  of  an  intercept  receiver  detecting  the  radio  fre¬ 
quency  (RF)  emissions  throughout  the  side-lobe  structure  of  the  antenna  pattern.  The 
second  antenna  characteristic  is  the  scan  pattern,  which  is  precisely  controlled  to  limit  the 
intercept  receiver  time  to  short  and  infrequent  intervals  (a  periodic  scan  cycle).  Scan 
methodologies  can  also  be  added  to  help  confuse  identifications  of  interest  if  they  occur. 
For  example,  scan  techniques  that  attempt  to  confuse  identification  by  an  intercept  re¬ 
ceiver  might  include  amplitude  modulation  of  a  monopulse  array  at  conical  scan  frequen¬ 
cies.  [1]. 

The  main  drawback  of  the  coherent  pulse  train  waveform  is  the  high  ratio  of  peak 
to  average  power  emitted  by  the  transmitter.  This  average  power  determines  the  detection 
characteristics  of  the  radar.  For  high  average  power,  the  short  pulse  (high  resolution) 
transmitter  must  have  a  high  peak  power,  which  necessitates  vacuum  tubes  and  high  volt¬ 
ages.  The  high  peak  power  transmissions  can  also  easily  be  detected  by  ES  receivers.  In 
modulated  continuous  wave  (CW)  signals,  however,  the  peak-to-average  power  ratio  is 
one  (100%  duty  cycle)  which  allows  a  considerable  lower  transmit  power  to  maintain  the 
same  detection  performance.  Also,  solid  state  transmitters  can  be  used  which  are  lighter 
in  weight.  A  comparison  is  shown  in  Figure  2. 


Power 


▲ 

Pulse  radar 
^  High  peak  power 

Small  duty  cycle 

^  CW  radar 

Low  continuous  power 
100%  duty  cycle 

/ 

^  Time 

Figure  2  Comparison  of  pulsed  radar  and  CW  Radar.  After  [  1  ] . 
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LPI  emitters  use  CW  signals  with  large  bandwidths  and  small  range  resolution 
cells,  resulting  in  a  moderate  size  target  return  being  spread  over  multiple  cells.  Periodi¬ 
cally  modulated  CW  signals  are  extensively  used  in  LPI  radar  and  are  ideally  suited  for 
pulse  compression.  They  achieve  a  unity  peak-to-average  power  ratio.  The  fact  that  the 
transmitted  signal  is  continuous  does  not  imply  that  the  portion  processed  in  order  to  de¬ 
tect  a  target  is  infinitely  long.  There  are  physical  constraints  (the  illumination  time,  for 
example)  and  processor  constraints  that  can  cause  this.  Fast  Fourier  Transform  (FFT) 
processors  (for  frequency-modulated  waveforms)  and  finite  duration  coherent  correlation 
processors  (for  phase  modulation  waveforms)  are  among  the  most  often  used  as  well  as 
combinations  of  both  [1]. 

There  are  many  wideband  modulation  techniques  available  from  the  transmitter 
function  that  provides  an  LPI  transmit  waveform.  Any  change  in  the  radar’s  signature 
(pulse  repetition  frequency  (PRF),  pulse  width  (PW),  carrier  frequency,  polarization,  scan 
modulation)  can  help  confuse  an  intercept  receiver  and  make  identification  difficult. 
Wideband  modulation  techniques  include: 

Linear  and  Non-Linear  frequency  modulation, 

Costas  arrays,  frequency  hopping, 

Phase  modulation  (polyphase  coding), 

Combined  phase  shift  keying,  frequency  shift  keying  (PSK,  FSK),  and 
Pseudo-noise  modulation. 

The  radiated  energy  is  spread  over  a  wide  frequency  range  in  a  noise-like  manner  with 
most  of  these  modulation  techniques. 

Another  feature  of  an  LPI  transmitter  is  power  management.  This  is  one  of  the 
benefits  to  using  solid-state  radars.  The  ability  to  control  the  signature  emitted  by  the  ar¬ 
ray  is  used  to  limit  the  emissions  to  the  appropriate  Range/Radar  Cross  Section  require¬ 
ment.  The  emissions  are  also  limited  in  time  (short  dwell  time).  With  the  use  of  wideband 
CW  emissions  it  is  only  necessary  to  transmit  a  few  watts  instead  of  tens  of  kilowatts  of 
peak  power  required  by  pulsed  radars  with  similar  perfonnance.  It  is  important  to  note 
that  the  radar’s  ability  to  detect  targets  depends  not  on  the  waveform  characteristics  but 
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on  the  transmitted  energy  density  returned  to  the  radar  from  the  target.  The  main  objec¬ 
tive  of  a  LPI  radar  is  to  operate  under  low  (SNR)  conditions  so  that  integration  of  the  sig¬ 
nal  over  several  contiguous  range  cells  can  be  used  to  detect  and  track  the  targets  of 
interest.  This  integration  leads  to  an  operative  resolution  cell  that  is  larger  than  FFT  bin 
[1]. 

Another  approach  to  achieving  a  lower  probability  of  interception  is  to  interleave 
the  LPI  radar  mode  with  an  infrared  sensor  (dual  mode  approach).  The  amount  of  time 
that  the  RF  transmitter  is  radiating  is  therefore  reduced. 

4.  Research  Objective 

This  thesis  developed  a  simulation  of  a  method  to  decompose  a  signal,  using  or¬ 
thogonal  basis  functions  and  a  quadrature  mirror  filter  bank  (QMFB)  tree  configuration. 
Detailed  information  about  processed  LPI  signals  were  then  extracted.  The  receiver  for 
this  approach  is  shown  in  Figure  3.  In  this  block  diagram  a  received  waveform  is  band¬ 
pass  filtered  and  sampled  at  the  Nyquist  rate.  The  digital  sequence  is  then  fed  to  the 
QMFB  tree  where  it  is  decomposed  and  matrices  of  weights  are  output  from  each  layer. 
These  weights  are  then  squared  to  produce  coefficients  representing  the  energy  in  each 
portion  of  the  waveform. 

The  information  is  then  analyzed  to  determine  all  the  relevant  parameters  from  the 
waveform.  Next,  the  output  of  the  analyzer  is  a  list  of  all  parameters  already  extracted. 
The  classifier  takes  the  list  from  the  analyzer,  determines  which  parameters  belong  to 
common  transmitters,  and  outputs  a  list  of  possible  transmitters,  their  types,  and  parame¬ 
ters.  Further  filtering  could  be  done  at  this  point,  eliminating  false  alanns  and  signals  that 
are  not  of  interest  to  the  interceptor.  The  classifier  may  even  be  adaptive  and  change  clas¬ 
sification  criteria  based  on  current  and  previous  input  and  results  [4]. 

The  research  in  this  thesis  was  primarily  concerned  with  the  receiver’s  QMFB 
tree  and  the  analyzer  block.  It  was  also  focused  on  developing  algorithms  to  detect  and 
determine  the  features  for  each  type  of  LPI  signals.  The  perfonnance  of  the  QMFB  under 
various  SNR  is  also  investigated  in  the  presence  of  additive  White  Gaussian  Noise 
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(WGN).  The  only  signals  assumed  to  be  present  are  single  LPI  signals  at  a  time.  The  in¬ 
put  bandpass  filter  and  sampler  in  the  receiver  shown  in  Figure  3  are  assumed  to  be  ideal. 

This  thesis  does  not  examine  the  classifier  function  in  any  detail.  The  primary  job 
of  the  classifier,  however,  is  to  take  the  list  of  parameters,  determine  the  most  likely 
number  of  transmitters,  and  group  the  parameters  to  the  transmitters.  This  specific  task 
was  not  part  of  this  thesis. 


Figure  3  LPI  receiver  block  diagram.  After  [4]. 


B.  PRINCIPAL  CONTRIBUTIONS 


The  first  phase  of  this  thesis  involved  doing  a  literature  research  to  reach  a  good 
understanding  on  LPI  radar  signals,  quadrature  mirror  filter  (QMF),  signal  processing  and 
MATLAB®  coding  [5].  This  was  found  through  different  books,  papers  and  courses  from 
databases  as  IEEE,  some  signal  processing  books,  and  courses  taken  during  the  studies  in 
Electrical  and  Computer  Department  of  NPS,  respectively.  The  next  phase  was  to  define 
the  QMFB  signal  processing  and  develop  a  new  program  to  examine  LPI  radar  signals. 
With  the  new  program  in  MATLAB®,  test  signals  were  used  to  verify  that  the  signal 


6 


processing  worked  correctly.  Next,  a  list  of  different  LPI  signals  were  analyzed  with  the 
program,  which  displays  time-frequency  output  plots  from  which  it  is  possible  to  extract 
the  different  parameters  of  the  LPI  signals.  A  comparison  of  the  results  of  the  processing 
done  by  the  QMFB  tree  on  all  the  signals  is  presented  as  a  QMFB  processing  perform¬ 
ance  on  LPI  radar  signals.  Finally,  a  summary  of  the  thesis,  conclusions  of  the  QMFB 
processing  and  recommendations  to  continue  with  this  research  are  indicated. 


C.  THESIS  OUTLINE 


Chapter  II  gives  a  brief  overview  on  the  theory  and  background  on  time- 
frequency  signal  processing  and  wavelets  decomposition.  Chapter  III  presents  the  LPI 
signals  processed  by  the  QMFB  tree  where  each  signal  is  analyzed  and  the  parameters  are 
extracted  from  the  most  relevant  QMFB  output  layer  finishing  with  an  analysis  and  com¬ 
parison  of  the  signal  processing  results.  Finally,  Chapter  IV  gives  conclusions  from  the 
QMFB  processing,  and  recommendations  for  future  research  are  indicated. 
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II.  QUADRATURE  MIRROR  FILTER  BANK 


This  thesis  investigates  an  LPI  intercept  receiver  based  on  a  linear  decomposition 
of  the  received  waveform  through  a  Quadrature  Mirror  Filter  Bank  (QMFB)  Tree  using 
wavelet  filters.  The  receiver  provides  good  detection  of  the  LPI  signal  parameters  in  an 
effort  to  distinguish  between  the  different  modulations.  This  chapter  presents  the  QMFB 
theory,  discussing  the  mathematical  background  for  signal  detection  and  the  waveform 
decomposition  using  Wavelets,  and  finishes  with  a  description  of  the  MATLAB®  codes 
used  to  implement  the  QMFB  tree. 

A.  BACKGROUND 

1.  Time-frequency  Decomposition 

Various  methods  of  decomposing  a  waveform  on  the  time-frequency  plane 
have  recently  been  investigated.  The  most  common  methods  can  be  divided  into  linear 
and  bilinear  transforms.  The  Short  Time  Fourier  and  Wavelet  Transform  are  examples  of 
linear  transforms,  and  the  Wigner  transform  is  example  of  a  bilinear  transfonn.  Wigner 
transforms  are  called  bilinear  because  the  input  waveform  appears  twice  in  the  develop¬ 
ment  of  the  transfonn.  Better  resolution  occurs  in  the  time-frequency  plane  than  with  lin¬ 
ear  techniques.  However,  the  computational  burden  is  greatly  increased  resulting  in  other 
side  effects.  At  the  same  time,  is  not  easy  to  find  studies  specifically  using  this  transform 
for  the  detection  of  LPI  waveforms.  With  the  computational  complexity  and  the  possibil¬ 
ity  of  confusing  cross-terms,  the  bilinear  transform  is  considered  further  in  [18],  Linear 
transfonns  of  the  continuous  time  signal  /'(/)  have  the  following  form 

ak=\fmk(t)dt  (2.1) 

where  is  the  basis  set,  t  is  the  time  index,  and  k  is  the  function  index 

The  Fourier  transform,  for  example,  has  a  basis  set  consisting  of  sines  and  co¬ 
sines  of  frequency  Ink  .  The  basis  functions  are  said  to  be  orthonormal  if, 
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jo  (t-k)®(t-l)dt 


1  k=l 
Ok*/. 


[6] 


(2.2) 


If  the  basis  functions  are  orthonormal,  there  is  no  redundancy  in  the  representation  of  the 
signal  f(t) .  It  is  possible  to  sample  the  input  waveform  at  the  Nyquist  rate  and  retain  all 
the  information.  In  that  case,  the  time  variable,  t ,  in  (2.1)  and  (2.2)  should  be  considered 
to  be  discrete  and  the  integral  should  be  replaced  with  summations.  The  basis  functions 
are  said  to  be  orthogonal  if 


j0(t)0(t-k)dt  =  ES(k) 


j  E  if  ik  =  0 
[0  otherwise 


(2.3) 


where  E  stands  for  energy  of  O(t)  [7]. 

Wavelet  basis  functions  are  effectively  non-zero  for  only  a  finite  time  interval, 
and  can  be  designed  to  satisfy  (2.2).  References  [4]  and  [6]  demonstrate  that  these  or¬ 
thogonal  Wavelets  can  be  implemented  using  Quadrature  Mirror  Filters:  filter  pairs  that 
are  designed  to  divide  the  input  signal  energy  into  two  orthogonal  components  based  on 
the  frequency.  The  basis  function  becomes  a  contracted  wavelet,  or  a  short  high  fre¬ 
quency  function  as  it  is  in  Figure  4  (a)  where  ,  and  the  wavelet  transform  divides  the 
time-frequency  plane  into  tiles  as  shown  in  Figure  4  (b).  Here  the  area  of  each  tile  repre¬ 
sents  (approximately)  the  energy  within  the  function  (rectangular  regions  of  the  fre¬ 
quency  plane).  Note  that  not  all  of  the  signal’s  energy  can  be  located  in  a  single  tile 
because  it  is  impossible  to  concentrate  the  function’s  energy  simultaneously  in  frequency 
and  time.  A  characteristic  of  the  Wavelet  transform  is  that  the  tiles  become  shorter  in 
time  and  occupy  a  larger  frequency  band  as  the  frequency  is  increased. 
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(a) 


Frequency 

A 


^  Time 

(b) 


Figure  4  Basis  functions  and  time-frequency  resolution  of  the  wavelet  transform. 

(a)  Basis  functions,  (b)  Coverage  of  time-frequency  plane. 

Using  Wavelet  techniques  to  develop  an  appropriate  basis  set  and  a  QMFB  for  imple¬ 
mentation,  it  is  possible  to  decompose  the  waveform  in  such  a  way  that  the  tiles  have  the 
same  dimensions  regardless  of  the  frequency  [4,  10].  Since  the  transform  is  linear,  there 
is  a  fundamental  limit  on  the  minimum  area  of  these  tiles.  However,  the  nature  of  the 
QMFB  configuration  is  such  that  each  layer  outputs  a  matrix  of  coefficients  for  tiles  that 
are  twice  as  long  (in  time)  and  half  as  tall  (in  frequency)  as  the  tile  in  the  previous  layer. 
By  properly  comparing  these  matrices,  it  is  possible  to  extract  signal  features  using  both 
fine  frequency  and  fine  time  resolutions.  Parameters  such  as  bandwidth,  center  frequency, 
energy  distribution,  phase  modulation,  signal  duration  and  location  in  the  time-frequency 
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plane  can  be  determined  using  these  techniques,  making  them  valuable  for  intercept  re¬ 
ceivers  in  order  that  they  can  detennine  the  number  of  transmitters  present  and  which 
types  are  in  operation. 

2.  Decomposition  of  Waveforms 

a.  Short  Time  Fourier  Transform 

Complex  sinusoids  are  used  by  the  Fourier  transform  to  perfonn  the 
analysis  of  signals  using  appropriate  basis  functions.  This  approach  is  difficult  due  to  the 
infinite  extent  of  the  basis  functions  as  any  time-local  infonnation,  such  as  an  abrupt 
change  in  the  signal,  is  spread  out  over  the  entire  frequency  axis.  This  problem  has  been 
addressed  by  introducing  windowed  complex  sinusoids  as  basis  functions  [10].  This  leads 
to  the  doubly  indexed  windowed  Fourier  transform: 

oo 

Xwf(g>,t)=  J  e~ia>,w(t  -  r)x(t)  dt  (2.4) 

-oo 

where  w(t  -  z)  constitutes  an  appropriate  window.  XWF(co,  r)  is  the  Fourier  transform  of 
x(t)  windowed  with  vv'(-)  shifted  by  r  .  The  major  advantage  of  the  windowed  or  short- 
time  Fourier  transfonn  (STFT)  is  that  if  a  signal  has  most  of  its  energy  in  a  given  time 
interval  [-T,r]  and  frequency  interval  [-0,0],  then  its  STFT  will  be  localized  in  the 

region  [-T,r]x  [-0,0]  and  will  be  close  to  zero  in  time  and  frequency  intervals  where 

the  signal  has  little  energy  [10],  A  limitation  of  the  STFT  is  that,  because  a  single  win¬ 
dow  is  used  for  all  frequencies,  the  resolution  of  the  analysis  is  the  same  at  all  locations 
in  the  time-frequency  plane.  The  possibility  of  having  arbitrarily  high  resolution  in  both 
time  and  frequency  is  thus  excluded. 

b.  Wavelet  Transform 

By  varying  the  window  used,  resolution  in  time  can  be  traded  for  resolu¬ 
tion  in  frequency.  To  isolate  discontinuities  in  signals  it  is  possible  to  use  some  basis 
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functions,  which  are  very  short,  while  longer  ones  are  required  to  obtain  a  fine  frequency 
analysis.  One  method  to  achieve  this  is  to  have  short  high-frequency  basis  functions,  and 
long  low-frequency  basis  functions  [10].  This  is  exactly  what  the  wavelet  transfonn 
achieves  where  the  basis  functions  are  obtained  from  a  single  prototype  wavelet  by  trans¬ 
lation  and  dilation/contraction: 


K,(0 


h 


t-b 
V  a  J 


(2.5) 


where  a  is  a  positive  real  number  and  b  is  a  real  number.  For  large  a ,  the  basis  function 
becomes  a  stretched  version  of  the  prototype  wavelet  (low  frequency  function);  while  for 
small  a  ,  the  basis  function  becomes  a  contracted  wavelet  (short  high  frequency  function) 
as  it  is  in  Figure  4  (a).  The  wavelet  transfonn  (WT)  is  defined  as 


(^’  ^0 


t  —  b ^ 

Cl  j 


x(t)  dt 


(2.6) 


The  time-frequency  resolution  of  the  WT  involves  a  different  tradeoff  than  that  used  by 
the  STFT.  At  high  frequencies,  the  WT  is  sharper  in  time;  while  at  low  frequencies,  the 
WT  is  sharper  in  frequency  [10].  If  the  parameters  (co,  r)  and  (a,b)  are  continuous,  the 
STFT  from  Equation  (2.4),  and  the  WT,  from  Equation  (2.6),  are  highly  redundant. 
Therefore,  the  transfonns  are  usually  evaluated  on  a  discrete  grid  on  the  time-frequency 
and  time-scale  planes,  respectively,  corresponding  to  a  discrete  set  of  continuous  basis 
functions.  At  this  point  a  grid  is  needed  such  that  the  set  of  basis  functions  constitutes  an 
orthonormal  basis  (no  redundancy).  Unfortunately,  for  the  STFT,  this  only  occurs  if 
w(t  -  r)  is  badly  localized  in  either  time  or  frequency.  This  is  the  reason  that  the  STFT  is 
usually  oversampled,  as  a  redundant  set  of  points  is  used,  so  that  better  behaved  window 
functions  can  be  used.  In  the  wavelet  transform  case,  however,  it  is  possible  to  design 
“practical”  functions  h(-)  such  that  the  set  of  translated  and  scaled  versions  of  /z(-)  forms 
an  orthonormal  basis.  “Practical”  means  that  the  function  should  be  at  least  continuous, 
and  perhaps  with  continuous  derivatives  also.  Making  the  translation  and  dila¬ 
tion/contraction  discrete  parameters  of  the  wavelets  in  Equation  (2.5): 

Kn (0  =  aom'2  ■  h(aomt  ~ nbo )  (2-7) 
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where  m  ,  n  are  integers,  a  >  1 ,  and  b0  ^  0 ,  which  correspond  to  a  =  a™  and  b  =  na™b0 . 
Note  that  the  translation  step  depends  on  the  dilation,  since  long  wavelets  are  advanced 
by  large  steps,  and  short  ones,  by  small  steps.  On  this  discrete  grid,  the  wavelet  transform 
is  thus 

oo 

Xw(m,n)  =  a0"'12  J  h(a0mt  - nb0)x(t)  dt  (2.8) 


Of  particular  interest  is  the  discretization  on  a  dyadic  grid,  which  occurs  for  a0  =  2, 
b0  =  1 .  It  is  possible  to  construct  functions  h(-)  so  that  the  set 

Kn (0  =  ««'" 2  •  h(a0mt  ~ nb0),  (2.9) 


where  m  ,  n  are  integers,  a0  =  2,  and  b0=  1  is  orthonormal  [10]. 


A  classic  example  is  the  Haar  basis,  which  is  not  continuous,  but  is  of  interest  because  of 
its  simplicity,  where: 


h(t)  =  { 


1 

-1 

0 


0  <  t  <  1/2 
1/2  <  ?  <  1 

otherwise. 


(2.10) 


The  orthonormality  is  easily  verified  since  at  a  given  scale,  translates  are 
nonoverlapping,  and  because  of  the  scale  change  by  2,  the  basis  functions  are  orthonor¬ 
mal  across  scale.  The  Haar  basis  is  shown  in  Figure  5.  However,  the  Haar  function  is  dis¬ 
continuous  and  is  not  generally  appropriate  for  signal  processing.  A  continuous  set  of 
basis  functions  must  be  used  [6]. 


h(t) 

1 

1  Time 

0 

-1 

Figure  5  Mother  wavelet  for  the  Haar  basis  set. 
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From  a  signal  processing  point  of  view,  a  wavelet  is  a  bandpass  filter.  Therefore  the 
wavelet  transform  can  be  interpreted  as  a  constant-Q  filtering  with  a  set  of  subband  fil¬ 
ters,  followed  by  a  sampling  at  the  respective  Nyquist  frequencies  corresponding  to  the 
bandwidth  of  the  particular  subband  [6]. 

3.  Discrete  Two- Channel  Quadrature  Mirror  Filter  Bank 

In  many  applications,  a  discrete-time  signal  x[n\  is  first  split  into  a  number  of 
subband  signals  {vt[n]}  by  means  of  an  analysis  filter  bank.  The  subband  signals  are  then 

processed  and  finally  combined  by  a  synthesis  filter  bank,  resulting  in  an  output  signal 
y\ri\ .  If  the  subband  signals  are  bandlimited  to  frequency  ranges  much  smaller  than  that 
of  the  original  input  signal,  they  can  be  down-sampled  before  processing.  Due  to  the 
lower  sampling  rate,  the  processing  of  the  down-sampled  signals  can  be  carried  out  more 
efficiently.  After  processing,  these  signals  are  up-sampled  before  being  combined  by  the 
synthesis  bank  into  a  higher-rate  signal.  The  combined  structure  employed  is  called  a 
QMFB.  If  the  down-sampling  and  the  up-sampling  factors  are  equal  to  or  greater  than  the 
number  of  bands  of  the  filter  bank,  then  the  output  y\n\  can  be  made  to  retain  some  or  all 
of  the  characteristics  of  the  input  x[n\  by  properly  choosing  the  filters  in  the  structure 
[11]. 


a.  The  Filter  Bank  Structure 

Figure  6  shows  the  basic  two-channel  QMFB.  Here,  the  input  signal  x[n] 
is  first  passed  through  a  two-band  analysis  filter  bank  containing  the  filters  H0(z)  and 
//,  (z) ,  which  typically  have  lowpass  and  highpass  frequency  responses,  respectively, 
with  a  cutoff  frequency  at  nil  as  indicated  in  [1 1]. 
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Figure  6  The  two-channel  quadrature  mirror  filter  (QMF)  bank.  After  [11]. 


The  subband  signals  { vk  [/?]}  are  then  down-sampled  by  a  factor  of  2  in  the  “signal  analy¬ 
sis  section”  to  be  transmitted  to  the  “signal  synthesis  section”  where  the  signals  will  be 
up-sampled  by  a  factor  of  2  and  passed  through  a  two-band  synthesis  filter  bank  com¬ 
posed  of  the  filters  G0(z)  and  Gt(z)  whose  outputs  are  then  added  yielding  y[n] .  It  fol¬ 
lows  from  the  figure  that  the  sampling  rates  of  the  input  signal  x[n\  and  output  signal 
y\n\  are  the  same.  The  analysis  and  the  synthesis  filters  in  the  QMFB  are  chosen  so  as  to 
ensure  that  the  reconstructed  output  y[/z]  is  a  reasonable  replica  of  the  input  x\n\ .  More 
over,  they  are  also  designed  to  provide  good  frequency  selectivity  to  ensure  that  the  sum 
of  the  power  of  the  sub-band  signals  is  reasonably  close  to  the  input  signal  power  [11]. 
For  the  signal  analysis  section,  H0(co )  is  a  lowpass  filter  and  Hx(co)  is  a  mirror-image 
highpass  filter  as  is  shown  by  Figure  7. 


Figure  7  Typical  frequency  response  of  the  analysis  filters.  After  [11], 
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4.  Filtering  and  Arbitrary  Tiling 


a.  Filtering  the  Lowpass  Component 

Finite  impulse  response  (FIR)  filters  and  decimators  can  be  arranged  in  the 
tree  structure  as  shown  in  Figure  8  to  effect  an  orthogonal  wavelet  decomposition  of  a 
signal  [12].  The  discrete  input  waveform  is  denoted  as  the  sequence  {c0}  and  the  output 

sequences  of  each  branch  are  as  shown  in  the  figure.  Since  each  branch  of  the  tree  down 
samples  by  2,  each  sequence  will  have  half  as  many  elements  as  the  preceding  sequence. 


Figure  8  Wavelet  filter  bank  tree  [12]. 


A  filter  tree  using  the  same  orthogonal  pair  of  filters  throughout,  and  with 
equal  length  branches,  as  in  Figure  8,  yields  a  rectangular  tiling  diagram.  The  time- 
frequency  tiling  diagram  shown  in  Figure  9,  is  one  method  that  can  be  used  to  describe 
this  decomposition.  Time-frequency  tile  is  “the  region  in  the  plane,  which  contains  most 
of  that  function’s  energy’’  [4].  However,  not  all  of  a  function’s  energy  can  be  located  in  a 
tile  because  it  is  impossible  to  fully  concentrate  energy  simultaneously  in  time  and  fre¬ 
quency. 
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Figure  9  Time-frequency  diagram  for  the  wavelet  filter  bank  tree  [12]. 

b.  Filtering  the  Highpass  Component  ( Arbitrary  Tiling ) 

The  last  section  demonstrated  that  by  cascading  filters  and  filtering  the 
lowpass  component  of  the  previous  output,  a  tiling  with  finer  frequency  resolution  at 
lower  frequencies  was  achieved.  In  many  cases  concerning  detection,  however,  this  is  not 
desirable.  Many  man-made  signals,  for  example,  will  have  a  constant  bandwidth  over  a 
wide  range  of  center  frequencies.  Instead,  an  arbitrary  tiling  is  desired  to  meet  specific 
requirements  for  the  type  of  signal  to  be  detected.  In  order  to  accomplish  this,  it  is  possi¬ 
ble  to  modify  the  filter  bank  presented  in  the  last  section. 

Consider  the  cascading  filter  diagram  in  Figure  10  where,  instead  of  filter¬ 
ing  the  lowpass  output  of  each  stage,  the  highpass  output  is  filtered.  Again,  the  input  se¬ 
quence  is  split  at  each  stage  into  high-frequency  and  low-frequency  orthogonal 
sequences,  and  the  tiling  diagram  shown  in  Figure  1 1  is  therefore  obtained  and  is  some¬ 
times  referred  to  as  “Wavelet  Packet  Tiling”  [4]. 
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Figure  10  Wavelet  filter  bank  tree. 


Notice  the  second  and  third  layer  seem  to  be  flipped  in  Figure  10.  The  fig¬ 
ure  is  drawn  so  that  the  output  sequence  at  the  top  of  the  drawing  contains  the  highest 
frequency  components  of  the  input  sequence.  This  will  be  important  later.  To  understand 
why  they  are  flipped,  consider  the  aliased  frequency  spectrum  of  the  filters,  shown  in 
Figure  12.  The  output  from  the  G  filter  in  the  first  layer  contains  the  higher  frequency 
components  of  the  original  sequence,  but  shifted,  so  it  is  actually  the  DC  component  of 
the  output  of  G.  The  result  is  that  the  output  of  G  is  frequency  reversed,  much  like  the 
lower  sideband  of  a  single  sideband  communications  system.  Of  course,  similarly  further 
down  the  cascade  will  “unflip”  the  signal.  A  simple  rule  exists  that  maintains  the  integrity 
of  the  output.  This  will  be  discussed  shortly. 
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Figure  1 1  Time-frequency  diagram  for  the  wavelet  filter  bank. 
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It  is  possible  to  create  another  tiling  scheme  by  combining  the  Wavelet  fil¬ 
ter  bank  and  Wavelet  tiling,  as  demonstrated  in  Figure  13.  The  construction  rule  for  this 
figure,  in  order  to  keep  the  higher  frequency  outputs  of  each  branch  above  the  lower  fre¬ 
quency  outputs,  is  to  count  the  number  of  G  filters  up  to  the  branch.  If  the  number  is 
even,  the  next  G  filter  will  output  the  high  frequencies.  If  odd,  the  next  H  filter  will  out¬ 
put  the  high  frequencies. 


20 


Figure  13  Combining  the  wavelet  filter  bank  and  wavelet  tiling.  (Decimation  by  2  is  in¬ 
cluded  in  each  filter  box). 


5.  Wavelet  Filters 

The  objective  of  this  section  is  to  find  a  FIR  wavelet  filter  that  best  approximates 
the  perfect  time-frequency  tiling  by  trading  off  and  minimizing  the  out-of  time  and  out-of 
frequency  energy  from  a  comparison  in  between  the  Haar  filter  and  the  Sine  filter. 


a.  The  Haar  Filter 

The  Haar  filter  was  discussed  briefly  at  the  end  of  Section  2. a.  It  has  two 

coefficients 


h(0)  =  h(\)  =  /^2 


(2.11) 
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The  energy  in  each  output  value  is  equal  to  the  low-pass  energy  from  the  two  correspond¬ 
ing  input  values  and  from  no  others.  Since  each  of  the  two  input  values  contributes 
equally  to  the  output,  the  pass  region  is  also  flat  along  the  time  dimension.  It  is  the  only 
Wavelet  FIR  filter  that  is  symmetric.  The  filter  perfectly  concentrates  (tiles)  the  input  en¬ 
ergy  in  time.  However,  the  Haar  filter  does  not  tile  well  in  frequency  [4,  15]. 


b.  The  Sine  Filter 


The  Sine  filter  perfectly  concentrates  energy  in  frequency.  It  has  however, 
an  infinite  number  of  coefficients.  This  condition  is  modified  in  this  thesis.  By  starting 
with  the  desired  frequency  response,  which  is  a  flat  passband,  an  infinitely  narrow  transi¬ 
tion,  and  a  zero  across  the  stop  band,  the  inverse  Fourier  Transform  is  taken.  This  will 
result  in  a  sine  function  in  the  time  domain  [4,  15].  Therefore,  the  sine  function  is 


sinc(k)  =  < 


sin(^k) 

nk 

1 


k*  0 
k  =  0 


(2.12) 


Since  the  passband  ranges  from  -n/2  <  co  <  n/2  or  -0.25  </<  0.25,  the  nulls  of  the  sine 
function  will  be  at  2 T  for  a  sampling  period  of  T  to  obtain  the  filter  coefficients.  The 
sine  function  will  be  sampled  at  the  nonnalized  sampling  rate  of  T  =  1  for  a  situation 
similar  to  that  shown  in  Figure  14.  One  way  to  sample  the  function  would  be  to  let  the 
main  tap  sample  occur  at  the  center  of  the  main  lobe.  However,  two  main  taps  are  needed 
and  their  sum  needs  to  be  as  large  as  possible.  This  occurs  for  the  sine  function  if  both 
main  tap  samples  are  equally  spaced  about  the  center  of  the  main  lobe.  The  sum  of  the 
square  of  the  coefficients  must  be  unity  also,  which  is  achieved  by  scaling  the  sine  by 
1/ v2  ,  giving 


where  n  is  an  integer. 


h(n )  =  — j=sinc 
V  2 


r  n  +  0.5^ 

V  2  , 


(2.13) 
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Figure  14  Sampling  under  a  sine  envelope. 

In  fact,  this  filter  meets  the  criteria  of  Wavelet  filters.  The  only  problem  is 
that  there  are  an  infinite  number  of  coefficients.  A  small  amount  of  non-orthogonality 
will  occur  when  trying  to  determine  how  to  truncate  this  filter  and  maintain  a  good  fre¬ 
quency  response.  Some  cross-correlation  will  take  place  between  filters. 

If  the  ends  of  the  filter  are  simply  truncated,  some  ripples  in  the  passband 
of  the  frequency  response  will  appear  which  is  sometimes  called  the  "Gibb's  Phenom¬ 
ena".  It  is  a  well-known  result  of  this  type  of  truncation  [12].  Multiplying  the  coefficients 
by  a  rectangular  shaped  window  in  the  time  domain  can  be  viewed  as  a  convolution  of 
the  perfect  filter  response  with  a  sine  function  (the  Fourier  Transform  of  the  rectangular 
window)  in  the  frequency  domain. 

The  solution  is  to  use  a  non-rectangular  window  and  one  whose  Fourier 
Transform  has  a  narrower  main  lobe  and  smaller  sidelobes  than  the  sine  function.  The 
Hamming  window  is  one  that  is  commonly  used  [11].  Multiplying  the  coefficients  from 
(2.13)  by  this  window,  and  using  the  results  in  a  FIR  filter,  the  frequency  response 
needed  is  generated.  Energy  will  be  lost  at  the  filter  transitions,  which  is  primarily  the 
result  of  the  loss  of  orthogonality  from  truncating  the  filter.  For  detection,  instead  of  los¬ 
ing  the  energy  at  those  frequencies,  a  better  trade-off  would  be  a  small  amount  of  cross¬ 
correlation  between  the  filters  so  that  some  energy  appears  in  more  than  one  tile. 
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To  achieve  this  type  of  prototype  filter  the  impulse  response  can  be  modi¬ 
fied  to  have  a  passband  that  is  slightly  greater  than  tc/2.  Thus,  the  H  and  G  filters  are 
squeezed  together  slightly.  This  can  be  achieved  by  compressing  the  sine  envelope  of 
(2.13)  slightly.  At  the  same  time,  it  will  be  desirable  to  rescale  the  coefficients  slightly, 
so  the  sum  of  the  squares  equals  one.  These  modifications  to  (2.13)  give 


h(n )  = 


Is  . 

f  n+0.5^1 

—sine 

h 

l  c  J 

w{n) 


(2.14) 


where  -N  /  2  <  n  <  ( N  -  2)  /  2 ,  C  is  the  compression  variable,  5  is  the  scaling  variable,  N 
is  the  number  of  coefficients,  and  w(n)  is  the  Hamming  window  to  suppress  the  Gibb’s 
phenomena.  An  infinite  number  of  coefficients  and  no  window  could  be  used.  C= 2  and 
5=1  would  create  orthogonal  wavelets  filters.  For  these  filters,  the  greatest  cross  correla¬ 
tion  occurs  between  tiles  in  the  same  frequency  band  and  adjacent  in  time  when  N=5 1 2 
(the  number  of  coefficients),  with  values  C=l. 99375872328059,  5=1.00618488680080, 
and  a  Hamming  window  giving  a  cross  correlation  of  less  than  0.001.  This  is  called  the 
“modified  sine  filter”  [12],  [14]. 


B.  QUADRATURE  MIRROR  FILTER  BANK  (QMFB)  TREE 


This  section  describes  how  the  QMFB  tree  is  implemented  using  the  theory  in  the 
previous  section. 

1.  The  Receiver 

Orthogonal  wavelet  decomposition  of  the  unknown  signal  can  be  implemented  us¬ 
ing  QMFs  by  designing  filter  pairs  to  divide  the  input  signal  energy  into  two  orthogonal 
components  based  on  frequency.  The  tiles  are  used  to  refer  to  the  rectangular  region  of 
the  time-frequency  plane  containing  the  basis  function’s  energy.  By  arranging  the  QMF 
pairs  in  a  fully  developed  tree  structure  it  is  possible  to  decompose  the  wavefonn  in  such 
a  way  that  the  tiles  have  the  same  dimensions  within  each  layer  [4,  15].  Thus,  every  filter 
output  is  connected  to  a  QMF  pair  in  the  next  layer  as  shown  in  Figure  15. 
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The  architecture  of  the  QMFB  tree  used  in  this  thesis  is  illustrated  in  Figure  15.  Each 
QMF  pair  divides  a  digital  input  waveform  into  its  high  frequency  and  low  frequency 
components  with  a  transition  centered  at  n  I 2  .  A  nonnalized  input  of  one  sample/second 
is  assumed,  with  a  signal  bandwidth  of  [0,/r] .  Since  each  fdter  output  signal  has  half  the 

bandwidth,  only  half  the  samples  are  required  to  meet  the  Nyquist  criteria,  therefore  these 
sequences  are  then  down  sampled  by  two.  The  same  number  of  output  samples,  as  were 
input  are  returned.  For  example  if  100  samples  appear  at  the  input  of  the  first  QMF  pair, 
100  samples  appear  at  the  output.  Each  of  the  two  resulting  sequences  are  then  fed  into 
QMF  pairs  forming  the  next  layer  where  the  process  is  repeated,  and  so  on  down  the  tree. 


Figure  15  Quadrature  mirror  filter  bank  (QMFB)  Tree.  After  [4], 


Considering  the  square  of  each  element  of  the  input  waveform  represents  the 
waveform’s  energy  for  that  sample,  each  represents  the  energy  contained  in  the  corre¬ 
sponding  tile  in  the  left  most  time-frequency  diagram  shown  in  Figure  15.  Similarly,  the 
outputs  from  each  layer  of  the  tree  form  a  matrix  whose  elements,  when  squared  ap- 
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proximately  represent  the  energy  contained  in  the  tiles  of  the  corresponding  time- 
frequency  diagrams  shown  in  the  figure.  As  shown  in  Section  A. 5,  it  is  not  possible  to 
find  a  filter  that  perfectly  divides  the  energy  into  tiles  in  both  time  and  frequency  [4,  13]. 

When  the  waveform  consists  entirely  of  WGN,  the  tile’s  energy  will  have  random 
values  with  a  Chi-squared  probability  distribution.  When  a  detenninistic  signal  is  added, 
tiles  containing  energy  from  the  signal  will  have  probability  distribution  that  is  Chi- 
squared  with  non-centrality  parameters  and  will,  therefore,  tend  to  have  larger  mean  val¬ 
ues  and  thus  make  threshold  detection  a  possibility  [4,  13]. 

Since  the  transform  is  linear,  a  fundamental  limit  on  the  minimum  area  of  each  of 
the  tiles  exists.  However,  looking  at  Figure  15,  it  can  be  noted  that  each  layer  outputs  a 
matrix  of  energy  values  for  tiles  that  are  twice  as  long  (in  time)  and  half  as  tall  (in  fre¬ 
quency)  as  the  tile  in  the  previous  layer.  By  properly  comparing  these  matrices,  it  is  pos¬ 
sible  to  find  concentrations  of  energy  and  estimate  their  position  and  sizes  with  high 
resolution  in  both  time  and  frequency.  Using  these  techniques,  then,  a  waveform  can  be 
decomposed  and  the  bandwidths,  the  time  widths,  and  locations  in  the  time-frequency 
plane  can  be  estimated.  All  this  infonnation,  of  course,  can  then  be  used  by  the  intercep¬ 
tor  to  decide  how  many  transmitters,  and  which  types,  are  in  operation  [4,  13]. 

The  receiver  block  diagram  is  shown  in  Figure  3  (see  p.  6).  A  received  waveform 
is  bandpass  filtered  and  sampled  at  the  Nyquist  rate.  The  digital  sequence  is  then  fed  to 
the  QMFB  tree  where  it  is  decomposed.  Matrices  of  values  are  output  from  each  layer, 
and  are  then  squared  to  produce  numbers  representing  the  energy  in  each  tile. 

2.  Filters 

When  considering  only  FIR  filters,  some  requirements  are  necessary.  First,  the  fil¬ 
ters  must  meet  the  wavelet  requirements  described  in  Section  2.  Basically,  these  restrict 
the  possible  filter  coefficient  values  to  ensure: 

1)  An  orthogonal  decomposition  so  that  the  energy  in  sequences  output  from  each 
QMF  pair  will  equal  the  energy  input. 
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2)  That  the  output  from  the  H  filter,  as  labeled  in  Figure  15,  consists  of  low  fre¬ 
quency  components  of  the  input,  while  the  output  from  the  G  filter  consists  of  high  fre¬ 
quency  components. 

As  it  turns  out,  one  practical  consequence  of  these  requirements  is  that,  when  a 
suitable  H  filter  is  found,  the  G  filter  is  obtained  by  negating  and  time  reversing  every 
other  coefficient  value. 

Second,  the  filters  approximately  collect  energy  in  tiles.  Essentially,  they  must 
pass  as  much  energy  from  inside  a  tile  as  possible,  while  rejecting  as  much  as  possible 
from  outside  a  tile  with  a  reasonably  flat  pass  region. 

Some  filters,  such  as  the  Haar  filter,  meet  the  wavelet  requirements  that  perfectly 
tile  the  input  energy  in  time,  but  unfortunately,  does  not  tile  well  in  frequency.  The  oppo¬ 
site  of  the  Haar  filter,  in  this  respect,  would  be  the  sine  filter.  Both  were  described  in  Sec¬ 
tion  A.  5.  The  correct  filter  is  the  “modified  sine  filter”,  which  will  return  a  good  tile  in 
time  and  frequency  [4,  13]. 


C.  SIMULATION  PROGRAMS 


Some  of  the  programs  used  in  the  simulations  for  this  thesis  are  described  in  this 
section.  The  simulations  are  used  to  verify  mathematically  derived  results  and  to  provide 
results  when  the  mathematics  becomes  extremely  difficult.  All  these  simulations  were 
done  on  a  Pentium  4-based  personal  computer,  with  a  2.0  GHz  CPU  and  1  Gb  of  RAM 
using  MATLAB ®  6. 1  for  Windows  2000. 

The  organization  of  the  simulation  programs  follows  the  receiver  block  diagram 
in  Figure  3.  The  wavefonn  is  assumed  to  be  properly  bandpass  filtered  and  sampled  at  the 
Nyquist  rate.  The  MATLAB®  code  for  generating  the  sampled  input  wavefoms  (BPSK, 
Frank  Code,  Polyphase,  FMCW,  and  Costas  signals)  can  be  found  in  Reference  [17].  The 
input  signals  are  two-row  data  containing  the  first  one  called  I,  which  is  the  real  part  of 
the  signal,  and  the  second  is  called  Q  and  is  the  imaginary  part  of  the  signal  generated. 
The  signals  generated  are  used  as  input  into  the  QMFB  program,  which  is  described  next 
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in  this  section  and  carries  out  the  function  of  the  QMFB  Tree.  The  output  is  then  written 
to  the  hard  drive.  The  QMFB  tree  output  data,  from  the  hard  drive,  can  then  be  analyzed 
using  code  written  to  execute  the  feature  extraction  from  every  processed  waveform. 

1.  Quadrature  Mirror  Fiter  Bank  Tree  Programs  and  Procedures 

Appendix  A  contains  the  m  file  startpoint.m  that  loads  the  signal,  which  will  be 
filtered  by  the  QMFB  (all  programs  and  program  variables  are  shown  in  bold).  The  code 
requires  the  sampling  frequency  fs  and,  in  case  the  signal  length  does  not  correspond  to  a 
power  of  two  in  the  number  of  samples,  the  signal  will  be  padded  with  zeros  until  the 
length  reaches  the  next  power  of  2.  In  this  manner  the  amount  of  layers  resulting  from  the 
processing  has  direct  relationship  with  the  length  of  the  signal  (number  of  samples)  in  the 
sense  that  the  number  of  layers  is  detennined  by  the  power  of  2  that  the  length  of  the  sig¬ 
nal  has.  For  example,  if  a  signal  row  data  has  a  length  of  2800  samples,  this  signal  is 
padded  with  zeros  until  reach  the  next  power  of  2,  which  is  212  =  4096  samples,  indicat¬ 
ing  that  the  number  of  layers  resulting  from  the  QMFB  processing  is  12  layers  (lay=  12). 
Afterwards,  the  code  filters  the  signal  applying  the  qmfb.m  function. 

Appendix  A  also  contains  the  m  file  function  qmfb.m  that  decomposes  the  input 
sequence  as  described  in  Section  A  (Background).  The  function  takes  the  incoming  signal 
f  formally  the  signal  tt  from  the  startpoint.m,  and  processes  it  with  a  QMFB  tree  struc¬ 
ture,  where  filter  specifies  the  file  containing  the  filter  coefficients  of  the  modified  sine 
filter  to  be  applied  through  the  QMFB  tree  over  the  signal.  Output  sequences  from  each 
layer  of  the  QMFB  are  written  to  files  in  the  same  directory  where  the  function  resides  on 
the  hard  drive  where  variable  lay  is  replaced  by  the  layer  number.  These  are  ASCII  files 
and  the  data  is  stored  in  a  two-matrix  fonnat.  The  first  matrix  is  called  R  and  contains  the 
real  part  of  the  data  and  the  second  matrix  is  called  Q  and  contains  the  imaginary  part  of 
the  data  in  each  column  representing  the  output  from  a  particular  filter  in  the  layer.  Fre¬ 
quency  is  represented  across  each  line  in  the  file  (lowest  frequency  to  the  left)  and  time  is 
represented  down  the  file  (lower  values  representing  later  time). 

The  number  of  layers  is  then  determined  from  the  length  of  the  input  sequence. 
The  input  matrix  of  each  layer  I  is  initially  set  equal  to  the  input  sequence  and  the  output 
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matrix  for  each  layer,  out,  is  set  equal  to  I  with  the  intention  to  save  memory.  The  code 
then  loops  once  for  each  layer  of  the  QMFB.  Inside  the  loop,  the  layer  number  is  first 
displayed  on  the  screen  to  indicate  to  the  user  what  the  program  is  doing  and  flag  is  ini¬ 
tialized  to  1 .  Then,  out  is  reshaped  to  match  the  dimensions  required  for  the  output  of  the 
layer. 

The  code  then  loops  again  inside  the  first  loop  once  for  each  column  of  the  cur¬ 
rent  input  matrix.  G  and  H  are  the  output  sequences  from  each  filter  pair.  As  described  in 
Sections  3  and  4,  these  sequences  represent  a  decomposition  of  the  input  into  lower  and 
upper  frequency  bands,  which  alternates  down  the  layer.  This  is  tracked  by  a  flag  and  the 
sequences  are  written  to  out  in  the  correct  order.  At  the  end,  the  output  for  the  current 
layer,  which  is  a  matrix,  is  used  as  the  input  for  the  next,  and  the  output  data  is  written  to 
disk. 

Appendix  A  additionally  contains  the  m  file  tsinc.m  that  uses  N  coefficients  gen¬ 
erated  by  the  m  file  tsinc_su.m  also  listed  in  Appendix  A.  This  file  generates  the  coeffi¬ 
cients  described  in  Equation  (2.14)  under  the  h.clat  file  (low  pass  filter)  and  g.dat  file 
(high  pass  filter)  for  a  modified  sine  filter  with  a  Hamming  window.  These  are  ASCII 
files  and  the  data  is  stored  in  double  precision  format  as  a  matrix.  Once  generated,  the 
coefficients  are  saved  to  the  hard  drive  and  called  upon  as  needed  by  tsinc.m,  which  uses 
them  to  decompose  the  column  vector  cO  into  a  column  vector  ell  (high  frequency)  and 
into  a  column  vector  cl  (low  frequency)  using  a  built-in  MATLAB8  function  called  filter 
which  is  designed  to  work  like  a  FIR  filter  and  i  is  used  to  accomplish  the  decimation  by 
2. 

The  startpoint.m  file,  contained  in  Appendix  A,  automatically  displays  the  time- 
frequency  plane  from  4  particular  layers  through  the  QMFB  processing.  The  first  figure 
displays  4  contour  plots  of  the  selected  layers.  The  second  figure  displays  the  same  4  lay¬ 
ers  using  the  pcolor  MATLAB®  command  (4  contour  plot  in  a  figure  and  4  color  plots  in 
other).  The  third  figure  displays  the  same  4  layers  in  a  3  dimensions  view  using  the  mesh 
MATLAB8  command.  The  fourth  and  last  figure  displays  the  4  layers  in  a  3  dimensional 
view  using  the  surf  MATLAB 8  command.  That  gives  a  total  of  four  figures  each  with 
four  plots. 
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In  the  case  of  all  output  matrix  layers  need  to  be  plotted,  four  more  files  are  pro¬ 
vided  in  Appendix  A.  The  multiple _colorplot.ni  file  displays  all  time-frequency  plane 
from  all  the  layers  using  the  pcolor  MATLAB®  command.  The  multiple _contourplot.m 
file  displays  all  time-frequency  plane  from  all  the  layers  using  the  contour  MATLAB1' 
command.  The  multiple _meshplot.m  file  displays  a  3-dimensional  plot  from  all  the  layers 
using  the  mesh  MATLAB8  command.  The  multiple _surfplot.m  file  displays  a  3- 
dimensional  plot  from  all  the  layers  using  the  surf  MATLAB8’  command. 


This  chapter  presented  the  QMFB  theory,  discussing  the  mathematical  back¬ 
ground  for  signal  detection  and  the  waveform  decomposition  using  Wavelets,  and  fin¬ 
ished  with  a  description  of  the  MATLAB 8  codes  used  to  implement  the  QMFB  tree.  The 
next  chapter  has  the  most  relevant  work  of  this  thesis.  There,  a  set  of  different  and  well 
known  LPI  radar  signals  are  processed  with  the  QMFB  tree  and  the  output  are  displayed 
and  analyzed  as  a  feature  extraction  procedure.  Chapter  III  finishes  with  an  analysis  and 
comparison  of  the  different  signals  processing  results,  where  the  performances  of  the 
QMFB  three  are  specified  for  every  class  of  signal  depending  on  the  results  obtained  in 
the  corresponding  signal  processing. 
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III.  OVERVIEW  OF  LPI  EMITTER  WAVEFORMS  AND 
PROCESSING  BY  THE  QMFB  TREE 


This  chapter  discusses  the  results  of  processing  the  LPI  signals  from  the  LPIG 
code  as  described  in  [17].  Chapter  three  finishes  with  an  analysis  and  comparison  of  the 
different  signals  processing  results,  where  the  perfonnances  of  the  QMFB  tree  are  speci¬ 
fied  for  every  class  of  signal  depending  on  the  results  obtained  in  the  corresponding  sig¬ 
nal  processing. 


In  this  chapter  the  signals  are  briefly  described.  Next,  the  parameters  from  each 
signal  are  extracted  from  the  QMFB  results.  To  standardize  the  manner  in  which  the  sig¬ 
nals  are  named,  Table  1  is  used.  This  table  applies  abbreviated  names  to  the  signal  with 
numbers  to  identify  the  frequencies  involved  and  some  of  the  parameters  such  as  the 
number  of  phases  and  code  periods  per  phase.  Table  1  shows  an  example  of  the  LPIG 
naming  convention  for  the  collection  of  signals.  The  example  shows  a  BPSK  signal  with 
a  2  kHz  as  carrier  frequency,  4  kHz  as  sampling  frequency,  13  Bits  Barker  Code,  and  1 
cycle  per  Barker  Bit.  The  list  of  signals  generated  by  LPIG  and  processed  by  the  QMFB 
tree  appears  in  Appendix  B. 
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B  2  4  13  1  s.mat 

t  t  t  t  t  t _ 


s  : 


SNR: 


0  : 
-6: 


signal  only 
OdB 
-6  dB 


Number  of  Cycles  per  Barker  Bit  (appears  only  in  BPSK) 

Bit  for  Barker  code  for  BPSK 

Modulation  bandwidth  for  FMCW 

Numbers  of  phase  code  for  Frank  and  Polyphase 

Sequence  for  Costas,  FSK/PSK  Costas  and  FSK/PSK  Target 

1  or  2  frequencies  for  T 

Sampling  frequency  (KHz) 


Carrier  frequency  (KHz) 


Code: 


B:  BPSK 
C:  Costas 
F:  FMCW 
FR:  Frank  Code 

FSK  PSK  C:  FSK/PSK  Costas 
FSKPSKT:  FSKPSK  Target 

P#:  Polyphase  Code  (#:1,2,3,4) 

T:  Test  Signal 


Table  1  Standard  named  signals  from  LPIG. 


A.  TONE  TEST 


To  probe  the  type  of  results  generated  by  the  QMFB  tree,  a  tone  test  was  applied 
using  a  single  frequency  input.  The  output  matrices  from  the  various  layers  were  exam¬ 
ined  to  verily  that  the  code  qmfb.m  is  correct,  and,  second,  to  provide  the  frequency  re¬ 
sponse  information  from  the  modified  sine  filter. 

1.  T_l_7_l_s 

The  first  of  these  tests  signals  uses  1  kHz  tone  sampled  at  7  kHz  with  amplitude 
of  unity  and  zero  phases  (T_l_7_l_s).  The  results  are  presented  in  Figure  16  showing 
five  of  the  seven  layers  obtained  by  the  colorplot.m  file.  Since  the  first  and  the  last  layers 
in  the  processing  always  are  going  to  be  a  row  data  (vector,  not  matrix)  they  can  not  be 
displayed  in  a  time-frequency  diagram. 
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T  1  7  1  s 


d).  e). 


Figure  16  Time-frequency  plots  for  T_l_7_l_s,  from:  a)  layer  2,  b)  layer  3,  c)  layer  4,  d) 
layer  5,  e)  layer  6. 


In  this  figure,  it  is  possible  to  see  how  each  layer  results  in  a  matrix  of  energy  values  for 
tiles  that  are  twice  as  long  (in  time)  and  half  as  tall  (in  frequency)  as  the  tile  in  the  previ¬ 
ous  layer.  By  properly  comparing  these  matrices,  it  is  possible  to  find  concentrations  of 
energy  and  to  estimate  their  position  and  sizes  with  high  resolution  in  both  time  and  fre¬ 
quency.  Figure  17  shows  the  output  matrix  from  layer  2.  The  total  number  of  layer  in  this 
example  is  7  ( lay=l )  This  results  in  very  narrow  tiles  in  time,  but  very  wide  tiles  in  fre¬ 
quency  and  indicates  that  the  resolution  of  layer  2  is  better  in  time  than  in  frequency.  The 
tiles  have  a  resolution  in  frequency  determined  by 


A/  = 


fs 


2(number  of  columns  of  the  output  matrix  at  this  layer)  -  1 

/,  _ L _ 

2(2number  of  ,his  layer  _ , )  2(number  of  tiles  in  frequency) 

7000 


(3.1) 


2(3) 


=  1166.67  Hz 
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and  since  the  sampling  period  of  the  signal  T  is  defined  as  1  /  fs ,  the  resolution  in  time  is 
determined  by 


^  _  (number  of  samples  of  input  signal  padded  with  zcros)(7) 

(number  of  rows  of  the  output  matrix  at  this  layer)- 1 
(number  of  samples  of  input  signal  padded  with  zeros)(  Vf  ) 
(number  of  tiles  in  time) 

(number  of  samples  of  input  signal  padded  with  zeros)(  Vf  ) 

(total  number  of  layers  -  number  of  actual  layer)  j\ 

=  (128X^°°°)  =  0.00059  =  590  jus. 
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(3.2) 


Thus,  from  Figure  17,  the  signal  is  contained  in  the  frequency  band  1166.67  Hz  to 
2333.34  Hz,  and  from  0  ms  to  14.75  ms. 


Figure  17  Time-frequency  plot  from  layer  2  for  T_l_7_l_s. 


Figure  18  shows  the  output  matrix  from  layer  4.  The  tiles  have  a  resolution  in  fre¬ 
quency  of  233.333  Hz  and  a  resolution  in  time  of  2.612  ms.  This  figure  shows  the  signal 
contained  in  the  frequency  band  9333.33  Hz  to  1 1666.7  Hz,  and  from  0  to  15.672  ms. 
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Figure  18  Time-frequency  plot  from  layer  4  for  T_l_7_l_s. 


Figure  19  shows  the  output  matrix  from  layer  6  with  tiles  being  very  defined  in 
frequency.  The  tiles  have  a  resolution  in  frequency  of  55.5556  Hz  and  a  resolution  in 
time  of  18.286  ms.  Therefore,  this  figure  shows  the  signal  contained  in  the  frequency 
band  944.445  Hz  to  1055.56  Hz  and  from  0  to  18.286  ms.  This  effectively  indicates  that 
at  a  higher  layer,  the  frequency  resolution  is  higher.  Thus,  it  could  be  seen  that  the  central 
frequency  corresponds  to  1000  Hz  approximately.  It  is  important  to  note  that  the  time  in¬ 
formation  is  not  missing.  The  only  difference  is  that  if  some  more  detailed  specific  time 
information  is  needed,  it  could  be  obtained  from  previous  layers  other  than  layer  6. 
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figure  19  Time-frequency  plot  from  layer  6  for  T_l_7_l_s. 


It  is  also  instructive  to  look  at  the  frequency  through  time.  Using  the  surf  plot. m 
file  in  Appendix  A,  the  frequency-energy  plots  of  the  tone  test  are  shown  in  Figure  20  to 
Figure  22.  Figure  20  and  Figure  21  show  a  surf  plot  (view  from  the  frequency  axis)  from 
layer  2  and  4,  respectively.  In  Figure  20,  the  frequency  resolution  is  1166.67  Hz,  but 
clearly  the  energy  concentrates  close  to  1000  Hz.  In  Figure  21,  the  frequency  resolution  is 
233.333  Hz,  been  higher  than  the  previous  case. 
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Figure  20  Frequency-energy  plot  from  layer  2  for  T_l_7_l_s. 


36 


Figure  2 1  Frequency-energy  plot  from  layer  4  for  T_l_7_l_s. 


Figure  22  shows  a  surf  plot  (view  from  the  frequency  axis)  from  layer  6.  The  fre¬ 
quency  resolution  is  55.5556  Hz,  been  higher  than  the  previous  cases  in  layers  2  and  4, 
which  corroborates  that  the  central  frequency  is  1000  Hz  for  the  test  signal.  The  last  three 
figures  reveal  that,  the  higher  the  layer,  the  higher  the  resolution  in  frequency.  Layer  6 
gives  the  best  resolution  and  is  shown  as  a  3 -dimensional  surf  plot  in  Figure  23  with  a 
frequency  resolution  of  55.5556  Hz. 
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Figure  23  Surf  plot  from  layer  6  for  T_l_7_l_s. 

Table  2  presents  a  summary  of  the  signal  processing,  showing  the  relevant  pa¬ 
rameters  extracted  from  the  QMFB  layers  2,  4  and  6. 


T_l_7_l_s 

Frequency 

Resolution 

Time  Resolution 

Comment 

Layer  2 

1166.67  Hz 

0.000590  ms 

Layer  4 

233.333  Hz 

2.612  ms 

Layer  6 

55.5556 

18.286  ms 

Carrier  Frequency 

1  kHz 

- 

Sampling  Fre¬ 

quency 

7  kHz 

SNR 

Only  Signal 

- 

Table  2  Signal  processing  summary  for  T_l_7_l_s. 


2.  T_12_7_2_s 

A  second  signal  test  is  applied  to  the  QMFB  tree.  The  signal  was  generated  with 
two  tones,  1  kHz  and  2  kHz,  sampled  at  7  kHz  with  amplitude  of  unity  and  zero  phase 
(T_12_7_2_s)  are  used  for  this  second  test.  The  amount  of  layers  obtained  when  process¬ 
ing  the  signal  with  the  QMFB,  as  it  was  seen  in  Chapter  II  Section  C.l,  depends  on  the 
length  of  the  signal  to  process,  and  corresponds  to  the  power  of  2  that  determine  this 
length  after  padding  with  zeros  if  it  is  necessary.  Since  the  first  and  the  last  layers  in  the 
processing  always  are  going  to  be  a  row  data  (vector,  not  matrix)  they  can  not  be  dis¬ 
played  in  a  time-frequency  diagram.  Thus,  the  rest  of  the  layers  obtained  are  presented  in 
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Figure  24.  In  this  figure,  it  is  possible  to  see  how  each  layer  outputs  a  matrix  of  energy 
values  for  tiles  with  an  increment  in  frequency  resolution  as  the  layers  increase.  By  prop¬ 
erly  comparing  these  matrices,  it  is  possible  to  find  concentrations  of  energy,  and  esti¬ 
mate  their  position  and  sizes  with  high  resolution  in  both  time  and  frequency. 


Figure  24  Time-frequency  plots  for  T_12_7_2_s,  from:  a)  layer  2,  b)  layer  3,  c)  layer  4, 
d)  layer  5,  e)  layer  6. 


Figure  25  shows  the  output  matrix  from  layer  2.  The  resolution  in  time  is  0.59  ms 
being  better  than  the  resolution  in  frequency  that  is  1 166.67  Hz.  In  fact,  the  signal  covers 
a  wide  band  in  frequency. 
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Contour  plot  from  layer  2  for  T_12_7_2_s. 


Figure  26  shows  the  output  matrix  from  layer  4.  The  resolution  in  frequency  is 
233.333  Hz  being  higher  than  in  layer  2.  The  plot  demonstrates  that  the  carrier  frequency 
is  close  to  the  test  signal,  which  are  1000  Hz  and  2000  Hz,  but  still  with  some  bandwidth, 
making  it  not  as  precise  as  desired.  However,  Figure  27  shows  the  output  matrix  from 
layer  6  with  frequency  resolution  of  55.5556  Hz  and  there  the  central  frequencies  corre¬ 
sponds  to  1000  Hz  and  2000  Hz. 
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Figure  26  Contour  plot  from  layer  4  for  T_12_7_2_s. 
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Figure  27  Contour  plot  from  layer  6  for  T_12_7_2_s. 

As  with  the  previous  test  signal,  the  frequency  through  time  is  considered.  Using 
the  surfplot.m  file,  a  frequency-energy  plot  of  the  tone  test  signal  is  shown  in  Figure  28  to 
Figure  30.  Figure  28  and  Figure  29  show  a  surf  plot  (view  from  the  frequency  axis)  from 
layer  2  and  4,  respectively.  In  Figure  28,  the  frequency  resolution  is  1 166.67  Hz  (not  high 
enough),  but  clearly  the  energy  concentrates  between  1000  Hz  and  2000  Hz.  In  Figure 
29,  the  frequency  resolution  is  233.333  Hz  (higher  than  the  previous  case),  approximating 
the  both  frequencies  at  1000  Hz  and  2000  Hz. 
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Figure  28  Frequency-energy  plot  from  layer  2  for  T_12_7_2_s. 

41 


Figure  29  Frequency-energy  plot  from  layer  4  for  T_12_7_2_s. 


Figure  30  shows  a  surf  plot  (view  from  the  frequency  axis)  from  layer  6.  The  fre¬ 
quency  resolution  is  55.5556  Hz  and  verifies  that  the  central  frequencies  are  1000  Hz  and 
2000  Hz  for  the  test  signal. 


The  last  three  figures  demonstrate  that  the  higher  the  layer,  the  higher  the  resolu¬ 
tion  in  frequency.  Thus,  the  best  resolution  is  that  presented  by  layer  6.  Figure  3 1  shows  a 
3 -dimensional  surf  plot  from  layer  6  with  a  frequency  resolution  of  55.5556  Hz. 
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Figure  3 1  Surf  plot  from  layer  6  for  T_  12_7_2_s. 


Table  3  presents  a  summary  of  the  signal  processing,  showing  the  relevant  pa¬ 
rameters  extracted  from  the  QMFB  layers  2,  4  and  6. 


T_l_7_l_s 

Frequency 

Resolution 

Time  Resolution 

Comment 

Layer  2 

1166.67  Hz 

0.000590  ms 

Layer  4 

233.333  Hz 

2.612  ms 

Layer  6 

55.5556 

18.286  ms 

Carrier  Frequency 

1  and  2  kHz 

- 

Sampling  Fre¬ 

quency 

7  kHz 

" 

SNR 

Only  Signal 

- 

Table  3  Signal  processing  summary  for  T_12_7_2_s. 


Up  to  this  point,  two  different  methods  were  used  to  show  the  output  matrix.  They 
are  the  color  and  contour  plots.  From  both  it  is  possible  to  conduct  the  feature  extraction, 
but  the  selected  plot  to  be  shown  in  the  following  signal  processing  will  depend  on  which 
will  display  a  better  resolution  from  any  layer. 
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B.  BINARY  PHASE  SHIFT  KEYING  (BPSK) 


1.  Brief  Description 

Phase  shift  keying  (PSK)  is  a  form  of  angle-modulated,  constant-amplitude  digital 
modulation.  PSK  is  similar  to  conventional  phase  modulation  except  that  with  PSK  the 
input  signal  is  a  binary  digital  signal  and  only  a  limited  number  of  output  phases  are  pos¬ 
sible.  BPSK  is  not  a  technique  employed  in  radar  modulation,  but  is  a  well-known  LPI 
signal  used  in  communication,  thus  constitutes  itself  as  an  excellent  test  signal  in  evaluat¬ 
ing  the  performance  of  the  QMFB  processing. 

With  binary  phase  shift  keying  (BPSK),  two  output  phases  are  possible  for  a  sin¬ 
gle  carrier  frequency  (binary,  meaning  2).  One  output  phase  represents  logic  1  and  the 
other  a  logic  0.  As  the  input  digital  signal  changes  state,  the  phase  of  the  output  carrier 
shifts  between  two  angles  that  are  180°  out  of  phase.  BPSK  is  a  fonn  of  suppressed- 
carrier,  square -wave  modulation  of  a  continuous  wave  (CW)  signal  [14]. 

In  a  BPSK  transmitter,  the  modulator  receives  two  input  signals,  the  carrier  and 
the  binary  data.  If  +  1V  is  assigned  to  a  logic  1  and  -IV  is  assigned  to  a  logic  0,  the  input 
carrier  (sin  coct )  is  multiplied  by  either  a  +  or  -1.  Consequently,  the  output  signal  is  either 

+  1  sin  a>ct  or  -1  sin  corl .  The  first  represents  a  signal  that  is  in  phase  with  the  reference  os¬ 
cillator  while  the  latter  is  a  signal  180°  out  of  phase  with  the  reference  oscillator.  Each 
time  the  input  logic  condition  changes,  the  output  phase  changes.  Mathematically,  the 
output  of  a  BPSK  modulator  is  proportional  to 

BPSK  output  =  [sin(2;r//)][sin(2;r/cr)]  (3.3) 

where  fa  =  maximum  fundamental  frequency  of  binary  input  (Hertz)  and  fc=  reference 

carrier  frequency  (Hertz).  Solving  for  the  trig  identity  for  the  product  of  two  sine  func¬ 
tions 

BPSK  output  =  ^cos[2x(fc-fa)t]-^cos[2x(fc+fa)t].  (3.4) 
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Figure  32  below  shows  a  basic  block  diagram  of  the  transmitter  design.  The  signal 
x{t )  is  a  sinusoidal  CW.  After  sampling  at  the  Nyquist  rate,  the  modulate  signal  is  cre¬ 
ated  by  adding  a  n  -bit  Barker  code.  This  code  has  been  widely  used  because  of  its  low- 
side  lobes  at  zero  Doppler.  Once  the  signal  has  been  modulated,  white  Gaussian  noise  is 
added  at  the  desired  level. 


Figure  32  BPSK  transmitter  block  diagram. 


The  modulating  signal  waveform  is  represented  by  the  dashed  lines  in  the  upper 
graph  in  Figure  33.  The  number  of  periods  per  Barker  Bit  is  set  to  one,  signifying  that 
one  full  period  of  the  sampled  signal  fits  within  one  bit  of  the  13-bit  Barker  code.  The 
first  five  bits  of  the  Barker  code  are  +1  and  the  next  2  bits  are  -1,  so  five  full  periods  un¬ 
der  the  first  +1  portion  of  the  modulating  wavefonn,  2  full  periods  under  the  -1  portion 
of  the  modulating  wavefonn,  and  so  forth  are  seen.  In  this  figure,  only  the  first  12  bits  of 
the  13-bit  Barker  code  are  represented.  Table  4  shows  Barker  codes  for  7,  1 1,  and  13  bits. 
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Figure  33  Sampled  signal  (upper  plot  in  blue  color),  modulating  signal  (upper  plot,  in 
red  color)  and  modulated  signal  (lower  plot,  in  blue  color). 


Number  of  bits 

Barker  Code 

7 

+  +  + - +  - 

11 

+  +  + - +  - 

13 

+++++--++-+-+ 

Table  4  Barker  Code  for  7,  1 1  and  13  Bits. 


2.  Processing  BPSK  Signals  with  QMFB  Tree 

The  BPSK  signals  to  be  worked  with  are  given  in  Table  5.  All  the  signals  were 
generated  with  a  carrier  frequency  (/c)  of  1000  Hz  and  sampling  frequency  (fs)  of  7000 

Hz.  The  number  of  bits  per  Barker  code,  the  number  of  cycles  per  bit  ( cpb ),  and  the  SNR 
are  variable  parameters. 
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BPSK 

Number  of  bits 
per  Barker  code 
(Bits) 

Number  of  cy¬ 
cles  per  bit  (cpb) 

SNR 

1 

B  1  7  7  1  s 

7 

1 

Signal  Only 

2 

B  1  7  7  1  0 

7 

1 

OdB 

3 

B  1  7  7  1  -6 

7 

1 

-6  dB 

4 

B  1  7  11  1  s 

11 

1 

Signal  Only 

5 

B  1  7  11  1  0 

11 

1 

OdB 

6 

B  1  7  11  1  -6 

11 

1 

-6  dB 

7 

B  1  7  7  5  s 

7 

5 

Signal  Only 

8 

B  1  7  7  5  0 

7 

5 

OdB 

9 

B  1  7  7  5  -6 

7 

5 

-6  dB 

10 

B  1  7  11  5  s 

11 

5 

Signal  Only 

11 

B  1  7  11  5  0 

11 

5 

OdB 

12 

B  1  7  11  5  -6 

11 

5 

-6  dB 

Table  5  BPSK  signals  to  be  processed  by  QMFB  Tree. 


From  the  list  of  signals  already  processed  by  the  QMFB  (in  Table  5),  only  one  of 
them  is  shown  next.  The  rest  of  the  signals  can  be  seen  in  a  technical  report  that  will  be 
published  soon  [21]. 


a.  B_l_7_7_5_s 

This  BPSK  signal  was  generated  with  the  parameters  described  in  Table  6. 
The  code  period  of  the  BPSK  signal  is 
T  (Number  of  cycles  per  Barker  bit)(Number  of  bits  per  Barker  code) 

J-  Q 

Carrier  Frequency 

_  (cpb)(Bits)  _  (5X7)  _  0  035  s  =  35  ms  <3  5) 

c  fc  1000 

The  bandwidth  of  the  signal  depends  on  the  cycles  per  bit  (or  chirp),  given 
by  the  LPIG  code  in  the  generation  of  the  signal.  Thus, 

BW  =  ^  = - 1000  Hz - =  200  Hz  =  0.2  KHz.  (3.6) 

cpb  5  cycles  per  bit 

Now,  and  according  to  the  description  of  the  code  given  in  Section  C.l  in 
Chapter  III,  the  signal  is  processed  as  an  input  to  the  QMFB  tree  to  get  the  different  out¬ 
put  layers.  In  this  manner,  Figure  34  was  created  which  shows  the  output  matrix  from 
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layer  6  in  a  contourplot  with  a  frequency  resolution  of  55.55  Hz  and  a  time  resolution  of 
9.677  ms.  The  values  of  the  carrier  frequency,  bandwidth,  and  code  period  are  extracted. 
Figure  35  shows  the  output  matrix  from  layer  2,  there  it  is  noted  the  phase  shifts  in  the 
signal,  the  Barker  code  sequence,  and  corroborate  the  code  period. 
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Figure  34  Output  matrix  from  layer  6  of  B_l_7_7_5_s  (contourplot). 
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Figure  35  Output  matrix  from  layer  2  of  B_l_7_7_5_s  (contourplot). 


Table  6  shows  a  summary  of  the  signal  processing  demonstrating  a  com¬ 
parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


B  1  7  7  5  s 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1  kHz 

Sampling  Frequency 

7  kHz 

7  kHz 

given 

Bits 

7  bits 

7 

cpb 

5 

- 

Bandwidth 

0.2  kHz 

0.222  kHz 

0.022  kHz  error 

Code  Period 

35  ms 

38.708  ms 

3.708  ms  error 

SNR 

Only  Signal 

- 

Table  6  Signal  processing  summary  for  B_l_7_7_5_s. 


b.  B_l_7_7_5_0 

This  BPSK  signal  was  generated  with  the  parameters  described  in  Table  7. 
The  different  output  layers  were  obtained  conducting  the  signal  processing  with  the 
QMFB  tree.  Figure  36  presents  the  output  matrix  from  layer  6  with  a  frequency  resolu¬ 
tion  of  55.55  Hz  and  a  time  resolution  of  9.677  ms.  There  the  values  of  the  carrier  fre¬ 
quency  and  the  bandwidth  are  extracted  matching  the  signal  values.  Since  it  still  is 
possible  to  find  a  pattern  or  a  frequency  in  the  repetition  of  tiles  and  the  resolution  in  time 
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is  acceptable,  it  was  not  necessary  to  examine  in  previous  layers  to  find  the  code  period. 
Consequently,  this  parameter  was  extracted  from  the  same  figure,  matching  the  value  in 
(3.5). 


‘Tins  -  Frequency  Plene  LsyerS" 


Figure  36  Output  matrix  from  layer  6  of  B_l_7_7_5_0  (contourplot). 


Table  7  presents  a  summary  of  the  signal  processing  demonstrating  a 
comparison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB 
tree. 


B  1  7  7  5  0 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1  kHz 

Sampling  Frequency 

7  kHz 

7  kHz 

given 

Bits 

7  bits 

- 

cpb 

5 

- 

Bandwidth 

0.2  kHz 

0.333  kHz 

0.133  kHz  error 

Code  Period 

35  ms 

38.708  ms 

3.708  ms  error 

SNR 

OdB 

- 

Table  7  Signal  processing  summary  for  B_l_7_7_5_0. 
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c. 


B_l_7_7_5_-6 


This  BPSK  signal  was  generated  with  the  parameters  described  in  Table  8. 
Processing  the  signal  with  the  QMFB  Tree  to  find  the  different  output  layers,  Figure  37 
was  obtained,  which  show  the  output  matrix  from  layer  6  with  a  frequency  resolution  of 
55.55  Hz  and  a  time  resolution  of  9.677  ms.  There  the  values  of  the  carrier  frequency  and 
the  bandwidth  are  extracted.  The  frequency  matches  the  signal  value,  but  the  bandwidth 
gives  a  higher  value.  It  was  not  possible  to  extract  the  code  period  from  the  plot  due  to 
the  noise  in  the  signal. 


‘Time  -  F<ieqtJBnty  Plane  LsysrS" 


Figure  37  Output  matrix  from  layer  6  of  B_l_7_7_5_-6  (contourplot). 


Table  8  presents  a  summary  of  the  signal  processing  indicating  a  compari¬ 
son  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 
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B  1  7  7  5  -6 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1  kHz 

Sampling  Frequency 

7  kHz 

7  kHz 

given 

Bits 

7  bits 

- 

cpb 

5 

- 

Bandwidth 

0.2  kHz 

0.444  kHz 

0.244  kHz  error 

Code  Period 

35  ms 

- 

SNR 

-6  dB 

- 

Table  8  Signal  processing  summary  for  B_l_7_7_5_-6. 


C.  FREQUENCY  MODULATED  CONTINUOUS  WAVE  (FMCW) 

1.  Brief  Description 

Linear  frequency  modulation  in  addition  to  being  an  LPI  radar  technique  has  the 
added  advantage  that  the  modulation  is  readily  compatible  with  solid-state  transmitters. 
The  most  popular  modulation  is  the  triangular  modulation  of  a  Frequency  Modulated 
Continuous  Wave  (FMCW).  The  linear  FMCW  emitter  uses  a  continuous  100%  duty- 
cycle  waveform  so  that  both  the  target  range  and  the  Doppler  infonnation  can  be  meas¬ 
ured  unambiguously  while  maintaining  a  low  probability  of  intercept.  The  FMCW  wave¬ 
form  represents  the  best  use  of  the  output  power  available  from  solid-state  devices. 
Linear  FMCW  is  easier  to  implement  than  phase  code  modulation  as  long  as  there  is  not 
strict  demand  on  linearity  over  the  modulation  bandwidth  [1]. 

The  linear  frequency  modulated  triangular  wavefonn  and  the  Doppler  shifted  sig¬ 
nal  are  shown  in  Figure  38.  The  triangular  modulation  consists  of  two  linear  frequency 
modulation  sections  with  positive  and  negative  slopes.  With  this  configuration,  the  range 
and  Doppler  frequency  of  the  detected  target  can  be  extracted  unambiguously  by  taking 
the  sum  and  the  difference  of  the  two  beat  frequencies  [1]. 
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Figure  38  Linear  frequency  modulated  triangular  wavefonn  and  the  Doppler 
shifted  signal. 


The  frequency  of  the  transmitted  signal  for  the  first  section  is 

,  if  Af 

/i(O  =  /0 — — +  — ■ t  (3.7) 

bn 

for  0  <  t  <  tm  and  zero  elsewhere.  Here  /0  is  the  RF  carrier,  A F  is  the  transmitted  modu¬ 
lation  bandwidth,  and  tm  is  the  modulation  period.  The  rate  of  the  frequency  change  or 
chirp  rate  F  is 


F  = 


A  F 


t 


m 


(3.8) 


The  phase  of  the  transmitted  RF  signal  is 


=  dx. 

0 


(3.9) 


Assuming  that  (/)0  —  0  at  t  =  0 , 
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<j)x  ( t )  =  2  n 


r 


AF 


f - 

v  2  j 


AF  2 
t  + — t 
2t„. 


(3.10) 


for  0  <  t  <  t  .  The  transmit  signal  is  given  by 


.s’,  (t)  =  a  sin  In 


r  A F" 

JO  rs 

v  ^  y 


AF  2 
t  + — t 
2 1 


(3.11) 


The  frequency  of  the  transmitted  waveform  of  the  second  section  is 

,  ,  ,  ,  AF  AF 

/2(0  =  /o+— - - < 


(3.12) 


for  0  <  t  <  t  .  The  transmitted  baseband  signal  is  given  by 


s2  ( t )  =  an  sin  2 n 


f  AF3  AF  2 

/0+ —  r - r 

2  J  2t„, 


(3.13) 


2.  Processing  FMCW  Signals  with  QMFB  Tree 

The  FMCW  signals  to  be  worked  with  are  given  in  Table  9.  All  the  signals  were 
generated  with  a  carrier  frequency  (/c)  of  1000  Hz  and  sampling  frequency  ( fs )  of  7000 

Hz.  The  modulation  bandwidth  (Hz),  the  modulation  period  (ms),  and  the  SNR  are  vari¬ 
able  parameters. 


FMCW 

Modulation 

Bandwidth 

Modulation 

Period 

SNR 

1 

F  1  7  250  20  s 

250  Hz 

20  ms 

Signal  Only 

2 

F  1  7  250  20  0 

250  Hz 

20  ms 

OdB 

3 

F  1  7  250  20  -6 

250  Hz 

20  ms 

-6  dB 

4 

F  1  7  250  30  s 

500  Hz 

30  ms 

Signal  Only 

5 

F  1  7  250  30  0 

500  Hz 

30  ms 

OdB 

6 

F  1  7  250  30  -6 

500  Hz 

30  ms 

-6  dB 

7 

F  1  7  500  20  s 

250  Hz 

20  ms 

Signal  Only 

8 

F  1  7  500  20  0 

250  Hz 

20  ms 

OdB 

9 

F  1  7  500  20  -6 

250  Hz 

20  ms 

-6  dB 

10 

F  1  7  500  30  s 

500  Hz 

30  ms 

Signal  Only 

11 

F  1  7  500  30  0 

500  Hz 

30  ms 

OdB 

12 

F  1  7  500  30  -6 

500  Hz 

30  ms 

-6  dB 

1 

fable  9  FMCW  signals  to  be  processed  by  QMFB  Tree. 
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From  the  list  of  signals  already  processed  by  the  QMFB  (in  Table  9),  only  one  of 
them  is  shown  next.  The  rest  of  the  signals  can  be  seen  in  a  technical  report  that  will  be 
published  soon  [21]. 


a.  F_l_7_500_20_s 


This  FMCW  signal  was  generated  with  the  parameters  described  in  Table 
10.  Therefore,  and  according  with  the  description  of  the  code  given  in  Section  C.l  in 
Chapter  III,  the  signal  is  processed  as  an  input  to  the  QMFB  Tree  to  obtain  the  different 
output  layers.  Thus,  Figure  39  was  created,  which  show  the  output  matrix  from  layer  5 
with  a  frequency  resolution  of  1 12.903  Hz  and  a  time  resolution  of  4.762  ms.  There,  the 
values  of  the  carrier  frequency,  the  modulation  bandwidth,  and  the  modulation  period  are 
extracted,  giving  almost  the  same  values  than  the  original  signal. 


"Time  -  Frequency  Plane  Layer  5' 


1000 

1600 

liCCl 


V  1200 

y  iooo 


300 

600 

400 

ao 

D 


- r 

Cunicr 
F  i  e  qncncjf 


005 


B.inilVd'htlh 


0  1 

Time  fenc} 


0. 15 


Figure  39  Output  matrix  from  layer  5  of  F_l_7_500_20_s  (contourplot). 


Table  10  show  a  summary  of  the  signal  processing  showing  a  comparison 
of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 
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F_l_7_500_20_s 

Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1016.13  kHz 

16.13  Hz  error 

Sampling  Frequency 

7  kHz 

7  kHz 

given 

Modulation  Band¬ 
width 

500  Hz 

564.515  Hz 

64.515  Hz  error 

SNR 

Only  Signal 

- 

Modulation  Period 

20  ms 

23.81  ms 

3.81  ms  error 

Table  10  Signal  processing  summary  for  F_l_7_500_20_s. 


b.  F_l_7_500_20_0 

This  FMCW  signal  was  generated  with  the  parameters  described  in  Table 
1 1 .  Figure  40  was  obtained,  which  provides  a  good  time-frequency  description  of  the 
evaluated  signal  from  the  output  matrix  at  layer  5  when  conducting  the  processing  with 
the  QMFB  tree  with  a  frequency  resolution  of  1 12.903  Hz  and  a  time  resolution  of  4.762 
ms.  There  the  values  of  the  carrier  frequency,  the  modulation  bandwidth,  and  the  modula¬ 
tion  period  are  extracted,  giving  almost  the  same  values  than  the  original  signal. 


Figure  40  Output  matrix  from  layer  5  of F_l_7_500_20_0  (contourplot). 


Table  1 1  shows  a  summary  of  the  signal  processing  demonstrating  a  com¬ 
parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 
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F_l_7_500_20_0 

Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1016.13  kHz 

16.13  Hz  error 

Sampling  Frequency 

7  kHz 

7  kHz 

given 

Modulation  Band¬ 
width 

500  Hz 

677.418  Hz 

177.418  Hz 

error 

SNR 

OdB 

- 

Modulation  Period 

20  ms 

23.81  ms 

3.81  ms  error 

Table  1 1  Signal  processing  summary  for  F_l_7_500_20_0. 


c.  F_l_7_500_20_-6 

This  FMCW  signal  was  generated  with  the  parameters  described  in  Table 
12.  Figure  41  was  obtained,  which  provides  a  good  time-frequency  description  of  the 
evaluated  signal  from  the  output  matrix  at  layer  6  when  conducting  the  processing  with 
the  QMFB  tree  with  a  frequency  resolution  of  55.5556  Hz  and  a  time  resolution  of  9.677 
ms.  There  the  values  of  the  carrier  frequency  and  the  modulation  bandwidth  are  ex¬ 
tracted.  The  extracted  carrier  frequency  matched  the  actual  value;  the  extracted  modula¬ 
tion  bandwidth  was  in  error  by  55.556  Hz.  The  modulation  period  could  not  have  been 
extracted  due  to  the  low  SNR. 


Figure  41  Output  matrix  from  layer  6  of  F_l_7_500_20_-6  (contourplot). 
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Table  12  shows  a  summary  of  the  signal  processing  demonstrating  a  com¬ 
parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


F_l_7_500_20_-6 

Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1  kHz 

Sampling  Frequency 

7  kHz 

7  kHz 

given 

Modulation  Band¬ 
width 

500  Hz 

555.556  Hz 

55.556  Hz 

error 

SNR 

-6  dB 

- 

Modulation  Period 

20  ms 

- 

Table  12  Signal  processing  summary  for  F_l_7_500_20_-6. 


D.  FRANK  CODES 

1.  Brief  Description 

Frank  codes  are  a  family  of  polyphase  codes  that  are  closely  related  to  the  Chirp 
and  Baker  codes  and  have  been  used  successfully  in  LPI  radars  signals.  In  the  case  of  a 
step  approximation  to  linear  frequency  modulation  with  N  frequency  steps,  N  samples  per 
frequency  are  obtained  and  (N)(N)  samples  result.  That  is,  the  Frank  codes  has  a  length 
N2  which  is  also  the  corresponding  pulse  compression  ratio.  In  the  case  of  single  side¬ 
band  detection,  the  result  is  the  Frank  code  [1].  For  example,  if  the  local  oscillator  is  at 
the  start  of  the  sweep  of  the  step  approximation  to  linear  frequency  wavefonn,  the  first  N 
samples  of  the  polyphase  code  are  0  phase.  The  second  N  samples  start  with  0  phase  and 
increase  with  phase  increments  of  (2 n  deg /  N )  from  sample  to  sample,  the  third  group  of 
N  samples  starts  with  0  phase  and  increase  with  {In  deg/  N )  degrees  increments  from 
sample  to  sample  and  so  on.  Figure  42  shows  the  relationship  between  the  index  in  the 
matrix  and  its  phase  shift  for  N2=  16.  If  i  is  the  number  of  the  sample  in  a  given  frequency 
and  j  is  the  number  of  the  frequency,  the  phase  of  the  i-th  sample  of  the  j- th  frequency  is 

$  j  -  ~~ 0  - 1)( 7  - 1)  (3-14) 

N 
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where  i  =  1,2,..., N ,  and  j  =  l,2,...,N  [1],  The  Frank  code  has  the  highest  phase  incre¬ 
ments  from  sample  to  sample  in  the  center  of  the  code.  The  Frank  polyphase  code  can 
also  be  described  by  a  matrix  as  follows 


"0 

0 

0 

0 

0 

1 

2 

...  (N- 1) 

0 

2 

4 

...  2(N  - 1) 

1 

o  • 

5; 

i  ' 

2(iV-l) 

...  (N -l)2  _ 

where  the  numbers  represent  multiplying  coefficients  of  the  basic  phase  angle  2n  I N  [1]. 


Figure  42  Phase  relationship  between  the  index  in  the  matrix  and  its  phase  shift  for 
iV2=T6. 


2.  Processing  Frank  Code  Signals  with  QMFB  Tree. 

The  Frank  code  signals  to  be  worked  with  are  given  in  Table  13.  All  the  signals 
were  generated  with  a  carrier  frequency  (/c)  of  1000  Hz  and  sampling  frequency  (  fs )  of 

7000  Hz.  The  number  of  phases  (N),  the  number  of  cycles  per  phase  (cpp)  and  the  SNR 
are  variable  parameters. 
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FRANK 

Number  of 
code  phases 
N 

Number  of  cycles 
per  phase 
cpp 

SNR 

1 

FR  1  7  4  1  s 

4 

1 

Signal  Only 

2 

FR  1  7  4  1  0 

4 

1 

OdB 

3 

FR  1  7  4  1  -6 

4 

1 

-6  dB 

4 

FR  1  7  4  5  s 

4 

5 

Signal  Only 

5 

FR  1  7  4  5  0 

4 

5 

OdB 

6 

FR  1  7  4  5  -6 

4 

5 

-6  dB 

7 

FR  1  7  8  1  s 

8 

1 

Signal  Only 

8 

FR  1  7  8  1  0 

8 

1 

OdB 

9 

FR  1  7  8  1  -6 

8 

1 

-6  dB 

10 

FR  1  7  8  5  s 

8 

5 

Signal  Only 

11 

FR  1  7  8  5  0 

8 

5 

OdB 

12 

FR  1  7  8  5  -6 

8 

5 

-6  dB 

Table  13  Frank  Code  signals  to  be  processed  by  QMFB  Tree. 


From  the  list  of  signals  already  processed  by  the  QMFB  (in  Table  13),  only  one  of 
them  is  shown  next.  The  rest  of  the  signals  can  be  seen  in  a  technical  report  that  will  be 
published  soon  [21]. 


a.  FR_l_7_4_5_s 


This  Frank  code  signal  was  generated  with  the  parameters  described  in 
Table  14.  The  code  period  of  the  Frank-coded  signal  is 


(number  of  cycles  per  phase)(number  of  phases2) 


(■ cpp){N 2) 

fc 


Carrier  Frequency 


<5M  =  80  ms 
1000 


(3.16) 


The  bandwidth  of  the  signal  depends  on  the  cycles  per  phase  (or  chirp)  as 

BW  =  Carrier  Frequency  =  A  =  =  200  Hz  .  (3.17) 

Cycles  per  phase  cpp  5 

Therefore,  and  according  to  the  description  of  the  code  given  in  Section 
C.l  in  Chapter  III,  the  signal  is  processed  as  an  input  to  the  QMFB  Tree  to  get  the  differ¬ 
ent  output  layers.  Thus,  Figure  43  was  created  which  shows  the  output  matrix  from  layer 
6  with  a  frequency  resolution  of  55.5556  Hz  and  a  time  resolution  of  9.286  ms.  There,  the 
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values  of  the  carrier  frequency,  the  bandwidth,  the  code  period  are  extracted.  The  number 
of  phases  is  extracted  from  Figure  44,  which  shows  a  zoom  on  contourplot  of  layer  2. 


Time  -  Frequency  Plane  LayerS" 


at  Q1  Q3  CU  0-5 
Time  (sbc) 


Figure  43  Output  matrix  from  layer  6  of  FR_l_7_4_5_s  (colorplot). 
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Figure  44  Zoom  on  output  matrix  from  layer  2  of  FR_l_7_4_5_s  (contourplot). 
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Table  14  shows  a  summary  of  the  signal  processing  demonstrating  a  com¬ 
parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


FR  1  7  4  5  s 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1  kHz 

Sampling  Frequency 

7  kHz 

7  kHz 

Given 

Number  of  Phases 

16 

16 

Cycles  per  phase 

5 

- 

SNR 

Only  Signal 

- 

Bandwidth 

200  Hz 

222.222  Hz 

22.222  Hz  error 

Code  Period 

80  ms 

83.574  ms 

3.574  ms  error 

Table  14  Signal  processing  summary  for  FR_l_7_4_5_s. 


b.  FR_1_7_4_5_0 

This  Frank  code  signal  was  generated  with  the  parameters  described  in 
Table  15.  Figure  45  was  obtained,  which  provides  a  good  time-frequency  description  of 
the  evaluated  signal  from  the  output  matrix  at  layer  6  when  conducting  the  processing 
with  the  QMFB  tree  with  a  frequency  resolution  of  55.5556  Hz  and  a  time  resolution  of 
9.286  ms.  There  the  values  of  the  carrier  frequency,  the  bandwidth,  and  the  code  period 
are  extracted,  giving  almost  the  same  values  than  the  original  signal. 


Figure  45  Output  matrix  from  layer  6  of  FR_1_7_4_5_0  (colorplot). 
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Table  15  shows  a  summary  of  the  signal  processing  demonstrating  a  com¬ 
parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


FR  1  7  4  5  0 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1  kHz 

Sampling  Frequency 

7  kHz 

7  kHz 

Given 

Number  of  Phases 

16 

- 

Cycles  per  phase 

5 

- 

SNR 

OdB 

- 

Bandwidth 

200  Hz 

222.222  Hz 

22.222  Hz  error 

Code  Period 

80  ms 

83.574  ms 

3.574  ms  error 

Table  15  Signal  processing  summary  for  FR_1_7_4_5_0. 


c.  FR_l_7_4_5_-6 

This  Frank  code  signal  was  generated  with  the  parameters  described  in 
Table  16.  Figure  46  was  obtained,  which  provides  a  good  time-frequency  description  of 
the  evaluated  signal  from  the  output  matrix  at  layer  6  when  conducting  the  processing 
with  the  QMFB  tree  with  a  frequency  resolution  of  55.5556  Hz  and  a  time  resolution  of 
9.286  ms.  There  the  values  of  the  carrier  frequency,  the  bandwidth,  and  the  code  period 
are  extracted,  giving  almost  the  same  values  than  the  original  signal,  never  the  less,  the 
strong  noise  condition  presented  by  the  signal  does  not  allow  do  it  easily. 
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Figure  46  Output  matrix  from  layer  6  of  FR_l_7_4_5_-6  (colorplot). 

Table  16  shows  a  summary  of  the  signal  processing  demonstrating  a  com¬ 
parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


FR  1  7  4  5  -6 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1  kHz 

Sampling  Frequency 

7  kHz 

7  kHz 

Given 

Number  of  Phases 

16 

- 

Cycles  per  phase 

5 

- 

SNR 

-6dB 

- 

Bandwidth 

200  Hz 

222.222  Hz 

22.222  Hz  error 

Code  Period 

80  ms 

83.574  ms 

3.574  ms  error 

Table  16  Signal  processing  summary  for  FR_l_7_4_5_-6. 


E.  PI  POLYPHASE  CODE 

1.  Brief  Description 

In  case  of  a  double  sideband  detection  (local  oscillator  is  at  band  center)  of  a  step 
approximation  of  a  linear  frequency  modulation,  a  PI  code  results.  The  PI  code  also  con- 
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sists  of  N2  elements  as  Frank  code,  that  way  PI  code  signal  with  N  =  4  produces  a  matrix 
of  16  different  phases,  if  N  =  8  produces  a  matrix  of  64  phases.  If  i  is  the  number  of  sam¬ 
ple  in  a  given  frequency  and  j  is  the  number  of  the  frequency,  the  phase  of  the  i- th  sam¬ 
ple  of  the  j-th  frequency  is 

^y-=^[^-(2;-l)][0-l)iV  +  (i-l)]  (3.18) 

N 

where  i  =  1,2,..., N ,  and  j  =  \,2,...,N  . 


2.  Processing  PI  Code  Signals  with  QMFB  Tree 

The  PI  code  signals  to  be  worked  with  are  given  in  Table  17.  All  the  signals  were 
generated  with  a  carrier  frequency  (/c)of  1000  Hz  and  sampling  frequency  (  fs)  of  7000 

Hz.  The  number  of  phases  (AO,  the  number  of  cycles  per  phase  ( cpp )  and  the  SNR  are 
variable  parameters. 


PI 

Number  of 
code  phases 
N 

Number  of  cycles 
per  phase 
cpp 

SNR 

1 

PI  1  7  4  1  s 

4 

1 

Signal  Only 

2 

PI  1  7  4  1  0 

4 

1 

OdB 

3 

PI  1  7  4  1  -6 

4 

1 

-6  dB 

4 

PI  1  7  4  5  s 

4 

5 

Signal  Only 

5 

PI  1  7  4  5  0 

4 

5 

OdB 

6 

PI  1745-6 

4 

5 

-6  dB 

7 

PI  1  7  8  1  s 

8 

1 

Signal  Only 

8 

PI  1  7  8  1  0 

8 

1 

OdB 

9 

PI  1  7  8  1  -6 

8 

1 

-6  dB 

10 

PI  1  7  8  5  s 

8 

5 

Signal  Only 

11 

PI  1  7  8  5  0 

8 

5 

OdB 

12 

PI  1  7  8  5  -6 

8 

5 

-6  dB 

Table  17  PI  signals  to  be  processed  by  QMFB  Tree. 


From  the  list  of  signals  already  processed  by  the  QMFB  (in  Table  17),  only  one  of 
them  is  shown  next.  The  rest  of  the  signals  can  be  seen  in  a  technical  report  that  will  be 
published  soon  [21]. 
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a.  Pl_l_7_4_5_s 

This  PI  code  signal  was  generated  with  the  parameters  described  in  Table 
18.  The  code  period  of  the  PI  code  signal  is 


(number  of  cycles  per  phase)(number  of  phases2) 


0 cpp)(N 2) 

fe 


Carrier  Frequency 
1000 


(3.19) 


The  bandwidth  of  the  signal  depends  on  the  cycles  per  phase  (or  chirp)  as 


BW  Carrier  Frequency  fe  1000  Hz  =  200  Hz 

Cycles  per  phase  cpp  5 


(3.20) 


Therefore,  and  according  to  the  description  of  the  code  given  in  Section 
C.l  in  Chapter  III,  the  signal  is  processed  as  an  input  to  the  QMFB  tree  to  obtain  the  dif¬ 
ferent  output  layers.  In  this  manner  Figure  47  was  created  which  shows  the  output  matrix 
from  layer  6  with  a  frequency  resolution  of  55.5556  Hz  and  a  time  resolution  of  9.286 
ms.  The  values  of  the  carrier  frequency,  the  bandwidth,  and  the  code  period  are  extracted, 
giving  almost  the  same  values  than  the  original  signal.  Figure  48  shows  a  zoom  of  the 
output  matrix  from  layer  2,  indicating  the  number  of  phases  in  the  signal. 
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Figure  47  Output  matrix  from  layer  7  of Pl_l_7_4_5_s  (colorplot). 
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Figure  48  Zoom  on  output  matrix  from  layer  2  of  Pl_l_7_4_5_s  (contourplot). 

Table  18  shows  a  summary  of  the  signal  processing  demonstrating  a  com¬ 
parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


"rime  -  F-rfu^jnn.;. ^  Plane  Lmfir2‘ 
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PI  1  7  4  5  s 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1  kHz 

Sampling  Frequency 

7  kHz 

7  kHz 

Given 

Number  of  Phases 

16 

16 

Cycles  per  phase 

5 

- 

SNR 

Signal  Only 

- 

Bandwidth 

200  Hz 

222.222  Hz 

22.222  Hz  error 

Code  Period 

80  ms 

83.574  ms 

3.574  ms  error 

Table  18  Signal  processing  summary  for  PI 

7  4  5  s. 

b.  Pl_l_7_4_5_0 


This  PI  code  signal  was  generated  with  the  parameters  described  in  Table 
19.  It  must  be  noted  that  this  signal  has  almost  the  same  parameters  as  the  previous  sig¬ 
nal,  the  only  difference  is  the  SNR  that  now  is  0  dB.  Figure  49  was  obtained,  which  pro¬ 
vides  a  good  time-frequency  description  of  the  evaluated  signal  from  the  output  matrix  at 
layer  6  when  conducting  the  processing  with  the  QMFB  tree  with  a  frequency  resolution 
of  55.5556  Hz  and  a  time  resolution  of  9.286  ms.  There  the  values  of  the  carrier  fre¬ 
quency,  the  bandwidth,  and  the  code  period  were  extracted,  giving  almost  the  same  val¬ 
ues  than  the  original  signal. 


ftid  -  frnfliJMirv  FlariH  Layar  6 


Time  fcee'l 


Figure  49  Output  matrix  from  layer  7  of Pl_l_7_4_5_0  (colorplot). 


68 


Table  19  shows  a  summary  of  the  signal  processing  indicating  a  compari¬ 
son  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


PI  1  7  4  5  0 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1  kHz 

Sampling  Frequency 

7  kHz 

7  kHz 

Given 

Number  of  Phases 

16 

- 

Cycles  per  phase 

5 

- 

SNR 

0  dB 

- 

Bandwidth 

200  Hz 

222.222  Hz 

22.222  Hz  error 

Code  Period 

80  ms 

83.574  ms 

3,574  ms  error 

Table  19  Signal  processing  summary  for  Pl_l_7_4_5_0. 


c.  Pl_l_7_4_5_-6 

This  Pi  code  signal  was  generated  with  the  parameters  described  in  Table 
20.  It  must  be  noted  that  this  signal  has  almost  the  same  parameters  as  the  previous  two 
signals,  the  only  difference  is  the  signal  to  noise  ratio  that  now  is  -6  dB.  Figure  50  was 
obtained,  which  provides  a  good  time-frequency  description  of  the  evaluated  signal  from 
the  output  matrix  at  layer  6  when  conducting  the  processing  with  the  QMFB  tree  with  a 
frequency  resolution  of  55.5556  Hz  and  a  time  resolution  of  9.286  ms.  There  the  values 
of  the  carrier  frequency,  the  bandwidth,  and  the  code  period  are  extracted,  giving  almost 
the  same  values  than  the  original  signal. 
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Figure  50  Output  matrix  from  layer  7  of Pl_l_7_4_5_-6  (colorplot). 

Table  20  shows  a  summary  of  the  signal  processing  demonstrating  a  com¬ 
parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


PI  1  7  4  5  -6 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1  kHz 

Sampling  Frequency 

7  kHz 

7  kHz 

Given 

Number  of  Phases 

16 

- 

Cycles  per  phase 

5 

- 

SNR 

-6dB 

- 

Bandwidth 

200  Hz 

222.222  Hz 

22.222  Hz  error 

Code  Period 

80  ms 

83.574  ms 

3.574  ms  error 

Table  20  Signal  processing  summary  for  Pl_l_7_4_5_-6. 


F.  P2  POLYPHASE  CODE 

1.  Brief  Description 

This  code  is  essentially  derived  in  the  same  way  as  the  PI  code.  The  P2  code  has 
the  same  phase  increments  within  each  group  as  the  PI  code,  except  that  the  starting 
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phase  is  different.  The  P2  code  is  valid  for  N  even,  and  each  group  of  the  code  is  sym¬ 
metric  about  0  phase.  These  phases  can  be  calculated  by 


n 


-[(iV-l)/iV]-[| 


f  n  ^ 
viVy 


(i-1)]  [iV  +  1-2 j\ 


(3.21) 


This  code  has  the  frequency  symmetry  of  the  PI  code.  The  P2  polyphase  code,  as 
well  as  the  PI,  has  more  of  a  symmetrical  spectrum  than  a  Frank-coded  signal  due  to  its 
symmetry  in  the  carrier. 


2.  Processing  P2  Code  Signals  with  QMFB  Tree 

The  P2  code  signals  to  be  worked  with  are  given  in  Table  21.  All  the  signals  were 
generated  with  a  carrier  frequency  (/c)  of  1000  Hz  and  sampling  frequency  (  fs )  of  7000 

Hz.  The  number  of  phases  (AO,  the  number  of  cycles  per  phase  ( cpp ),  and  the  SNR  are 
variable  parameters. 


P2 

Number  of 
code  phases 
N 

Number  of  cycles 
per  phase 
cpp 

SNR 

1 

P2  1  7  4  1  s 

4 

1 

Signal  Only 

2 

P2  1  7  4  1  0 

4 

1 

OdB 

3 

P2  1  7  4  1  -6 

4 

1 

-6  dB 

4 

P2  1  7  4  5  s 

4 

5 

Signal  Only 

5 

P2  1  7  4  5  0 

4 

5 

OdB 

6 

P2  1  7  4  5  -6 

4 

5 

-6  dB 

7 

P2  1  7  8  1  s 

8 

1 

Signal  Only 

8 

P2  1  7  8  1  0 

8 

1 

OdB 

9 

P2  1  7  8  1  -6 

8 

1 

-6  dB 

10 

P2  1  7  8  5  s 

8 

5 

Signal  Only 

11 

P2  1  7  8  5  0 

8 

5 

OdB 

12 

P2  1  7  8  5  -6 

8 

5 

-6  dB 

Table  21  P2  signals  to  be  processed  by  QMFB  Tree. 


From  the  list  of  signals  already  processed  by  the  QMFB  (in  Table  21),  only  one  of 
them  is  shown  next.  The  rest  of  the  signals  can  be  seen  in  a  technical  report  that  will  be 
published  soon  [21], 
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a.  P2_l_7_4_5_s 

This  P2  code  signal  was  generated  with  the  parameters  described  in  Table 
22.  The  code  period  of  the  PI  code  signal  is 


(number  of  cycles  per  phase)(number  of  phases2) 
Carrier  Frequency 

(cppkn  yiw  =  80m 


(3.22) 


The  bandwidth  of  the  signal  depends  on  the  cycles  per  phase  (or  chirp)  as 


BW  Carrier  Frequency  fe  1000  Hz  =  200  Hz 

Cycles  per  phase  cpp  5 


(3.23) 


Now,  and  according  with  the  description  of  the  code  given  in  Section  C.l 
in  Chapter  III,  the  signal  is  processed  as  an  input  to  the  QMFB  Tree  to  get  the  different 
output  layers.  In  this  manner  Figure  51  was  created  which  show  the  output  matrix  from 
layer  7  with  a  frequency  resolution  of  27.777  Hz  and  a  time  resolution  of  18.870  ms.  The 
values  of  the  carrier  frequency,  the  bandwidth,  and  the  code  period  are  extracted,  giving 
almost  the  same  values  than  the  original  signal.  Figure  52  shows  the  output  matrix  layer  2 
in  a  contourplot  indicating  the  number  of  phases  in  the  signal. 
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Figure  52  Zoom  of  output  matrix  from  layer  2  of P2_l_7_4_5_s  (contourplot). 


Table  22  shows  a  summary  of  the  signal  processing  indicating  a  compari¬ 
son  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 
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P2  1  7  4  5  s 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

999.972  Hz 

0.028  Hz  error 

Sampling  Frequency 

7  kHz 

7  kHz 

Given 

Number  of  Phases 

16 

16 

Cycles  per  phase 

5 

- 

SNR 

Only  Signal 

- 

Bandwidth 

200  Hz 

194.439  Hz 

5.561  Hz  error 

Code  Period 

80  ms 

75.48  ms 

4.52  ms  error 

Table  22  Signal  processing  summary  for  P2_l_7_4_5_s. 


b.  P2_l_7_4_5_0 

This  P2  code  signal  was  generated  with  the  parameters  described  in  Table 
23.  It  must  be  noted  that  this  signal  has  almost  the  same  parameters  as  the  previous  sig¬ 
nal,  the  only  difference  is  the  SNR  that  now  is  0  dB.  Figure  53  was  obtained,  which  pro¬ 
vides  a  good  time-frequency  description  of  the  evaluated  signal  from  the  output  matrix  at 
layer  7  when  conducting  the  processing  with  the  QMFB  tree  with  a  frequency  resolution 
of  27.777  Hz  and  a  time  resolution  of  18.870  ms.  There  the  values  of  the  carrier  fre¬ 
quency,  the  bandwidth,  and  the  code  period  are  extracted,  giving  almost  the  same  values 
than  the  original  signal. 


Figure  53  Output  matrix  from  layer  7  of P2_l_7_4_5_0  (colorplot). 
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Table  23  shows  a  summary  of  the  signal  processing  indicating  a  compari¬ 
son  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


P2  1  7  4  5  0 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

999.972  Hz 

0.028  Hz  error 

Sampling  Frequency 

7  kHz 

7  kHz 

Given 

Number  of  Phases 

16 

- 

Cycles  per  phase 

5 

- 

SNR 

0  dB 

- 

Bandwidth 

200  Hz 

194.439  Hz 

5.561  Hz  error 

Code  Period 

80  ms 

75.48  ms 

4.52  ms  error 

Table  23  Signal  processing  summary  for  P2_l_7_4_5_0. 


c.  P2_l_7_4_5_-6 

This  P2  code  signal  was  generated  with  the  parameters  described  in  Table 
24.  It  must  be  noted  that  this  signal  has  almost  the  same  parameters  as  the  previous  two 
signals,  the  only  difference  is  the  signal  to  noise  ratio  that  now  is  -6  dB.  Figure  54  was 
obtained,  which  provides  a  good  time-frequency  description  of  the  evaluated  signal  from 
the  output  matrix  at  layer  7  when  conducting  the  processing  with  the  QMFB  tree  with  a 
frequency  resolution  of  27.777  Hz  and  a  time  resolution  of  18.870  ms.  There  the  values 
of  the  carrier  frequency,  the  bandwidth,  and  the  code  period  are  extracted,  giving  almost 
the  same  values  than  the  original  signal. 
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Figure  54  Output  matrix  from  layer  7  of P2_l_7_4_5_-6  (colorplot). 


Table  24  shows  a  summary  of  the  signal  processing  demonstrating  a  com¬ 
parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


P2  1  7  4  5  -6 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

0.999  kHz 

0.028  Hz  error 

Sampling  Frequency 

7  kHz 

7  kHz 

Given 

Number  of  Phases 

16 

- 

Cycles  per  phase 

5 

- 

SNR 

-6dB 

- 

Bandwidth 

200  Hz 

166.662  Hz 

33.338  Hz  error 

Code  Period 

80  ms 

75.48  ms 

4.52  ms  error 

Table  24  Signal  processing  summary  for  P2_l_7_4_5_-6. 


G.  P3  POLYPHASE  CODE 

1.  Brief  Description 

This  code  is  derived  by  converting  a  linear-frequency  modulation  waveform  to 
baseband  using  a  local  oscillator  on  one  end  of  the  frequency  sweep  and  sampling  the  In- 
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phase  I  and  Quadrature  Q  video  at  the  Nyquist  rate.  If  it  is  assumed  that  the  waveform 
has  a  pulse  length  T  in  frequency  f  =  fo+kt ,  where  k  is  a  constant,  the  bandwidth  of  the 
signal  will  be  approximately  B=kT.. 

The  bandwidth  will  support  a  compressed  pulse  length  of  about  tc=\/  B  and  the 
waveform  will  provide  a  pulse  compression  ratio  of  pc  —  T  /  tc  =  BT  . 

Assuming  that  the  first  sample  of  I  and  Q  is  taken  at  the  leading  edge  of  the  wave¬ 
form,  the  phases  of  successive  samples  taking  tc  apart  are: 

o-i  )tc 

ti=2  tt  J  [(f0+kt)-f0]dt  =  7Tk(i-\)2t2c  (3.24) 
0 


where  i=  1,2,  . . .,  N.  Substituting  k=B/T  and  tc=  \/B,  the  equation  can  be  written  as 

,  _  7i (i  - 1)2  _  7i (i  - 1)2 

®  ~  BT  N 


(3.25) 


2.  Processing  P3  Code  Signals  with  QMFB  Tree 

The  P3  code  signals  to  be  worked  with  are  given  in  Table  25.  All  the  signals  were 
generated  with  a  carrier  frequency  (/c)  of  1000  Hz  and  sampling  frequency  (  fs )  of  7000 

Hz.  The  number  of  phases  (AO,  the  number  of  cycles  per  phase  ( cpp )  and  the  SNR  are 
variable  parameters. 
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P3 

Number  of 
code  phases 
N2 

Number  of  cycles 
per  phase 
cpp 

SNR 

1 

P3  1  7  16  1  s 

16 

1 

Signal  Only 

2 

P3  1  7  16  1  0 

16 

1 

OdB 

3 

P3  1  7  16  1  —6 

16 

1 

-6  dB 

4 

P3  1  7  16  5  s 

16 

5 

Signal  Only 

5 

P3  1  7  16  5  0 

16 

5 

OdB 

6 

P3  1  7  16  5  -6 

16 

5 

-6  dB 

7 

P3  1  7  64  1  s 

64 

1 

Signal  Only 

8 

P3  1  7  64  1  0 

64 

1 

OdB 

9 

P3  1  7  64  1  -6 

64 

1 

-6  dB 

10 

P3  1  7  64  5  s 

64 

5 

Signal  Only 

11 

P3  1  7  64  5  0 

64 

5 

OdB 

12 

P3  1  7  64  5  -6 

64 

5 

-6  dB 

Table  25  P3  signals  to  be  processed  by  QMFB  Tree. 


From  the  list  of  signals  already  processed  by  the  QMFB  (in  Table  25),  only  one  of 
them  is  shown  next.  The  rest  of  the  signals  can  be  seen  in  a  technical  report  that  will  be 
published  soon  [21]. 


a.  P3_l_7_l  6_5_s 

This  P3  code  signal  was  generated  with  the  parameters  described  in  Table 
26.  The  code  period  of  the  P3  code  signal  is 


T  _  (number  of  cycles  per  phase)(number  of  phases2) 

3  c  | 

Carrier  Frequency 
_  (cpp)(N2)  _  (5X16)  _  8Q  ms 
fc  1000 


(3.26) 


The  bandwidth  of  the  signal  depends  on  the  cycles  per  phase  (or  chirp)  as 

BW  =  Carrier  Frequency  =  =  1000  Hz  =  ^  Rz  (3  27) 

Cycles  per  phase  cpp  5 

Now,  and  according  with  the  description  of  the  code  given  in  Section  C.l 
in  Chapter  III,  the  signal  is  processed  as  an  input  to  the  QMFB  tree  to  get  the  different 
output  layers.  In  this  manner,  Figure  55  was  created  which  shows  the  output  matrix  from 
layer  6  with  a  frequency  resolution  of  55.5556  Hz  and  a  time  resolution  of  9.286  ms.  The 
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values  of  the  carrier  frequency,  the  bandwidth,  and  the  code  period  are  extracted,  giving 
almost  the  same  values  than  the  original  signal.  Figure  56  shows  a  zoom  of  the  output 
matrix  layer  2  indicating  the  number  of  phases  in  the  signal. 
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Figure  55  Output  matrix  from  layer  7  of  P3_l_7_16_5_s  (colorplot). 
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Figure  56  Zoom  of  output  matrix  layer  2  of  P3_l_7_16_5_s  (contourp 
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ot). 


Table  26  shows  a  summary  of  the  signal  processing  demonstrating  a  com¬ 
parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


P3  1  7  16  5  s 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1  kHz 

Sampling  Frequency 

7  kHz 

7  kHz 

Given 

Number  of  Phases 

16 

16 

Cycles  per  phase 

5 

- 

SNR 

Only  Signal 

- 

Bandwidth 

200  Hz 

222.222  Hz 

22.222  Hz  error 

Code  Period 

80  ms 

83.574  ms 

3.574  ms  error 

Table  26  Signal  processing  summary  for  P3_l_7_16_5_s. 


b.  P3_l_7_16_5_0 

This  P3  code  signal  was  generated  with  the  parameters  described  in  Table 
27.  It  must  be  noted  that  this  signal  has  almost  the  same  parameters  than  the  previous 
signal,  the  only  difference  is  the  signal  to  noise  ratio  that  now  is  0  dB.  Figure  57  was  ob¬ 
tained,  which  provides  a  good  time-frequency  description  of  the  evaluated  signal  from 
the  output  matrix  at  layer  6  when  conducting  the  processing  with  the  QMFB  tree  with  a 
frequency  resolution  of  55.5556  Hz  and  a  time  resolution  of  9.286  ms.  The  values  of  the 
carrier  frequency,  the  bandwidth,  and  the  code  period  are  extracted,  giving  almost  the 
same  values  than  the  original  signal. 
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Figure  57  Output  matrix  from  layer  7  of  P3_l_7_16_5_0  (colorplot). 


Table  27  shows  a  summary  of  the  signal  processing  demonstrating  a  com¬ 
parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


P3  1  7  16  5  0 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1  kHz 

Sampling  Frequency 

7  kHz 

7  kHz 

Given 

Number  of  Phases 

16 

- 

Cycles  per  phase 

5 

- 

SNR 

0  dB 

- 

Bandwidth 

200  Hz 

222.222  Hz 

22.222  Hz  error 

Code  Period 

80  ms 

92.86  ms 

12.86  ms  error 

Table  27  Signal  processing  summary  for  P3_l_7_16_5_0. 


c.  P3_l_7_l  6_5_-6 

This  P3  code  signal  was  generated  with  the  parameters  described  in  Table 
28.  It  must  be  noted  that  this  signal  has  almost  the  same  parameters  as  the  previous  two 
signals,  the  only  difference  is  the  signal  to  noise  ratio  that  now  is  -6  dB.  Figure  58  was 
obtained,  which  provides  a  good  time-frequency  description  of  the  evaluated  signal  from 
the  output  matrix  at  layer  6  when  conducting  the  processing  with  the  QMFB  tree  with  a 
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frequency  resolution  of  55.5556  Hz  and  a  time  resolution  of  9.286  ms.  There  the  values 
of  the  carrier  frequency,  the  bandwidth,  and  the  code  period  are  extracted,  giving  almost 
the  same  values  than  the  original  signal. 


Figure  58  Output  matrix  from  layer  7  of P3_l_7_ 165-6  (colorplot). 


Table  28  shows  a  summary  of  the  signal  processing  demonstrating  a  com¬ 
parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


P3  1  7  16  5  -6 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1  kHz 

Sampling  Frequency 

7  kHz 

7  kHz 

Given 

Number  of  Phases 

16 

- 

Cycles  per  phase 

5 

- 

SNR 

-6  dB 

- 

Bandwidth 

200  Hz 

222.222  Hz 

22.222  Hz  error 

Code  Period 

80  ms 

92.86  ms 

12.86  ms  error 

Table  28  Signal  processing  summary  for  P3_l_7_16_5_— 6. 
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H.  P4  POLYPHASE  CODE 


1.  Brief  Description 

Conceptual  coherent  double  sideband  detection  of  the  linear  frequency  modula¬ 
tion  waveform  and  sampling  at  the  Nyquist  rate  yields  a  polyphase  code  named  the  P4. 
The  P4  code  consists  of  discrete  phases  of  the  linear  chirp  wavefonn  taken  at  specific 
time  intervals  and  exhibits  the  same  range  Doppler  coupling  associate  with  the  chirp 
waveform.  However,  the  peak  sidelobe  levels  are  lower  than  those  of  the  unweighted 
chirp  wavefonn.  Various  weighting  techniques  can  be  applied  to  reduce  the  sidelobe  lev¬ 
els  further.  Phase  code  elements  of  the  P4  code  are  given  by  [1] 

*=£^-,(,■-1)  (3.28) 

for  i  =  1  to  N. 


2.  Processing  P4  Code  Signals  with  QMFB  Tree 

The  P4  code  signals  to  be  worked  with  are  given  in  Table  29.  All  the  signals  were 
generated  with  a  carrier  frequency  (/c)  of  1000  Hz  and  sampling  frequency  (  fs )  of  7000 

Hz.  The  number  of  phases  (AO,  the  number  of  cycles  per  phase  ( cpp )  and  the  SNR  are 
variable  parameters. 


83 


P4 

Number  of 
code  phases 
N2 

Number  of  cycles 
per  phase 
cpp 

SNR 

1 

P4  1  7  16  1  s 

16 

1 

Signal  Only 

2 

P4  1  7  16  1  0 

16 

1 

OdB 

3 

P4  1  7  16  1  -6 

16 

1 

-6  dB 

4 

P4  1  7  16  5  s 

16 

5 

Signal  Only 

5 

P4  1  7  16  5  0 

16 

5 

OdB 

6 

P4  1  7  16  5  -6 

16 

5 

-6  dB 

7 

P4  1  7  64  1  s 

64 

1 

Signal  Only 

8 

P4  1  7  64  1  0 

64 

1 

OdB 

9 

P4  1  7  64  1  -6 

64 

1 

-6  dB 

10 

P4  1  7  64  5  s 

64 

5 

Signal  Only 

11 

P4  1  7  64  5  0 

64 

5 

OdB 

12 

P4  1  7  64  5  -6 

64 

5 

-6  dB 

Table  29  P4  signals  to  be  processed  by  QMFB  Tree. 


From  the  list  of  signals  already  processed  by  the  QMFB  (in  Table  29),  only  one  of 
them  is  shown  next.  The  rest  of  the  signals  can  be  seen  in  a  technical  report  that  will  be 
published  soon  [21]. 


a.  P4_l_7_l  6_5_s 

This  P4  code  signal  was  generated  with  the  parameters  described  in  Table 
30.  The  code  period  of  the  P4  code  signal  is 


T  _  (number  of  cycles  per  phase)(number  of  phases2) 

4  c  | 

Carrier  Frequency 
_  (cpp)(N2)  _  (5X16)  _  8Q  ms 
fc  1000 


(3.29) 


The  bandwidth  of  the  signal  depends  on  the  cycles  per  phase  (or  chirp)  as 

BW  =  Carrier  Frequency  =  =  1000  Hz  =  ^  Rz  (J  JQ) 

Cycles  per  phase  cpp  5 

Now,  and  according  with  the  description  of  the  code  given  in  Section  C.l 
in  Chapter  III,  the  signal  is  processed  as  an  input  to  the  QMFB  Tree  to  obtain  the  differ¬ 
ent  output  layers.  In  this  manner,  Figure  59  was  created  which  show  the  output  matrix 
from  layer  6  with  a  frequency  resolution  of  55.5556  Hz  and  a  time  resolution  of  9.286 
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ms.  The  values  of  the  carrier  frequency,  the  bandwidth,  and  the  code  period  are  extracted, 
giving  almost  the  same  values  than  the  original  signal.  Figure  60  shows  a  zoom  of  output 
matrix  layer  2  indicating  the  number  of  phases  in  the  signal. 


"Trim  -  fraquenn}'  Fljnn  Laynr6‘ 
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Figure  59  Output  matrix  from  layer  7  of P4_l_7_16_5_s  (colorplot). 
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Figure  60  Zoom  of  output  matrix  layer  2  of P4_l_7_16_5_s  (colorplot). 
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Table  30  shows  a  summary  of  the  signal  processing  demonstrating  a  com¬ 
parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


P4  1  7  16  5  s 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1  kHz 

Sampling  Frequency 

7  kHz 

7  kHz 

Given 

Number  of  Phases 

16 

16 

Cycles  per  phase 

5 

- 

SNR 

Only  Signal 

- 

Bandwidth 

200  Hz 

222.222  Hz 

22.222  Hz  error 

Code  Period 

80  ms 

83.574  ms 

3.574  ms  error 

Table  30  Signal  processing  summary  for  P4_l_7_16_5_s. 


b.  P4_l_7_l  6_5_0 

This  P4  code  signal  was  generated  with  the  parameters  described  in  Table 
3 1 .  It  must  be  noted  that  this  signal  has  almost  the  same  parameters  as  the  previous  sig¬ 
nal,  the  only  difference  is  the  SNR  that  now  is  0  dB.  Figure  61  was  obtained,  which  pro¬ 
vides  a  good  time-frequency  description  of  the  evaluated  signal  from  the  output  matrix  at 
layer  6  when  conducting  the  processing  with  the  QMFB  tree  with  a  frequency  resolution 
of  55.5556  Hz  and  a  time  resolution  of  9.286  ms.  There  the  values  of  the  carrier  fre¬ 
quency,  the  bandwidth,  and  the  code  period  are  extracted,  giving  almost  the  same  values 
than  the  original  signal. 
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Figure  61  Output  matrix  from  layer  7  of P4_l_7_16_5_0  (colorplot). 


Table  31  shows  a  summary  of  the  signal  processing  demonstrating  a  com¬ 
parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


P4  1  7  16  5  0 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1  kHz 

Sampling  Frequency 

7  kHz 

7  kHz 

Given 

Number  of  Phases 

16 

16 

Cycles  per  phase 

5 

- 

SNR 

0  dB 

- 

Bandwidth 

200  Hz 

222.222  Hz 

22.222  Hz  error 

Code  Period 

80  ms 

83.574  ms 

3.574  ms  error 

Table  3 1  Signal  processing  summary  for  P4_l_7_16_5_0. 


c.  P4_l_7_l  6_5_-6 

This  P4  code  signal  was  generated  with  the  parameters  described  in  Table 
32.  It  must  be  noted  that  this  signal  has  almost  the  same  parameters  as  the  previous  two 
signals,  the  only  difference  is  the  signal  to  noise  ratio  that  now  is  -6  dB.  Figure  62  was 
obtained,  which  provides  a  good  time-frequency  description  of  the  evaluated  signal  from 
the  output  matrix  at  layer  6  when  conducting  the  processing  with  the  QMFB  tree  with  a 
frequency  resolution  of  55.5556  Hz  and  a  time  resolution  of  9.286  ms.  The  values  of  the 
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carrier  frequency,  the  bandwidth,  and  the  code  period  are  extracted,  giving  almost  the 
same  values  than  the  original  signal. 


Figure  62  Output  matrix  from  layer  7  of P3_l_7_ 165-6  (colorplot). 


Table  32  shows  a  summary  of  the  signal  processing  demonstrating  a  com¬ 
parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


P4  1  7  16  5  -6 
Parameters 

Generation 

Detection 

Comment 

Carrier  Frequency 

1  kHz 

1  kHz 

Sampling  Frequency 

7  kHz 

7  kHz 

Given 

Number  of  Phases 

16 

- 

Cycles  per  phase 

5 

- 

SNR 

-6dB 

- 

Bandwidth 

200  Hz 

222.222  Hz 

22.222  Hz  error 

Code  Period 

80  ms 

83.574  ms 

3.574  ms  error 

Table  32  Signal  processing  summary  for  P3_l_7_16_5_— 6. 
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I. 


COSTAS  CODE 


1.  Brief  Description 

In  a  frequency  hopping  system,  the  signal  consists  of  one  or  more  frequencies  be¬ 
ing  chosen  from  a  set  j/1,/2,.. .,//«}  of  available  frequencies,  for  transmission  at  each  of 

a  set  |Vl,r2,...,m}of  consecutives  time  intervals.  For  modeling  purposes,  it  is  reasonable 
to  consider  the  situation  in  which  m  =  n ,  and  a  different  one  of  n  equally  spaced  fre¬ 
quencies  j/1,/2, ...,/n}  is  transmitted  during  each  of  the  n  equal  duration  time  intervals 

j?l,f2,...,tn}  .  Such  a  signal  is  represented  by  a  nxn  permutation  matrix  A,  where  the  n 

rows  correspond  to  the  n  frequencies,  the  n  columns  correspond  to  the  n  intervals,  and 
the  entry  atj  equals  1  means  transmission  and  0  otherwise  [7]. 


Figure  63  Binary  matrix  representation  of  (a)  quantized  linear  FM  and  (b)  Costas  Signal. 


This  signifies  that,  at  any  given  time,  a  slotone  frequency  is  transmitted,  and  each 
frequency  is  transmitted  only  one  (Figure  63(a)).  Other  possible  frequency-hopping 
sequences  that  belong  to  this  family.  This  hopping  order  stongly  effects  the  ambiguity 
function  of  these  signals.  Frequency-hopping  signals  allow  a  simple  procedure  that 
results  in  a  rough  approximation  of  their  ambiguity  function.  This  is  possible  because  the 
cross  correlation  signals  at  different  frequencies  approaches  zero  when  the  frequence 
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difference  is  large  relative  to  the  inverse  of  the  signal  duration.  The  ambiguity  function, 
at  any  given  coordinates,  is  an  integral  of  the  product  between  the  original  signal  and  a 
replica  of  it,  which  is  shifted  in  time  and  frequency  according  to  the  delay  and  the 
Doppler  coordinates  of  the  function. 

Performing  an  exercise  on  the  matrix  in  Figure  63(b),  results  show  that  except  for 
the  zero-shift  cases,  when  the  number  of  coincidences  is  N,  finding  a  combination  of 
shifts  yielding  more  than  one  coincidence  is  not  possible.  This  is  actually  the  criteria  of 
Costas  sequences,  where  sequences  of  frequency  hopping  yield  no  more  than  one 

coincidence.  For  example:  if  jay.}=  4,7, 1,6, 5, 2, 3  is  a  Costas  sequence,  then  its  coding 

matrix  and  difference  matrix  are  shown  in  Figure  64. 
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Figure  64  The  coding  matrix,  different  matrix  and  ambiguity  sidelobes  matrix  of  a 
Costas  signal. 
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2.  Processing  Costas  Code  Signals  with  QMFB  Tree 

The  Costas  code  signals  to  be  worked  with  are  given  in  Table  33.  All  the  signals 
were  generated  with  a  carrier  frequency  (/c)  of  1000  or  3000  Hz,  a  sampling  frequency 

(fs)  of  15000  Hz,  and  the  frequency  sequence  4-7-1-6-5-2-3  (each  of  this  numbers  repre¬ 
sents  frequencies  in  kHz),  with  transmission  time  in  each  frequency  equal  to  10  or  20  ms. 
The  number  of  cycles  per  phase  ( cpp )  and  the  SNR  are  variable  parameters. 


COSTAS 

Sequence 

Number  of  cycles 
per  phase 
cpp 

SNR 

1 

C  1  15  10  s 

4716523 

10 

Signal  Only 

2 

C  1  15  10  0 

4716523 

10 

OdB 

3 

C  1  15  10  -6 

4716523 

10 

-6  dB 

4 

C  1  15  20  s 

4716523 

20 

Signal  Only 

5 

C  1  15  20  0 

4716523 

20 

OdB 

6 

C  1  15  20  -6 

4716523 

20 

-6  dB 

7 

C  2  17  10  s 

2638751 

10 

Signal  Only 

8 

C  2  17  10  0 

2638751 

10 

OdB 

9 

C  2  17  10  -6 

2638751 

10 

-6  dB 

10 

C  2  17  20  s 

2638751 

20 

Signal  Only 

11 

C  2  17  20  0 

2638751 

20 

OdB 

12 

C  2  17  20  -6 

2638751 

20 

-6  dB 

Table  33  Costas  signals  to  be  processed  by  QMFB  Tree. 


From  the  list  of  signals  already  processed  by  the  QMFB  (in  Table  33),  only  one  of 
them  is  shown  next.  The  rest  of  the  signals  can  be  seen  in  a  technical  report  that  will  be 
published  soon  [21]. 


a.  C_l_15_10_s 

This  Costas  code  signal  was  generated  with  the  parameters  described  in 
Table  34.  The  frequency  sequence  is  the  sequence  number  1  (4-7-1-6-5-2-3  kHz)  with  10 
cycles  per  phase.  The  transmission  time  in  each  frequency  of  the  sequence  is  equal  to  10 
ms.  Now,  and  according  with  the  description  of  the  code  given  in  Section  C.l  in  Chapter 
III,  the  signal  is  processed  as  an  input  to  the  QMFB  Tree  to  obtain  the  different  output 
layers.  In  this  manner  Figure  65  was  created  which  show  the  output  matrix  from  layer  6 
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with  a  frequency  resolution  of  124.754  Hz.  A  zoom  of  a  previous  layer  output  matrix  is 
shown  in  Figure  66;  the  values  of  the  Costas  frequency  sequence  are  identified,  as  well  as 
the  code  period  and  bandwidth,  giving  the  same  values  than  the  original  signal. 
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Figure  65  Output  matrix  from  layer  6  of  C_  11 5_  1  O  s  (colorplot). 
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Figure  66  Zoom  of  output  matrix  from  layer  5  of  C_  11 5_  1  O  s  (colorplot). 
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Table  34  presents  a  summary  of  the  signal  processing  indicating  a  com¬ 


parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


C_l_15_10_s 

Parameters 

Generation 

Detection 

Comment 

Costas  Sequence 

4-7-1-6-5-2-3  kHz 

4-7-1-6-5-2-3  kHz 

Sampling  Frequency 

15  kHz 

15  kHz 

Given 

SNR 

Only  Signal 

- 

Transmission  time  per 
Frequency 

10  ms 

10  ms 

Code  Period 

70  ms 

70  ms 

Bandwidth 

6000  Hz 

6112.978  Hz 

1 12.978  Hz  error 

Table  34  Signal  processing  summary  for  C_l_l 51  O  s. 


b.  C_l_15_10_0 

This  Costas  code  signal  was  generated  with  the  parameters  described  in 
Table  35.  It  must  be  noted  that  this  signal  has  almost  the  same  parameters  as  the  previous, 
the  only  difference  is  the  SNR  that  now  is  0  dB.  Figure  67  was  obtained,  which  provides 
a  good  time-frequency  description  of  the  evaluated  signal  from  the  output  matrix  at  layer 
6  when  conducting  the  processing  with  the  QMFB  tree  with  a  frequency  resolution  of 
124.754  Hz.  A  zoom  of  a  previous  layer  output  matrix  is  shown  in  Figure  68;  the  values 
of  the  Costas  frequency  sequence  are  identified,  as  well  as  the  code  period  and  band¬ 
width,  giving  almost  the  same  values  than  the  original  signal. 
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Figure  67  Output  matrix  from  layer  7  of  C_  11 5_  1 0  0  (colorplot). 
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Figure  68  Zoom  of  output  matrix  from  layer  5  of  C_  11 5_  1 0  0  (colorplot). 
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Table  35  presents  a  summary  of  the  signal  processing  demonstrating  a 
comparison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB 


tree. 


C_l_15_10_0 

Parameters 

Generation 

Detection 

Comment 

Costas  Sequence 

4-7-1-6-5-2-3  kHz 

4-7-1-6-5-2-3  kHz 

Sampling  Frequency 

15  kHz 

15  kHz 

Given 

SNR 

0  dB 

- 

Transmission  time  per 
Frequency 

10  ms 

10  ms 

Code  Period 

70  ms 

70  ms 

Bandwidth 

6000  Hz 

6112.978  Hz 

1 12.978  Hz  error 

Table  35  Signal  processing  summary  for  C_  11 5_  1 0  0 . 


c.  C_l_15_10_-6 

This  Costas  code  signal  was  generated  with  the  parameters  described  in 
Table  36.  It  must  be  noted  that  this  signal  has  almost  the  same  parameters  as  the  previous 
two  signals,  the  only  difference  is  the  signal  to  noise  ratio  that  now  is  -6  dB.  Figure  69 
was  obtained,  which  provides  a  good  time-frequency  description  of  the  evaluated  signal 
from  the  output  matrix  at  layer  6  when  conducting  the  processing  with  the  QMFB  tree 
with  a  frequency  resolution  of  124.754  Hz.  A  zoom  of  a  previous  layer  output  matrix  is 
shown  in  Figure  70;  the  values  of  the  Costas  frequency  sequence  are  hardly  identified,  as 
well  as  the  code  period  and  bandwidth,  giving  roughly  the  same  values  than  the  original 
signal. 
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Figure  69  Output  matrix  from  layer  7  of  C_  11 5_  1 0  —6  (colorplot). 
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Figure  70  Zoom  of  output  matrix  from  layer  5  of  C_l_l  51 0  -6  (colorplot). 


Table  36  shows  a  summary  of  the  signal  processing  demonstrating  a  com¬ 
parison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 
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C_l_15_10_-6 

Parameters 

Generation 

Detection 

Comment 

Costas  Sequence 

4-7-1-6-5-2-3  kHz 

4-7-1-6-5-2-3  kHz 

Sampling  Frequency 

15  kHz 

15  kHz 

Given 

SNR 

-6  dB 

- 

Transmission  time  per 
Frequency 

10  ms 

10  ms 

Code  Period 

70  ms 

70  ms 

Bandwidth 

6000  Hz 

6112.978  Hz 

1 12.978  Hz  error 

Table  36  Signal  processing  summary  for  C_  1  _  1 5_  1 0  —6 . 


J.  FREQUENCY  SHIFT  KEYING/PHASE  SHIFT  KEYING  COMBINED 

WITH  COSTAS  CODE  (FSK/PSK  COSTAS) 

1.  Brief  Description 

This  modulation  technique  is  the  result  of  a  combination  of  frequency  shift  keying 
based  on  a  Costas  frequency  hopping  matrix  and  phase  shift  keying  using  Barker  se¬ 
quences  of  different  lengths.  In  a  Costas  frequency  hopped  signal,  the  firing  order  defines 
what  frequencies  will  appear  and  with  what  duration.  During  each  sub-period,  as  the  sig¬ 
nal  stays  in  one  of  the  frequencies,  a  binary  phase  modulation  occurs  according  to  a 
Barker  sequence  of  length  five,  seven,  eleven  or  thirteen  bits. 

The  final  waveform  may  be  seen  as  a  binary  phase  shifting  modulation  within 
each  frequency  hop,  resulting  in  25  phase  slots  equally  distributed  in  each  frequency.  As 
illustrated  in  Figure  71,  if  we  consider  NF  as  the  number  of  frequency  hops  and  Np  as 

the  number  of  phase  slots  of  duration  Tp  in  each  sub-period  TF ,  the  total  number  of  phase 
slots  in  the  FSK/PSK  waveform  is  given  by: 

N  =  NfNp.  (3.31) 
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Figure  7 1  General  FSK/PSK  Signal  containing  Nf  frequency  hops  with  Np  phase  slots 
per  frequency. 


2.  Processing  FSK/PSK  Costas  Code  Signals  with  QMFB  Tree 

The  FSK/PSK  Costas  code  signals  to  be  worked  with  are  given  in  Table  33.  All 
the  signals  were  generated  with  a  Costas  sequence  of  7  frequencies  (4-7-1-6-5-2-3  kHz), 
sampling  frequency  equal  to  15000  Hz,  the  number  of  bits  Barker  code,  the  cycles  per  bit, 
and  the  SNR  are  variables  parameters. 


FSK/PSK  Costas 

Sequence 

Barker 

Bits 

Number  of 
cycles  per  bit 
cpp 

SNR 

1 

FSK  PSK  C  1  15  5  1  s 

4716523 

5 

1 

Signal  Only 

2 

FSK  PSK  C  1  15  5  1  0 

4716523 

5 

1 

OdB 

3 

FSK  PSK  C  1  15  5  5  s 

4716523 

5 

5 

Signal  Only 

4 

FSK  PSK  C  1  15  5  5  0 

4716523 

5 

5 

OdB 

5 

FSK  PSK  C  1  15  11  1  s 

4716523 

11 

1 

Signal  Only 

6 

FSK  PSK  C  1  15  11  1  0 

4716523 

11 

1 

OdB 

7 

FSK  PSK  C  1  15  11  5  s 

4716523 

11 

5 

Signal  Only 

8 

FSK  PSK  C  1  15  11  5  0 

4716523 

11 

5 

OdB 

Table  37  Costas  signa 

s  to  be  processed  by  QMFB  Tree. 

From  the  list  of  signals  already  processed  by  the  QMFB  (in  Table  37),  only  one  of 
them  is  shown  next.  The  rest  of  the  signals  can  be  seen  in  a  technical  report  that  will  be 
published  soon  [21]. 
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a.  FSK_PSK_C_1_1 5_ll_5_s 

This  FSK/PSK  Costas  code  signal  was  generated  with  the  parameters  de¬ 
scribed  in  Table  38.  The  frequency  sequence  is  the  sequence  number  1  (4-7-1-6-5-2-3 
kHz),  1 1  bits  Barker  code  with  5  cycles  per  bit.The  bandwidth  used  in  each  one  of  the 
carry  frequencies  is  kept  constant;  therefore,  the  number  of  cycle  per  bit  must  vary  as 
shown  in  Equations  (3.32)  and  (3.34).  As  a  consequence  of  keeping  the  same  bandwidth 
for  each  frequency,  the  code  period  and  the  time  spent  by  the  signal  in  each  frequency  is 
constant  as  demonstrated  in  (3.33)and  (3.35). 

For  example  the  bandwidth  for  the  signal  with  carrier  frequency  1  kHz  is 


BW, 


1  KHz 


Carrier  frequency 
Cycles  per  bit 


fc_ 

cpp 


1000  Hz 
5 


=  200  Hz 


(3.32) 


and  the  Barker  code  period  is 


T  (Cycles  per  bit)(No  of  bits) 

A  C 

Carrier  frequency 


(5X11) 

1000 


=  55  ms. 


(3.33) 


The  bandwidth  for  the  signal  with  carrier  frequency  7  kHz  is 


bw7KHz  = 


Carrier  Frequency  _  fc  _  7000  Hz 
Cycles  per  bit  cpp  35 


200  Hz 


(3.34) 


and  the  Barker  code  period  is 


T  _  (Cycles  per  bit)(No  of  bits) 

A  c  \ 

Carrier  frequency 


(35X11)=55 

7000 


ms. 


(3.35) 


It  can  be  inferred  that  since  the  code  period  is  55  ms,  the  transmission  time 
in  each  frequency  of  the  sequence  is  five  times  the  code  period  (due  to  the  signal  was 
generated  with  5  periods  Barker  code  for  each  frequency  [15]),  giving  a  transmission 
time  in  each  frequency  of  275  ms.  Therefore,  since  the  Costas  sequence  has  7  different 
frequencies,  the  Costas  code  period  is  7  times  the  transmission  time  per  frequency,  result¬ 
ing  in  a  total  of  1.925  s. 

Now,  and  according  with  the  description  of  the  code  given  in  Section  C.l 
in  Chapter  III,  the  signal  is  processed  as  an  input  to  the  QMFB  tree  to  obtain  the  different 
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output  layers.  In  this  manner,  Figure  72  was  created  which  show  the  output  matrix  from 
layer  9  with  a  frequency  resolution  of  14.734  Hz.  A  zoom  of  the  same  layer  output  matrix 
is  shown  in  Figure  73,  there  the  values  of  the  Costas  frequency  sequence  are  identified,  as 
well  as  the  Costas  code  period  and  bandwidth,  giving  the  same  values  than  the  original 
signal. 


'TimF! .  f  rnqi.encv  Flam  LsjntrB' 
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Figure  72  Output  matrix  from  layer  9  of  FSK_PSK_C_1_15_1  l_5_s  (colorplot). 
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Figure  73  Zoom  of  output  matrix  from  layer  9  of  FSK_PSK_C_1_15_1  l_5_s  (color- 
plot). 

Examining  a  previous  output  matrix  layer  it  is  possible  to  achieve  a  better 
resolution  in  time  and  to  see  more  details  of  the  signals  in  each  frequency  of  the  se¬ 
quence.  Figure  74  shows  the  matrix  output  layer  8  in  a  contourplot  marking  the  frequency 
sequence  to  be  analyzed  (4  kHz).  Figure  75  shows  a  zoom  of  the  same  layer  focusing  in 
the  4  kHz  frequency  of  the  sequence  for  the  second  Costas  code  period.  There  the  band¬ 
width  (200  kHz),  transmission  time  in  the  frequency  (275  ms),  and  the  Barker  code  pe¬ 
riod  (55  ms)  can  be  extracted.  The  1 1  Barker  bits  are  not  visible  in  the  plot  because  the 
resolution  is  not  as  high  as  desirable.  Also,  in  the  plot  appears  a  modulation  of  5  periods 
in  the  same  carrier  frequency  indicating  that  the  Barker  code  is  repeated  five  times  been 
this  a  setting  in  the  Low  Probability  of  Intercept  Generator  Program  (LPIG  [15]). 
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Figure  74  Output  matrix  from  layer  8  of  FSK_PSK_C_1_15_1  l_5_s  (contourplot). 
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Figure  75  Zoom  of  output  matrix  from  layer  8  of  FSK_PSK_C_1_15_1  l_5_s  (contour- 
plot). 
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Table  38  presents  a  summary  of  the  signal  processing  showing  a  compari¬ 


son  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


FSK  PSK  C  1  15  11  5  s 
Parameters 

Generation 

Detection 

Comment 

Costas  Sequence 

4-7-1-6-5-2-3  kHz 

4-7-1-6-5-2-3  kHz 

Sampling  Frequency 

15  kHz 

15  kHz 

Given 

SNR 

Only  Signal 

- 

Cycles  per  bit 

5 

- 

Barker  Code  Period 

55  ms 

55  ms 

Number  of  Bits  per  Barker 
Code 

11 

" 

Transmission  time  per  Fre¬ 
quency 

275  ms 

275  ms 

Bandwidth  per  Carrier  fre¬ 
quency  (Barker) 

200  Hz 

200  Hz 

Costas  Code  Period 

1.925  s 

1.925  s 

Costas  Bandwidth 

6200  Hz 

6306  Hz 

106  Hz  error 

Table  38  Signal  processing  summary  for  FSK_PSK_C_1_15_1  l_5_s. 


b.  FSK_PSK_C_1_15_11_5_0 

This  FSK/PSK  Costas  code  signal  was  generated  with  the  parameters  de¬ 
scribed  in  Table  39.  It  must  be  noted  that  this  signal  has  almost  the  same  parameters  as 
the  previous,  the  only  difference  is  the  SNR  that  now  is  0  dB.  Figure  76  was  obtained, 
which  provides  a  good  time-frequency  description  of  the  evaluated  signal  from  the  output 
matrix  at  layer  9  when  conducting  the  processing  with  the  QMFB  tree  with  a  frequency 
resolution  of  14.734.  A  zoom  of  the  same  layer  output  matrix  is  shown  in  Figure  77;  the 
values  of  the  Costas  frequency  sequence  were  identified,  as  well  as  the  Costas  code  pe¬ 
riod  and  bandwidth,  giving  the  same  values  than  the  original  signal. 
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Figure  76  Output  matrix  from  layer  9  of  FSK_PSK_C_1_15_1 15  0  (colorplot). 
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Figure  77  Zoom  output  matrix  from  layer  9  of  FSK_PSK_C_1_15_1 15  0  (colorplot). 
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Examining  a  previous  output  matrix  layer  it  is  possible  to  achieve  a  better 
resolution  in  time  and  see  more  details  of  the  signals  in  each  frequency  of  the  sequence. 
Figure  78  shows  the  matrix  output  layer  8  in  a  contourplot  marking  the  frequency  se¬ 
quence  to  be  analyzed  (4  kHz).  Figure  79  shows  a  zoom  of  the  same  layer  focusing  in  the 
4  kHz  frequency  of  the  sequence  for  the  second  Costas  code  period.  The  bandwidth  (200 
kHz),  transmission  time  in  the  frequency  (275  ms),  and  the  Barker  code  period  (55  ms) 
can  be  extracted.  The  1 1  Barker  bits  are  not  visible  in  the  plot  because  the  resolution  is 
not  as  high  as  desirable.  Also,  in  the  plot  appears  a  modulation  of  5  periods  in  the  same 
carrier  frequency  indicating  that  the  Barker  code  is  repeated  five  times,  been  a  setting  in 
the  FPIG  [15]. 
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Figure  78  Output  matrix  from  layer  8  of  FSK_PSK_C_1_15_1 15  0  (contourplot). 
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Figure  79  Zoom  of  output  matrix  from  layer  8  of  FSK_PSK_C_1_15_1 15  0  (contour- 
plot). 


Table  39  presents  a  summary  of  the  signal  processing  showing  a  compari¬ 


son  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB  tree. 


FSK  PSK  C  1  15  11  5  0 
Parameters 

Generation 

Detection 

Comment 

Costas  Sequence 

4-7-1-6-5-2-3  kHz 

4-7-1-6-5-2-3  kHz 

Sampling  Frequency 

15  kHz 

15  kHz 

Given 

SNR 

OdB 

- 

Cycles  per  bit 

5 

- 

Barker  Code  Period 

55  ms 

55  ms 

Number  of  Bits  per  Barker 
Code 

11 

" 

Transmission  time  per  Fre¬ 
quency 

275  ms 

275  ms 

Bandwidth  per  Carrier  fre¬ 
quency  (Barker) 

200  Hz 

200  Hz 

Costas  Code  Period 

1.925  s 

1.925  s 

Costas  Bandwidth 

6200  Hz 

6306  Hz 

106  Hz  error 

Table  39  Signal  processing  summary  for  FSK_PSK_C_1_15_1 15  0. 
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K.  FSK/PSK  COMBINED  WITH  TARGET-MATCHED  FREQUENCY 
HOPPING 


1.  Brief  Description 

Instead  of  spreading  the  energy  of  the  signal  equally  over  a  broad  bandwidth,  this 
type  of  technique  concentrates  the  signal  energy  in  specific  spectral  locations  of  most  im¬ 
portance  for  the  radar  and  its  typical  targets,  within  the  broad-spectrum  bandwidth.  The 
produced  signals  have  a  pulse  compression  characteristic  and  therefore  they  can  achieve 
low  probability  of  intercept. 

The  implementation  starts  with  a  simulated  target  time  radar  response.  This  data 
is  then  Fourier  transformed  and  the  correspondent  frequencies  and  initial  phases  are  cal¬ 
culated  using  built-in  Matlab  functions  in  LPIG  [17].  A  random  selection  process  chooses 
each  frequency  (between  64  different  frequencies)  with  a  probability  distribution  function 
defined  by  the  spectral  characteristics  of  the  target  of  interest  obtained  from  the  FFT,  so 
that  frequencies  at  the  spectral  peaks  of  the  target  (highest  magnitudes)  are  transmitted 
more  often.  Each  frequency  hop,  transmitted  during  a  specific  period  of  time,  is  also 
modulated  in  phase,  having  its  initial  phase  value  modified  by  a  pseudo-random  spread¬ 
ing  phase  sequence  code  of  values  equally  likely  to  be  zero  or  n  radians. 

The  matched  FSK/PSK  radar  will  then  use  a  correlation  receiver  with  a  phase 
mismatched  reference  signal  instead  of  a  perfectly  phase  matched  reference.  This  allows 
the  radar  to  generate  signals  that  can  match  a  target’s  spectral  response  in  both  magnitude 
and  phase. 

Figure  80  illustrates  the  block  diagram  for  the  generation  of  FSK/PSK.  This  dia¬ 
gram  and  a  deep  description  of  the  signal  can  be  found  in  [17]  and  [18].  In  addition  Fig¬ 
ures  81,  82,  and  83  show  the  FSK/PSK  target  simulated  response,  the  frequency 
components  and  frequency  components  histogram  with  number  of  occurrences  per  fre¬ 
quency  for  256  frequency  hops,  respectively. 
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Figure  80  Block  diagram  of  the  implementation  of  the  FSK/PSK  Target  matched  wave¬ 
form  starting  from  the  target  radar  response. 


Figure  8 1  FSK/PSK  Target  simulated  response. 
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Figure  82  FSK/PSK  Target  frequency  components  and  frequency  probability  distribu¬ 
tion. 


Figure  83  FSK/PSK  Target  frequency  components  histogram  with  number  of  occur¬ 
rences  per  frequency  for  256  frequency  hops. 

2.  Processing  FSK/PSK  Target  Code  Signals  with  QMFB  Tree 

The  FSK/PSK  Target  code  signals  to  be  worked  with  are  given  in  Table  40.  All 
the  signals  were  generated  with  5  or  10  different  phases  and  128  or  256  random  frequen¬ 
cies  (from  the  original  64  frequencies),  the  sampling  frequency  was  15000  Hz,  the  num- 
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ber  of  cycles  per  phase  and  the  number  of  phases  have  the  same  value  in  the  signal  (since 
the  cycle  per  phase  is  1,  the  total  amount  of  cycles  is  directly  related  with  the  number  of 
phases),  and  the  SNR  are  variable  parameters. 


FSK/PSK  Target 

Number  of 
random  hops 
in  frequency 

Number  of 
cycles  and 
phases 

SNR 

1 

FSK  PSK  T  15  128  5  s 

128 

5 

Signal  Only 

2 

FSK  PSK  T  15  128  5  0 

128 

5 

OdB 

3 

FSK  PSK  T  15  128  10  s 

128 

10 

Signal  Only 

4 

FSK  PSK  T  15  128  10  0 

128 

10 

OdB 

5 

FSK  PSK  T  15  256  5  s 

256 

5 

Signal  Only 

6 

FSK  PSK  T  15  256  5  0 

256 

5 

OdB 

7 

FSK  PSK  T  15  256  10  s 

256 

10 

Signal  Only 

8 

FSK  PSK  T  15  256  10  0 

256 

10 

OdB 

Tal 

ble  40  FSK/PSK  Target  signals  to  be  processed  by  QMFB  Tree. 

From  the  list  of  signals  already  processed  by  the  QMFB  (in  Table  40),  only  one  of 
them  is  shown  next.  The  rest  of  the  signals  can  be  seen  in  a  technical  report  that  will  be 
published  soon  [21]. 


a.  FSK_PSK_T_15_128_1 0_s 

This  FSK/PSK  Target  code  signal  was  generated  with  the  parameters  de¬ 
scribed  in  Table  41,  with  10  cycles,  10  phases,  and  128  random  frequencies.  Now,  and 
according  with  the  description  of  the  code  given  in  Section  C.l  in  Chapter  III,  the  signal 
is  processed  as  an  input  to  the  QMFB  Tree  to  obtain  the  different  output  layers.  In  this 
manner  Figure  84  and  Figure  85  were  created  which  show  the  output  matrix  from  layer  7 
in  color  and  contour  plot,  respectively.  As  a  result  of  the  particularly  efficient  LPI  charac¬ 
teristics  of  this  modulation,  there  is  not  much  information  that  can  be  extracted  from  the 
plots  in  addition  to  the  bandwidth  and  the  concentration  of  the  power  in  frequency.  The 
random  presentation  of  frequencies  in  the  signal  gives  the  impression  of  containing  only 
noise.  That  way  it  can  be  observed  a  kind  of  noise  in  the  signal  and  the  power  is  mostly 
concentrated  between  3  and  6  kHz  as  expected. 
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Figure  84  Output  matrix  from  layer  7  of  FSK_PSK_T_15_128_10_s  (colorplot). 
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Figure  85  Output  matrix  from  layer  7  of  FSK_PSK_T_15_128_10_s  (contourplot). 
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Table  41  presents  a  summary  of  the  signal  processing  demonstrating  a 
comparison  of  the  real  signal  parameters  versus  the  parameters  extracted  by  the  QMFB 
tree. 


FSK_PSK_T_15_128_10_s 

Parameters 

Generation 

Detection 

Comment 

Sampling  Frequency 

15  kHz 

15  kHz 

Given 

Number  of  Phases  and  cy¬ 
cles 

10 

" 

Number  of  Frequencies 

64 

- 

Number  of  Hops 

128 

- 

Bandwidth 

3000  Hz 

3000  Hz 

SNR 

Only  Signal 

- 

Table  41  Signal  processing  summary  for  FSK_PSK_T_15_128_10_s. 


b.  FSK_PSK_T_1 5_1 28_1 0_0 

This  FSK/PSK  Target  code  signal  was  generated  with  the  parameters  de¬ 
scribed  in  Table  42.  This  signal  has  almost  the  same  parameters  as  the  previous  signal. 
The  only  difference  is  the  SNR  that  now  is  0  dB.  As  the  signal  is  too  noisy,  the  noise 
hides  the  signal  and  spreads  the  bandwidth  making  it  extremely  difficult  to  detect  the 
FSK_PSK_T_15_128_10_0  signal  parameters  and  extract  some  features.  It  is  even  harder 
to  try  to  extract  the  code  period. 

In  Table  42  a  summary  of  the  signal  processing  is  shown,  demonstrating  a 
comparison  of  the  real  signal  parameters  versus  the  parameters  detected  by  the  QMFB 
tree. 


FSK_PSK_T_15_128_10_0 

Parameters 

Generation 

Detection 

Comment 

Sampling  Frequency 

15  kHz 

- 

Given 

Number  of  Phases  and  cy¬ 
cles 

10 

" 

Number  of  Frequencies 

64 

- 

Number  of  Hops 

128 

- 

Bandwidth 

3000  Hz 

- 

SNR 

0  dB 

- 

Table  42  Signal  processing  summary  for  FSK_PSK_T_15_128_10_0. 
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L.  ANALYSIS  AND  SUMMARY  OF  THE  DIFFERENT  SIGNALS 
PROCESSING  RESULTS 


After  been  processing  all  the  LPI  signals  listed  in  Appendix  B  with  the  QMFB 
tree  and  analyzing  the  results  it  is  possible  to  summarize  the  performance  of  the  QMFB 
in  detecting,  displaying,  and  extracting  features  from  the  signals. 

1.  Binary  Phase  Shift  Keying 

As  was  expected,  the  QMFB  tree  revealed  a  very  good  performance  in  detecting 
the  carrier  frequency  of  the  BPSK  signals.  The  processing  had  no  major  difficulty  in  de¬ 
tecting  the  carrier  frequency  of  signals  with  0  or  -6  dB  of  noise.  With  respect  to  the 
bandwith  the  QMFB  tree  shows  that  the  extraction  of  this  feature  presents  the  same  per¬ 
formance  as  the  carrier  frequency  detection.  The  cycles  per  phase  were  not  detected  by 
the  processing  due  to  the  resolution  of  the  diagrams  not  being  high  enough  to  do  this.  Ex¬ 
tracting  the  Code  period  the  processing  showed  a  good  performance,  but  as  for  signals 
with  high  levels  of  noise  (-6  dB)  the  resolution  did  not  allow  extracting  this  feature  (no 
yellow  bar  for  this  feature  in  the  bar  plot  below).  Figure  86  shows  a  summary  chart  of  the 
performance  of  the  QMFB  to  detect  carrier  frequency,  bandwidth,  code  period  and  bits  in 
Barker  code  in  the  BPSK  signals  analyzed  previously.  This  figure  presents  a  comparison 
of  the  parameters  in  three  different  environments:  signal  only  (blue),  SNR=  0  dB  (red) 
and  SNR=  -6  dB  (yellow).  The  percentage  in  the  chart  describes  an  average  of  how  close 
the  extracted  values  are  to  the  theoretical  values. 
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Carrier  freq. (Hz) 

Bandwidth  (Hz) 

Code  period  (ms) 

Bits/code 

Only  signal 
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106% 
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0  dB 
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0% 

(-)  6  dB 
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0% 

0% 
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Figure  86  Performance  of  the  QMFB  processing  detecting  BPSK  signals. 


2.  Frequency  Modulation  Continuous  Wave  (FMCW) 

For  this  signal  the  QMFB  tree  show  a  very  good  resolution  in  detecting  the  carrier 
frequency  under  any  SNR  condition.  The  modulation  bandwidth  extraction  gives  an  error 
incrementing  the  real  value  as  the  SNR  goes  to  negative.  The  modulation  period  extrac¬ 
tion  shows  a  better  detection  for  signals  generated  with  longer  modulation  period  (30  ms) 
than  shorter  ones  (20  ms)  and  higher  modulation  bandwidth  (500  Hz),  but  always  has  a 
problem  extracting  this  feature  from  signals  with  SNR=  -6  dB.  Figure  87  shows  a  sum¬ 
mary  chart  of  the  performance  of  the  QMFB  to  detect  carrier  frequency,  modulation 
bandwidth  and  modulation  period  in  the  FMCW  signals  analyzed  previously.  This  figure 
presents  a  comparison  of  the  parameters  in  three  different  environments:  signal  only 
(blue),  SNR=  0  dB  (red)  and  SNR=  -6  dB  (yellow).  The  percentage  in  the  chart  describes 
an  average  of  how  close  the  extracted  values  are  from  the  theoretical  values. 
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Figure  87  Performance  of  the  QMFB  processing  detecting  FMCW  signals 


3.  Frank  Code 

The  processing  of  this  kind  of  signal  with  the  QMFB  tree  gives  good  resolution 
extracting  the  carrier  frequency,  code  period,  and  the  number  of  phases  for  signals  with¬ 
out  noise,  but  the  extraction  of  these  parameters  for  noisy  signals  is  difficult,  specially 
obtaining  the  number  of  phases.  Extracting  the  bandwidth  of  the  signal  was  more  accu¬ 
rate  giving  good  results  in  most  of  the  processing,  except  for  the  -6  dB  condition.  Figure 
88  shows  a  summary  chart  of  the  performance  of  the  QMFB  to  detect  carrier  frequency, 
bandwidth,  code  period  and  number  of  phases  in  Frank-coded  signals  analyzed  previ¬ 
ously.  This  figure  presents  a  comparison  of  the  parameters  in  three  different  environ¬ 
ments:  signal  only  (blue),  SNR=  0  dB  (red)  and  SNR=  -6  dB  (yellow).  The  percentage  in 
the  chart  describes  an  average  of  how  close  the  extracted  values  are  from  the  theoretical 
values. 
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Carrier  freq. (Hz) 

Bandwidth  (Hz) 

Code  period  (ms) 

Phases 
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102% 

1 1 6% 

1 04% 

100% 

0  dB 

102% 

107% 

1 04% 

0% 

(-)  6  dB 

75% 

77% 

78% 

0% 

Frank  code  effectiveness  on  detection 


150% 
130%  - 


1 1 0% 


§>  90% -  -  —  - 

2 

|  70% -  JT-J  -  -  -  -  ----- 

CL 

50% -  -  -  - 

30% -  -  -  -  -  - 

10%  -I - - ^ - - T - - Ml - - T - - - - - T - - 

Carrier  freq. (Hz)  Bandw  idth  (Hz)  Code  period  (ms)  Phases 

Parameters 


□  Only  signal 
■  0  dB 

□  (-)  6  dB 


Figure  88  Performance  of  the  QMFB  processing  detecting  Frank-coded  signals. 


4.  PI  Polyphase  Code 

Processing  this  signals,  the  QMFB  detected  the  carrier  frequency,  the  bandwidth, 
and  the  code  period  with  very  good  resolution  for  the  case  of  signal  only  and  SNR=0  dB, 
but  under  SNR=  -6  dB,  the  efficiency  of  the  QMFB  goes  down  in  a  50%  approximately. 
The  number  of  phases  can  be  extracted  just  in  signals  without  noise.  Figure  89  shows  a 
summary  chart  of  the  performance  of  the  QMFB  to  detect  carrier  frequency,  bandwidth, 
code  period  and  number  of  phases  in  PI -coded  signals  analyzed  previously.  This  figure 
presents  a  comparison  of  the  parameters  in  three  different  environments:  signal  only 
(blue),  SNR=  0  dB  (red)  and  SNR=  -6  dB  (yellow).  The  percentage  in  the  chart  describes 
an  average  of  how  close  the  extracted  values  are  from  the  theoretical  values.  After 
processing  the  signals,  the  plots  provide  the  values  for  carrier  frequency,  bandwidth  and 
code  period  in  any  environment  analyzed. 
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Figure  89  Performance  of  the  QMFB  processing  detecting  PI -coded  signals. 


5.  P2  Polyphase  Code 

Processing  this  signals,  the  QMFB  detect  the  carrier  frequency,  the  bandwidth, 
and  the  code  period  with  very  good  resolution  for  the  case  of  signal  only  and  SNR=  0  dB, 
but  under  SNR=  -6  dB,  the  efficiency  of  the  QMFB  goes  down  in  a  50%  approximately. 
The  number  of  phases  can  be  extracted  just  in  signals  without  noise.  Figure  90  shows  a 
summary  chart  of  the  performance  of  the  QMFB  to  detect  carrier  frequency,  bandwidth, 
code  period  and  number  of  phases  in  P2-coded  signals  analyzed  previously.  This  figure 
presents  a  comparison  of  the  parameters  in  three  different  environments:  signal  only 
(blue),  SNR=  0  dB  (red)  and  SNR=  -6  dB  (yellow).  The  percentage  in  the  chart  describes 
an  average  of  how  close  the  extracted  values  are  from  the  theoretical  values.  After  proc¬ 
essing  the  signals,  the  plots  provide  very  accurate  values  for  carrier  frequency,  and  band¬ 
width,  but  not  code  period. 
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Figure  90  Performance  of  the  QMFB  processing  detecting  P2-coded  signals 


6.  P3  Polyphase  Code 

Processing  this  signals  with  the  QMFB  gave  almost  the  same  results  as  the  P2  and 
PI  signals  showing  an  excellent  job  detecting  the  carrier  frequency,  bandwidth,  and  code 
period  for  signals  without  noise  or  SNR=  0  dB,  but  a  very  bad  under  SNR=  -6  dB.  The 
number  of  phases  can  be  extracted  just  in  signals  without  noise.  Figure  91  shows  a  sum¬ 
mary  chart  of  the  performance  of  the  QMFB  to  detect  carrier  frequency,  bandwidth,  code 
period  and  number  of  phases  in  P3-coded  signals  analyzed  previously.  This  figure  pre¬ 
sents  a  comparison  of  the  parameters  in  three  different  environments:  signal  only  (blue), 
SNR=  0  dB  (red)  and  SNR=  -6  dB  (yellow).  The  percentage  in  the  chart  describes  an  av¬ 
erage  of  how  close  the  extracted  values  are  from  the  theoretical  values.  After  processing 
the  signals,  the  plots  provide  very  accurate  values  for  carrier  frequency,  bandwidth  and 
code  period  in  any  environment  analyzed  with  a  lower  performance  when  SNR  is  less 
than  -6  dB.  The  number  of  phases  can  be  only  extracted  when  noise  is  not  presented  and 
a  number  of  cycles  per  phase  is  greater  than  5. 
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Figure  91  Performance  of  the  QMFB  processing  detecting  P3-coded  signals. 


7.  P4  Polyphase  Code 

This  signals  demonstrated  that  the  QMFB  gave  almost  the  same  results  as  the 
other  polyphase  code  signals  showing  an  excellent  job  detecting  the  carrier  frequency, 
bandwidth,  and  code  period  for  signals  without  noise  or  SNR=  0  dB,  but  a  very  bad  under 
SNR=  -6  dB.  The  number  of  phases  can  be  extracted  just  in  signals  without  noise.  Figure 
92  shows  a  summary  chart  of  the  performance  of  the  QMFB  to  detect  carrier  frequency, 
bandwidth,  code  period  and  number  of  phases  in  P4-coded  signals  analyzed  previously. 
This  figure  presents  a  comparison  of  the  parameters  in  three  different  environments:  sig¬ 
nal  only  (blue),  SNR=0  dB  (red)  and  SNR=  -6  dB  (yellow).  The  percentage  in  the  chart 
describes  an  average  of  how  close  the  extracted  values  are  from  the  theoretical  values. 
After  processing  the  signals,  the  plots  provide  very  accurate  values  for  carrier  frequency, 
bandwidth  and  code  period  in  any  environment  analyzed  with  a  lower  performance  when 
SNR  is  less  than  -6  dB.  The  number  of  phases  can  be  only  extracted  when  noise  is  not 
presented  and  a  number  of  cycles  per  phase  is  greater  than  5. 
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Figure  92  Performance  of  the  QMFB  processing  detecting  P4-coded  signals. 


8.  Costas  Code 

This  specific  kind  of  signal  when  processed  by  the  QMFB  tree  gave  excellent  fre¬ 
quency  sequence  detection  with  the  best  resolution  measuring  the  transmission  time  per 
frequency.  The  code  period  was  detected  with  a  very  good  resolution  for  all  the  signals 
and  the  same  happened  for  the  bandwidth.  Figure  93  shows  a  summary  chart  of  the  per¬ 
formance  of  the  QMFB  to  detect  carrier  frequency,  bandwidth,  and  code  period  in  Co- 
stas-coded  signals.  This  figure  presents  a  comparison  of  the  parameters  in  three  different 
environments:  signal  only  (blue),  SNR=  0  dB  (red)  and  SNR=  -6  dB  (yellow).  The  per¬ 
centage  in  the  chart  describes  an  average  of  how  close  the  extracted  values  are  from  the 
theoretical  values.  After  processing  the  signals,  the  plots  provide  very  accurate  values  for 
carrier  frequency,  bandwidth  and  code  period  in  any  environment  analyzed. 
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Figure  93  Performance  of  the  signal  processing  detecting  Costas-coded  signals. 


9.  Frequency  Shift  Keying/Phase  Shift  Keying  Combined  with  Costas 
Code  (FSK/PSK  Costas) 

This  very  specific  kind  of  signal  has  very  important  characteristics  and  applica¬ 
tions  as  an  LPI  radar  signal  and  the  QMFB  tree  had  a  very  good  performance  analyzing 
it.  The  frequency  sequence  was  detected  in  every  case  with  excellent  resolution.  The 
same  situation  happened  in  the  detection  of  the  Barker  code  period  and  the  transmission 
time  per  frequency.  The  bandwidth  per  carrier  frequency  was  detected  with  a  good  reso¬ 
lution,  harder  to  see  than  in  other  signals,  but  clear  enough  to  be  defined.  Figure  94 
shows  a  summary  chart  of  the  performance  of  the  QMFB  to  detect  carrier  frequency, 
bandwidth  and  code  period  in  FSK/PSK  Costas-coded  signals  analyzed  previously.  This 
figure  presents  a  comparison  of  the  parameters  in  two  different  environments:  signal  only 
(blue)  and  SNR=  0  dB  (red),  hence  there  are  no  yellow  bars.  The  percentage  in  the  chart 
describes  an  average  of  how  close  the  extracted  values  are  from  the  theoretical  values. 
After  processing  the  signals,  the  plots  provide  very  precise  values  for  carrier  frequency, 
bandwidth  and  code  period  in  any  environment  analyzed  except  for  SNR  lower  than  0  dB 
where  none  of  the  parameters  can  be  extracted. 
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Figure  94  Performance  of  the  QMFB  processing  detecting  Costas-coded  signals. 

10.  Frequency  Shift  Keying/Phase  Shift  Keying  Combined  with  Target- 
Matched  Frequency  Hopping  (FSK/PSK  Target) 

Processing  this  signals  with  the  QMFB  tree  revealed  that  it  was  very  hard  to  ex¬ 
tract  features  other  than  the  bandwidth.  But  since  these  signals  were  generated  with  128 
or  256  random  jumps  in  frequency  between  64  original  and  different  ones  from  a  target 
radar  response,  it  is  extremely  hard  to  determine  the  frequency  sequence.  Therefore,  the 
plots  generated  by  the  processing  looks  like  a  signal  with  a  strong  concentration  of  noise 
in  a  specific  bandwidth,  making  it  really  hard  to  identify  it  as  a  FSK/PSK  Target  signal. 
Figure  95  shows  a  summary  chart  of  the  performance  of  the  QMFB  to  detect  carrier  fre¬ 
quency,  bandwidth  and  code  period  in  FSK/PSK  Target  signals  analyzed  previously.  This 
figure  presents  a  comparison  of  the  parameters  in  three  different  environments:  signal 
only  (blue),  SNR=  0  dB  (red)  and  SNR=  -6  dB  (yellow).  The  percentage  in  the  chart  de¬ 
scribes  an  average  of  how  close  the  extracted  values  are  from  the  theoretical  values.  After 
processing  the  signals,  the  plots  provide  very  precise  values  for  bandwidth  for  signal-only 
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condition.  The  sequence  and  code  period  in  any  environment  analyzed  cannot  be  ex¬ 
tracted. 
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Figure  95  Performance  of  the  QMFB  processing  detecting  FSK/PSK  target  signals. 


M.  COMPARISON  OF  DIFFERENT  POLYPHASE-CODED  SIGNALS 

The  objective  of  comparing  different  polyphase-coded  signals  is  to  demonstrate 
the  effectiveness  of  the  proposed  signal  processing  to  discriminate  among  different  poly¬ 
phase  modulations.  One  signal  per  each  modulation  was  selected,  having  the  same  carrier 
frequency  (1  kHz)),  sampling  frequency  (7  kHz),  number  of  phases  (16)  and  cycles  per 
phase  (5).  The  most  important  deviation  among  the  signal  is  related  to  the  phase  shift. 
The  increment  between  phases  produces  diverse  performance  in  the  modulations,  facili¬ 
tating  an  accurate  identification  of  each  polyphase  modulation. 
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Table  43  shows  each  one  of  the  analyzed  signals  and  their  most  important  charac¬ 
teristics.  In  addition,  this  table  presents  the  resulting  phases  for  the  generation  of  16 
phases  (N=  4  or  N2=  16).  The  resulting  figures  after  QMFB  are  provided  along  with  the 
phase  shift  plots.  The  analysis  of  these  plots  reveals  a  close  relationship  where  the  con¬ 
ducts  of  the  phases  approximate  to  the  frequency-time  description  of  the  signal. 


Signal 

Phase  shift  values 

Comment 

Frank 

0  0  0  0 

0  1.5708  3.1416  4.7124 

0  3.1416  6.2832  9.4248 

0  4.7124  9.4248  14.1372 

Phase  changes  between  adja¬ 
cent  codes  are  the  smallest. 

PI 

0  -3.1416  6.2832  28.2743 

-2.3562  -3.9270  7.0686  30.6305 
-4.7124  -4.7124  7.8540  32.9867 
-7.0686  -5.4978  8.6394 

Lowest  code  increments  from 
one  code  element  to  code 
element  in  the  center  of  the 
waveform 

P2 

3.5343  1.1781  -1.1781  -3.5343 

1.1781  0.3927  -0.3927  -1.1781 

-1.1781  -0.3927  0.3927  1.1781 
-3.5343  -1.1781  1.1781  3.5343 

Symmetric  at  center  fre¬ 
quency 

P3 

0  0.0123  0.0491  0.1104  0.1963  0.3068 

0.4418  0.6013  0.7854  0.9940  1.2272  1.4849  1.7671 

2.0739  2.4053  2.7612  (one  row) 

Largest  code  element  to  code 
element  are  on  middle  of  the 
P3  code 

P4 

0  -2.9452  -5.4978  -7.6576  -9.4248  -10.7992 

-11.7810-12.3700  -12.5664-12.3700  -11.7810  -10.7992 
-9.4248  -7.6576  -5.4978  -2.9452  (one  row) 

Largest  code  element  to  code 
element  are  on  the  two  ends 
of  the  P4  code 

- - • - . - . - - 2 - 

Table  43  Different  Polyphase-coded  signals  and  differences  for  N  =16. 


1.  Frank  Code 

The  actual  Frank-coded  signal  consists  of  a  carrier  (1  kHz),  the  place  of  which  is 
modulated  according  to  the  indicated  baseband  waveform  sequence.  For  each  frequency 
or  section  of  the  step  chirp,  a  phase  group  consisting  of  N  phase  samples  is  obtained  and 
the  total  number  of  code  phases  is  N2.  Note  that  the  phase  increments  within  the  four 
phase  groups  are  0°,  90°,  180°,  and  270°.  However,  the  phases  of  the  last  group  are  am¬ 
biguous  (>180°)  and  appear  as  -90°  phase  steps.  The  last  phase  group,  because  of  the 
ambiguity,  appears  to  complete  one  360°  (Figure  96). 
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Figure  96  Frank-coded  signal  resulting  plot  after  QMFB  and  phase  shift. 

2.  PI 

This  code  has  the  lowest  code  increments  from  one  code  element  to  code  element 
in  the  center  of  the  waveform  as  shown  in  Figure  97.  For  an  odd  number  of  N,  the  result¬ 
ing  phases  are  the  same  as  the  Frank  code  except  the  phase  groups  are  rearranged.  For  N 
even,  PI  code  has  the  same  phase  increments,  within  each  phase  group,  as  the  P2  except 
the  starting  phases  are  different. 
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Figure  97  PI -coded  signal  resulting  plot  after  QMFB  and  phase  shift. 


3.  P2 

This  code  is  valid  for  even  numbers  of  phases,  and  each  group  of  the  code  is 
symmetric  about  zero  phase  as  shown  in  Figure  98. 
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Figure  98  P2-coded  signal  resulting  plot  after  QMFB  and  phase  shift. 


4.  P3 

P3  only  differs  from  Frank  code  by  180°  phase  shifts  every  N  code  elements 
(one  frequency  group)  and  by  added  phase  increments  that  repeats  every  N  samples 
(every  frequency  group).  The  largest  phase  increments  from  code  element  to  code  ele¬ 
ment  are  in  the  middle  of  the  P3  code  as  presented  in  Figure  99. 
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Figure  99  P3 -coded  signal  resulting  plot  after  QMFB  and  phase  shift. 


5.  P4 

P4  code  is  very  similar  to  PI  code  except  that  the  phase  samples  are  those  of  a 
sampled  chirped  waveform  rather  than  step-chirp  waveform.  It  is  noted  that  the  largest 
phase  increment  s  from  code  element  are  on  the  two  ends  of  the  P4  code.  Figure  100  de¬ 
scribes  the  resulting  signal  and  its  phase  shift. 
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Figure  100  P4-coded  signal  resulting  plot  after  QMFB  and  phase  shift. 


This  chapter  discussed  the  results  of  processing  the  LPI  signals  from  the  LPIG 
code  as  described  in  [15].  It  finishes  with  an  analysis  and  comparison  of  the  different  sig¬ 
nals  processing  results,  where  the  performances  of  the  QMFB  tree  are  specified  for  every 
class  of  signal  depending  on  the  results  obtained.  The  next  chapter  gives  conclusions  and 
recommendations  generated  from  this  research. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  goal  of  this  research  is  to  use  a  method  to  decompose  waveforms,  using  or¬ 
thogonal  basis  functions  and  a  quadrature  mirror  filter  bank  (QMFB)  tree,  and  to  extract 
detailed  information  about  embedded  LPI  (spread  spectrum)  radar  signals. 

A  list  of  different  and  well-known  LPI  radar  signals  were  processed  with  the 
QMFB  tree  and  the  output  layers  were  displayed  and  analyzed  in  a  feature  extraction  pro¬ 
cedure.  Specifically,  estimating  the  signal  parameters  such  as  bandwidth,  center  fre¬ 
quency,  duration  and  position  in  time  of  the  waveform  were  accomplished.  The 
processing  started  by  analyzing  a  tone  signal  and  binary  phase  shift  keying  signals  to 
show  how  the  processing  works,  what  kind  of  results  were  expected,  and  how  to  extract 
the  parameters  from  the  outputs  of  the  QMFB  tree. 

Polyphases  signals  were  processed  obtaining  very  effective  feature  extraction  re¬ 
sults,  where  different  parameters  such  as  carrier  frequency,  bandwidth,  and  code  period 
were  extracted  with  good  resolution  for  the  case  of  signal  only  and  SNR=  0  dB.  In  the 
case  of  SNR=  -6  dB,  the  extraction  of  this  parameters  was  difficult  due  to  the  noisy  con¬ 
dition  of  the  signals,  but  still  gave  results  close  enough  to  estimate  the  parameters.  Re¬ 
lated  to  the  number  of  phases  and  the  Barker  code  sequence,  the  detection  of  these 
parameters  was  possible  for  the  case  of  signal  only. 

Costas  Code  signals  were  processed  giving  the  best  results  as  a  feature  extraction 
procedure.  Parameters  such  as  the  sequence  of  the  frequency  hopping,  the  transmission 
time  per  frequency,  and  the  code  period  were  extracted  with  an  effectiveness  of  100% 
unaffected  by  the  SNR  condition  of  the  signals. 

FSK/PSK  Costas  Code  signals  were  processed  with  the  QMFB  giving  a  very  good 
effectiveness  in  detecting  all  the  important  parameters  related  to  this  signal  (the  sequence 
of  the  frequency  hopping,  the  bandwidth,  and  the  code  period)  in  the  case  of  signals 
without  noise  and  SNR=  0  dB. 
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The  FSK/PSK  Target  signal  processing  gave  a  poor  result  as  a  feature  extraction 
procedure.  The  bandwidth  was  the  only  parameter  that  could  be  extracted  with  a  good 
resolution,  but  just  for  the  case  of  signal  only.  This  is  due  mainly  to  the  complexity  of  the 
generation  of  this  kind  of  signal,  where  the  plots  generated  by  the  QMFB  processing 
looks  like  a  signal  with  a  strong  concentration  of  noise  in  a  specific  bandwidth,  making  it 
really  hard  to  identify  it. 

The  overall  architecture  of  the  QMFB  presented  in  this  thesis,  as  shown  in  Figure 
3,  is  a  good  alternative  to  an  intercept  receiver,  particularly  for  detecting  hopped  LPI 
spread  spectrum  radar  signals,  as  the  FSK/PSK  Costas  processed  in  Chapter  III. 

The  feature  extraction  gives  a  time-frequency  plane  with  good  estimation  of  posi¬ 
tion,  center  frequency  and  time  duration  of  the  signals  processed,  although  the  estimates 
of  bandwidth  degrade  in  poor  SNR.  With  respect  to  the  detection  of  phases  and  cycles  per 
phase  (or  per  bit,  depending  on  the  signal),  there  is  much  to  do  to  improve  the  resolution, 
particularly  in  signals  with  some  levels  of  noise. 

This  technique  alone  is  not  sufficient  to  process  the  multiplicity  of  available  LPI 
waveforms,  but  the  combined  use  of  this  technique  with  others,  such  as  Wigner  distribu¬ 
tion  [18],  Cyclostationary  processing  [16],  and  Higher  Order  Statistics  [15]  will  provide 
the  desired  capability. 

The  QMFB  tree  is  demonstrated  as  a  good  tool  to  process  any  basic  LPI  radar 
signal  and  it  can  be  used  on  any  other  type  of  signals  to  do  the  feature  extraction  and  that 
way  obtain  a  better  understanding  of  the  signal  characteristics. 

B.  RECOMMENDATIONS 

The  QMFB  tree  must  be  improved  to  obtain  a  better  resolution  in  detecting  the 
phases  and  cycles  per  phase  of  the  intercepted  signals. 

An  algorithm  to  determine  relevant  signal  cell  energy  features  can  analyze  the 
output  matrices  from  any  layer  of  the  QMFB.  A  list  of  these  features  could  be  sent  to  a 
classifier,  which  detennines  which  cells  belong  to  common  transmitters,  in  order  to  out¬ 
put  a  list  of  transmitter  types  and  parameters.  Further,  filtering  could  be  done  at  this  point 
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and  thus  eliminate  probable  false  alarms  and  signals  that  are  of  no  interest  to  the  intercep¬ 
tor.  Including  a  signal  classifier  in  the  analysis  is  an  obvious  next  step.  Especially  in 
some  particular  scenarios,  it  would  be  interesting  to  see  how  well  an  adaptive  classifier, 
such  as  an  artificial  neural  network,  could  use  the  information  output  from  an  analyzer 
algorithm  to  detect  and  classify  signals. 

It  is  also  instructive  to  compare  the  capabilities  and  performance  of  other  methods 
to  the  one  discussed  here.  It  may  be  possible  to  combine  some  of  the  concepts  presented 
here  with  some  of  these  methods.  For  example,  to  look  for  a  solution  to  the  problem  to 
extracting  specific  parameters  in  signals  with  added  noise,  a  High  Order  Statistic  (HOS) 
process  at  the  output  of  every  layer  in  the  QMFB  tree  structure  could  be  added  and  that 
way  try  to  filter  the  noisy  signals  getting  a  better  performance  and  resolution. 

MATFAB®  has  been  shown  to  be  a  very  useful  tool  to  develop  the  codes  to  work 
with  the  QMFB  tree  structure,  but  it  takes  some  time  to  process  the  signals,  particularly 
when  the  signals  loaded  are  long.  In  consequence,  some  efforts  must  be  focused  on  de¬ 
veloping  a  more  efficient  code  that  takes  less  time  to  do  the  processing. 
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Appendix  A:  Matlab  Programs  and  Functions 


%  startpoint.m 

%  Main  Control  Program 

% 

%  Use:  This  is  an  m  file  to  load  the  signal  that  will  be  filtered  by  the  QMFBT. 

%  The  signal  file  name  need  to  be  written  after  load  command  and  the 

%  sampling  frequency  needs  to  be  specified.  The  program  considers  the 

%  signals  storage  in  a  folder  called  "LPISignalGeneratorU". 

% 

%  Author:  Pedro  Jarpa.  Captain  Chilean  Air  Force. 

%  Date:  July  30,  2002. 

% 

clear  all;clc; 

% _ 

%  —Setting  the  signal  to  analize— 

%  Write  the  signal  directory  and  file  name  after  load  command  in  the  next  line 

load  C:\Jarpa\Thesis\LPI_Signal_Generator_II\FR_l_7_4_5_s;  %  This  is  to  load  the  signal  from 

another  directory 

fs=input('What  is  the  sampling  frequency  of  the  loadded  signal  (in  Hz)?');  %  This  parametrer  must  be  setted 
from  the  file  loaded  above 

[v,u]=size(I);  %  Since  the  QMFBT  work  with  incoming  signals  as  columns,  this  step  is  to  transform 
%  the  incoming  signal  to  a  column  vector  if  its  a  row  vector 

if  v=l 

1=1';  %  This  apply  if  the  incoming  signal  is  a  row  vector 

Q=Q';  %  This  apply  if  the  incoming  signal  is  a  row  vector 
end 

% _ 

%  —  Formating  the  Signal  —  (  Note:  This  section  apply  only  if  the  next  section  is  comented/not  used) 

%  This  section  must  be  used  only  in  the  case  of  analizing  the  signal  received  without  padding  zeros  to 
make  it  a  power  of  two. 

%  NUM_SAMPLES  =  2  ,A(log(length(I))/log(2));  %  To  generate  a  Test  Signal  of  same  number  of  samples 
as  input  Signal 
%  T=l/fs; 

%  signal=I+j*Q; 

%  tt=signal; 

% _ 

%  —  Formatting  the  Signal  — 

% 

T=l/fs;  %  Period  of  the  sampled  signal 

signal=I+j*Q;  %  Forming  the  signal 

L=length(signal);  %  Length  of  the  signal 

Lc=ceil(log2(L));  %  It  will  give  the  number  of  layers. 

z=((2.ALc)  -  L);  %  In  the  case  that  the  signal  length  does  not  correspond  to  a  power  of  two 

%  the  signal  will  be  padded  with  zeros  until  the  next  power  of  2  (2A(number  of  layers)) 
zz=zeros(l,z);  %  Vector  of  zeros  to  be  padded  to  incoming  signal  to  do  it  a  power  of  two  in  length 

tt=cat(l,  signal,  zz');  %  New  signal  vector  to  be  applyed  to  the  QMFB  resulting  from  concatenate  the  origi¬ 
nal  signal  with  the  zeros 

NUM_SAMPLES  =  length(tt);%  Length  of  the  new  signal 

% _ 

%  —  Filtering  the  Signal  with  the  QMFB  —  (qmfb  function) 
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%  qmfb(tt,'tsinc');  %  Here  the  function  qmfb  is  used,  it  will  apply  the  filter  bank  to  the  signal  with  the 
sine  modified  filter 

lay=qmfb(tt,'tsinc');  %  Here  the  function  qmfb  is  used,  it  will  apply  the  filter  bank  to  the  signal  with  the 
sine  modified  filter 

%  "lay"  is  set  equal  to  the  qmfb  function  to  get  the  number  of  the  last  layer  processed  by 

qmfb  function 

%  and  that  way  use  it  to  run  the  multiple  set  of  plot. 

%  clear  tt;  %  Optional.  Just  to  clear  the  signal. 

% _ 

%  --  Plots  — 

% 

%  Loading  the  Layers 
for  yy=4:7 

load(['c:\Jarpa\Thesis\data\layer',num2str(yy)])  ;  %  Setting  the  file  directory  and  layer's  number  in  the 
load  command 

%  Formating  the  Layer 

layer  =  abs(R  +  j*Q);  %  Taking  the  absolute  value  of  the  layer 

M  =  layer  ,A2; 

[m,n]=size(M); 

w=linspace(0,fs/2,n);  %  Setting  the  frequency  axis 

time=linspace(0,NUM_SAMPLES*T,m);  %  Setting  the  time  axis 

clear  layer  ;  %  Take  off  the  layer  matrix  from  the  workspace 

%  Ploting  the  Layer 

%  Contour  Plots 

figure(l) 

subplot(2,2,yy-3) 

contour(  time, w,M') 

title(["'Time  -  Frequency  Plane  Layer  ',num2str(yy),""]) 

xlabel('Time  (sec)'),ylabel('Frequency  (Hz)'),grid 

colorbar,colormap(jet) 

hold  on 

%  Color  Plots 

figure(2) 

subplot(2,2,yy-3) 

pcolor(time,w,M') 

title(["'Time  -  Frequency  Plane  Layer  ',num2str(yy),""]) 

xlabel('Time  (sec)'),ylabel('Frequency  (Hz)') 

colorbar,colormap(jet) 

hold  on 

%  Mesh  Plots 

figure(3) 

subplot(2,2,yy-3) 

mesh(time,w,M') 

title(['Cells  Energy  for  Layer',num2str(yy),""J) 

xlabel('Time'),ylabel('Frequency'),zlabel('Energy') 

colormap(jet) 

hold  on 

%  Surf  Plots 

figure(4) 

subplot(2,2,yy-3) 

surf(time,w,M') 

title(['Cells  Energy  for  Layer', num2str(yy),""]) 
xlabel('Time  (sec)'),ylabel('Frequency  (Hz)'),zlabel('Energy') 
colormap(jet) 
end 

% 


136 


%  qmfb.m 

%  Quadrature  Mirror  Filter  Bank  Function 

%  qmfb(f,filter,N) 

% 

%  Use: 

% 

% 

% 

% 

% 

%  Author: 

%  Date: 

% 

function  I  ay=qmfb(f,  fi  I  ter.N )  %  Quadrature  Mirror  Filter  Bank  function 

% _ 

%  —  Number  of  Layers  to  Compute  - 

n  =  floor(log2(length(f)));  %  Determine  the  amount  of  layers  from  the  length  of  the  signal 
if  nargin  <  3 ;  %  number  of  layers 

N  =  n; 

end; 

% _ 

%  —  Formatting  the  Signal  —  (to  pass  the  signal  through  the  filter  bank) 

I([l:2  .An],l)  =  f([l:2 

out  =  I; 

% _ 

%  —  Generating  the  Output  Layers  — 

%  Decompose  the  function 
for  lay  =  1  :N  %  layer 
disp(lay)  %  Show  what  layer  is  been  generated 
flag  =  1;  %  Flag  used  to  set  up  the  columns  in  the  output  matrix 

%  Reshape  the  output  matrix 
[r,c]  =  size(out); 
out  =  zeros(r./2,  c.*2); 

for  i  =  1 :2  ,A(lay-l)  %  column  of  1  (low  to  high) 

[G,F1]  =  feval(filter,  I(:,i));  %  Evaluate  the  filter  (sine  modified)  over  the  signal  going  by  columns 
%  Setting  the  output  matrix  (Layer)  by  High  Pass  (G)  and  Low  Pass  (FI)  filtering  output 


if  flag 

out(:,  i.*2-l)  =  H; 

%  Low  Pass  Branch  of  QMFB  Tree 

out(:,  i.*2)  =  G; 

%  High  Pass  Branch  of  QMFB  Tree 

else 

out(:,  i.*2-l)  =  G; 

%  High  Pass  Branch  of  QMFB  Tree 

out(:,  i.*2)  =  H; 

%  Low  Pass  Branch  of  QMFB  Tree 

end; 

flag  =  -flag; 

%  Change  in  the  flag  value 

end; 

%  Output  Matrix  and  data  to  work  with 
I  =  out;  %  Output  matrix  (signal  already  filtered) 

R  =  real(I);  %  Real  part  of  the  filtered  signal 
Q  =  imag(I);  %  Imaginary  part  of  the  filtered  signal 

save(['layer',int2str(lay)],'R','Q');  %  Saving  the  data  to  the  same  directory 

%  eval(['save  c:\Jarpa\Thesis\data\layer', int2str(lay),'.dat  I  /ascii  /double']);  %  Used  to  save  the  data  in  a 
different  directory 
end; 

% 


Quadrature  Mirror  Filter  Bank  function.  It  takes  the  input  waveform  f  (formally  the  signal 
"tt"),  and  processes  it  with  a  QMF  Bank,  'filter'  specifies  the  file  containing  the  filter 
coefficents  of  the  modified  sine  filter  to  be  applyed  through  the  QMFB  over  the  signal. 

N  is  an  optional  input  and  specifies  the  last  layer  to  be  computed.  This  can  save 
time  if  the  higher  layers  are  not  needed.  The  default  is  to  compute  all  of  the  layers. 

Pedro  Jarpa.  Captain  Chilean  Air  Force. 

July  30,  2002. 
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%  tsinc.m 

%  Modified  Sine  Filter  Function 

%  [dl,  cl]  =  tsinc(cO) 

% 

%  Use:  Modified  Sine  Filter  function.  It  use  truncated  modified  sine  filter 

%  coefficients  to  decompose  the  column  vector  cO  into  a  (high  frequency) 

%  column  vector  dl  and  (low  frequency)  column  vector  cl. 

%  tsinc  su  must  have  been  run  previously  to  create  h.dat  and 

%  g.dat  files  that  contain  the  filter  coefficients. 

% 

%  Author:  Pedro  Jarpa.  Captain  Chilean  Air  Force. 

%  Date:  July  30,  2002. 

% 

function  [dl,  cl]  =  tsinc(cO)  %  Modified  Sine  Filter  function 

% _ 

%  —  Load  Filter  Coefficients  (Sine  Modified)  — 

% 

load  h.dat;  h  =  h'; 
loadg.dat;  g  =  g'; 

% _ 

%  —  Setting  the  Signal  — 

N  =  length(cO);  %  Length  of  the  signal  (from  qmfb  function  "cO"  is  a  column) 

pad3  =  length(h)./2  +1;  %  pad  with  zeros  to  clear  out  filter 

cO  =  [cO;  zeros(pad3,l)];  %  padding  with  zeros 

% _ 

%  Decompose  the  column  vector  cO 

i  =  pad3:2:(N+pad3-2);  %  i  will  decimate  by  2 

%  cl  low  frequency  column  vector 

cl  =  filter(fliplr(h),l,cO);  %  compute  cl 

cl  =  cl (i);  %  decimate 

%  dl  high  frequency  column  vector  dl 

dl  =  filter) fliplr(g), l,c0);  %  compute  dl 

dl=dl(i);  %  decimate 

% 
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%  tsinc_su.m 

%  Tmncted  Sine  Function 

%  tsinc_su(N) 

% 

%  Use:  Truncated  Sine  function  finds  the  coefficients  for  the  truncated  sine  modified  filters. 

%  These  are  saved  in  the  data  directory  under  h.dat  (lpf)  and  g.dat  (hpf).  The  Hamming 

%  window  is  used  to  suppress  the  Gibb's  phenomena. 

%  N  is  the  number  of  coefficients  filter  pair.  For  best  results,  N  should  be  a  power  of  two. 

%  If  N  is  512  and  some  specific  values  of  S  and  C  in  the  FIR  filter  formula  yields  nearly 

%  orthogonal  filters  with  cross  correlation  of  less  than  0.001.  Due  to  that  the  FIR  filter  found 

%  is  called  "Modified  Sine  Filter". 

% 

%  Author:  Pedro  Jarpa.  Captain  Chilean  Air  Force. 

%  Date:  July  30,  2002. 

% 

function  tsinc_su(N);  %  Tmncted  Sine  Function 

% _ 

%  —  Sine  Modified  Filter  Formula  Application  — 

N  =  512;  %  Number  of  filters  coefficients  pairs. 

C  =  1.99375872328059;  %  Best  value  to  get  nearly  orthogonal  filters  with  cross  correlation  of  less  than 

0.001 

S  =  1.00618488680080;  %  Best  value  to  get  nearly  orthogonal  filters  with  cross  correlation  of  less  than 

0.001 

x  =  [-floor(N./2):floor(N./2)-l]';  %  Vector  coefficients 

h  =  sqrt(S./2).*sinc((x  +  ,5)./C);  %  Coefficients  filter  formula 

w  =  hamming(N);  %  Hamming  Window 

h  =  w.*h;  %  Low  Pass  Filter  Coefficients 

g  =  h;  %  High  Pass  Filter  Coefficients 

g([2:2:N])  =  -h([2:2:N]);  %  Alternating  the  sign  of  the  coefficients 

% _ 

%  —  Saving  the  Data  — 

save  c:\Jarpa\Thesis\data\h.dat  h  /ascii  /double;  %  Low  Pass  Filter  Coefficients 
save  c:\Jarpa\Thesis\data\g.dat  g  /ascii  /double;  %  High  Pass  Filter  Coefficients 

% 
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%  multiple_colorplot.m 

%  Time-frequency  Plot  using  “pcolor”  command 

% 

%  Use:  Script  file  to  display  a  color  plot  of  the  output  layers  from  the  QMFB  tree. 

%  Needs  to  write  the  name  of  the  signal  in  the  “save”  command  at  the  end  of  the  code. 

% 

%  Author:  Pedro  Jarpa.  Captain  Chilean  Air  Force. 

%  Date:  July  30,  2002. 

% 

% _ 

%  —  Load  Layer  — 

%  Write  the  file  directory  and  layer's  number  after  load  command  in  the  next  line 

yy=0; 

counter- 1 ; 
for  yy=2:lay-l 

load(['c:\Jarpa\Thesis\data\layer',num2str(yy)])  ;  %  Setting  the  file  directory  and  layer's  number  in  the 
load  command 

%  Formating  the  Layer 

layer  =  abs(R  +  j*Q);  %  Taking  the  absolute  value  of  the  layer 

layer  =  layer/max(max(layer));  %  Normalizing  the  layer 

M  =  layer  ,A2; 

[m,n]=size(M); 

w=linspace(0,fs/2,n);  %  Setting  the  frequency  axis 

time=linspace(0,NUM_SAMPLES*T,m);  %  Setting  the  time  axis 

clear  layer;  %  Take  off  the  layer  matrix  from  the  workspace 

%  Color  Plots 

figure(counter) 

pcolor(time,w,M') 

title(['"Time  -  Frequency  Plane  Layer  ',num2str(yy),""]) 
xlabel('Time  (sec)'),ylabel('Frequency  (Hz)') 
colorbar,colormap(jet) 

saveas(gcf,['c:\Jarpa\Thesis\plots\FSK_PSK_T_15_256_10_s_colorplot_layer_',num2str(yy)],'fig'); 

%  CAREFULL! ! !  . .  ..The  name  of  the  signal  needs  to  be  changed! ! ! ! ! 

count  er=counter+ 1 ; 
end 
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%  multiple_contourplot.m 

%  Time-frequency  Plot  using  “contour”  command 

% 

%  Use:  Script  file  to  display  a  contour  plot  of  the  output  layers  from  the  QMFB  tree. 

%  Needs  to  write  the  name  of  the  signal  in  the  “save”  command  at  the  end  of  the  code. 

% 

%  Author:  Pedro  Jarpa.  Captain  Chilean  Air  Force. 

%  Date:  July  30,  2002. 

% 

% _ 

%  —  Load  Layer  — 

%  Write  the  file  directory  and  layer's  number  after  load  command  in  the  next  line 

yy=0; 

counter- 1 ; 
for  yy=2:lay-l 

load(['c:\Jarpa\Thesis\data\layer',num2str(yy)])  ;  %  Setting  the  file  directory  and  layer's  number  in  the 
load  command 

%  Formating  the  Layer 

layer  =  abs(R  +  j*Q);  %  Taking  the  absolute  value  of  the  layer 

layer  =  layer/max(max(layer));  %  Normalizing  the  layer 

M  =  layer  ,A2; 

[m,n]=size(M); 

w=linspace(0,fs/2,n);  %  Setting  the  frequency  axis 

time=linspace(0,NUM_SAMPLES*T,m);  %  Setting  the  time  axis 

clear  layer;  %  Take  off  the  layer  matrix  from  the  workspace 

%  Contour  Plots 

figure(counter) 

contour(time,w,M') 

title(['"Time  -  Frequency  Plane  Layer  ',num2str(yy),""]) 
xlabel('Time  (sec)'),ylabel('Frequency  (Hz)'),grid 
colorbar,colormap(jet) 

saveas(gcf, ['c:\Jarpa\Thesis\plots\FSK_PSK_T_15_256_10_s_contourplot_layer_', num2str(yy)], 'fig'); 

%  CAREFULL...! ! .  ..the  name  of  the  signal  needs  to  be  changed! ! ! 

count  er=counter+ 1 ; 
end 
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%  multiple_  meshplot.m 

%  3  Dimension  View  Plot  using  “mesh”  command 

% 

%  Use:  Script  file  to  display  a  mesh  plot  of  the  output  layers  from  the  QMFB  tree. 

%  Needs  to  write  the  name  of  the  signal  in  the  “save”  command  at  the  end  of  the  code. 

% 

%  Author:  Pedro  Jarpa.  Captain  Chilean  Air  Force. 

%  Date:  July  30,  2002. 

% 

yy=0; 

counter=l; 
for  yy=2:lay-l 

load(['c:\Jarpa\Thesis\data\layer',num2str(yy)])  ;  %  Setting  the  file  directory  and  layer's  number  in  the 
load  command 

%  Formating  the  Layer 

layer  =  abs(R  +  j*Q);  %  Taking  the  absolute  value  of  the  layer 

layer  =  layer/max(max(layer));  %  Normalizing  the  layer 

M  =  layer  ,A2; 

[m,n]=size(M); 

w=linspace(0,fs/2,n);  %  Setting  the  frequency  axis 

time=linspace(0,NUM_SAMPLES*T,m);  %  Setting  the  time  axis 

clear  layer;  %  Take  off  the  layer  matrix  from  the  workspace 

%  Mesh  Plots 

figure(counter) 

%  mesh(time,w,M') 

mesh(time,w,M'),view(90,0)  %  To  see  the  frequency  axis 

title(['Cells  Energy  for  Layer',num2str(yy),""]) 
xlabel('Time  (sec)'),ylabel('Frequency  (Hz)'),zlabel('Energy') 
colormap(jet) 

saveas(gcf, ['c:\Jarpa\Thesis\plots\FSK_PSK_C_l_15_l  l_5_0_meshplot_layer_',num2str(yy)], 'fig'); 

%  CAREFULL. . . ! ! . .  .the  name  of  the  signal  needs  to  be  changed ! ! ! ! ! 

count  er=counter+ 1 ; 
end 
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%  multiple_  surfplot.m 

%  3  Dimension  View  Plot  using  “surf”  command 

% 

%  Use:  Script  file  to  display  a  surf  plot  of  the  output  layers  from  the  QMFB  tree. 

%  Needs  to  write  the  name  of  the  signal  in  the  “save”  command  at  the  end  of  the  code. 

% 

%  Author:  Pedro  Jarpa.  Captain  Chilean  Air  Force. 

%  Date:  July  30,  2002. 

% 

yy=0; 

counter=l; 
for  yy=2:lay-l 

load(['c:\Jarpa\Thesis\data\layer',num2str(yy)])  ;  %  Setting  the  file  directory  and  layer's  number  in  the 
load  command 

%  Formating  the  Layer 

layer  =  abs(R  +  j*Q);  %  Taking  the  absolute  value  of  the  layer 

layer  =  layer/max(max(layer));  %  Normalizing  the  layer 

M  =  layer  ,A2; 

[m,n]=size(M); 

w=linspace(0,fs/2,n);  %  Setting  the  frequency  axis 

time=linspace(0,NUM_SAMPLES*T,m);  %  Setting  the  time  axis 

clear  layer;  %  Take  off  the  layer  matrix  from  the  workspace 

%  Mesh  Plots 

figure(counter) 

surf(time,w,M') 

%  surf(time,w,M'),view(90,0)  %  To  see  the  frequency  axis 

title(['Cells  Energy  for  Layer',num2str(yy),""]) 
xlabel('Time  (sec)'),ylabel('Frequency  (Hz)'),zlabel('Energy') 
colormap(jet) 

saveas(gcf, ['c:\Jarpa\Thesis\plots\FSK_PSK_C_l_15_l  l_5_0_surfplot_layer_',num2str(yy)], 'fig'); 

%  CAREFULL. . . ! ! . .  .the  name  of  the  signal  needs  to  be  changed ! ! ! ! ! 

count  er=counter+ 1 ; 
end 
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Appendix  B:  List  of  Signals  Generated  by  LPIG  and  Processed  by 

QMFB. 


TEST  SIGNALS 


Test  Signal 

Carrier  Frequeny 

SNR 

1 

T  1  7  1  s 

1000  Hz 

Signal  Only 

2 

T  2  7  1  s 

1000  Hz,  2000  Hz 

Signal  Only 

BINARY  PHASE  SHIFT  KEYING  (BPSK) 


BPSK 

Number  of  bits 
per  Barker  code 
(Bits) 

Number  of  cy¬ 
cles  per  bit  (cpb) 

SNR 

1 

B  1  7  7  1  s 

7 

1 

Signal  Only 

2 

B  1  7  7  1  0 

7 

1 

OdB 

3 

B  1  7  7  1  -6 

7 

1 

-6  dB 

4 

B  1  7  11  1  s 

11 

1 

Signal  Only 

5 

B  1  7  11  1  0 

11 

1 

OdB 

6 

B  1  7  11  1  -6 

11 

1 

-6  dB 

7 

B  1  7  7  5  s 

7 

5 

Signal  Only 

8 

B  1  7  7  5  0 

7 

5 

OdB 

9 

B  1  7  7  5  -6 

7 

5 

-6  dB 

10 

B  1  7  11  5  s 

11 

5 

Signal  Only 

11 

B  1  7  11  5  0 

11 

5 

OdB 

12 

B  1  7  11  5  -6 

11 

5 

-6  dB 

FREQUENCY  MODULATED  CONTINUOUS  WAVE  (FMCW) 


FMCW 

Modulation 

Bandwidth 

Modulation 

Period 

SNR 

1 

F  1  7  250  20  s 

250  Hz 

20  ms 

Signal  Only 

2 

F  1  7  250  20  0 

250  Hz 

20  ms 

OdB 

3 

F  1  7  250  20  -6 

250  Hz 

20  ms 

-6  dB 

4 

F  1  7  250  30  s 

500  Hz 

30  ms 

Signal  Only 

5 

F  1  7  250  30  0 

500  Hz 

30  ms 

OdB 

6 

F  1  7  250  30  -6 

500  Hz 

30  ms 

-6  dB 

7 

F  1  7  500  20  s 

250  Hz 

20  ms 

Signal  Only 

8 

F  1  7  500  20  0 

250  Hz 

20  ms 

OdB 

9 

F  1  7  500  20  -6 

250  Hz 

20  ms 

-6  dB 

10 

F  1  7  500  30  s 

500  Hz 

30  ms 

Signal  Only 

11 

F  1  7  500  30  0 

500  Hz 

30  ms 

OdB 

12 

F  1  7  500  30  -6 

500  Hz 

30  ms 

-6  dB 
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FRANK  CODES 


FRANK 

Number  of 
code  phases 
N 

Number  of  cycles 
per  phase 
cpp 

SNR 

1 

FR  1  7  4  1  s 

4 

1 

Signal  Only 

2 

FR  1  7  4  1  0 

4 

1 

OdB 

3 

FR  1  7  4  1  -6 

4 

1 

-6  dB 

4 

FR  1  7  4  5  s 

4 

5 

Signal  Only 

5 

FR  1  7  4  5  0 

4 

5 

OdB 

6 

FR  1  7  4  5  -6 

4 

5 

-6  dB 

7 

FR  1  7  8  1  s 

8 

1 

Signal  Only 

8 

FR  1  7  8  1  0 

8 

1 

OdB 

9 

FR  1  7  8  1  -6 

8 

1 

-6  dB 

10 

FR  1  7  8  5  s 

8 

5 

Signal  Only 

11 

FR  1  7  8  5  0 

8 

5 

OdB 

12 

FR  1  7  8  5  -6 

8 

5 

-6  dB 

PI  POLYPHASE  CODE 


PI 

Number  of 
code  phases 
N 

Number  of  cycles 
per  phase 
cpp 

SNR 

1 

PI  1  7  4  1  s 

4 

1 

Signal  Only 

2 

PI  1  7  4  1  0 

4 

1 

OdB 

3 

PI  1  7  4  1  -6 

4 

1 

-6  dB 

4 

PI  1  7  4  5  s 

4 

5 

Signal  Only 

5 

PI  1  7  4  5  0 

4 

5 

OdB 

6 

PI  1  7  4  5  -6 

4 

5 

-6  dB 

7 

PI  1  7  8  1  s 

8 

1 

Signal  Only 

8 

PI  1  7  8  1  0 

8 

1 

OdB 

9 

PI  1  7  8  1  -6 

8 

1 

-6  dB 

10 

PI  1  7  8  5  s 

8 

5 

Signal  Only 

11 

PI  1  7  8  5  0 

8 

5 

OdB 

12 

PI  1  7  8  5  -6 

8 

5 

-6  dB 
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P2  POLYPHASE  CODE 


P2 

Number  of 
code  phases 
N 

Number  of  cycles 
per  phase 
cpp 

SNR 

1 

P2  1  7  4  1  s 

4 

1 

Signal  Only 

2 

P2  1  7  4  1  0 

4 

1 

OdB 

3 

P2  1  7  4  1  -6 

4 

1 

-6  dB 

4 

P2  1  7  4  5  s 

4 

5 

Signal  Only 

5 

P2  1  7  4  5  0 

4 

5 

OdB 

6 

P2  1  7  4  5  -6 

4 

5 

-6  dB 

7 

P2  1  7  8  1  s 

8 

1 

Signal  Only 

8 

P2  1  7  8  1  0 

8 

1 

OdB 

9 

P2  1  7  8  1  -6 

8 

1 

-6  dB 

10 

P2  1  7  8  5  s 

8 

5 

Signal  Only 

11 

P2  1  7  8  5  0 

8 

5 

OdB 

12 

P2  1  7  8  5  -6 

8 

5 

-6  dB 

P3  POLYPHASE  CODE 


P3 

Number  of 
code  phases 
N2 

Number  of  cycles 
per  phase 
cpp 

SNR 

1 

P3  1  7  16  1  s 

16 

1 

Signal  Only 

2 

P3  1  7  16  1  0 

16 

1 

OdB 

3 

P3  1  7  16  1  —6 

16 

1 

-6  dB 

4 

P3  1  7  16  5  s 

16 

5 

Signal  Only 

5 

P3  1  7  16  5  0 

16 

5 

OdB 

6 

P3  1  7  16  5  -6 

16 

5 

-6  dB 

7 

P3  1  7  64  1  s 

64 

1 

Signal  Only 

8 

P3  1  7  64  1  0 

64 

1 

OdB 

9 

P3  1  7  64  1  -6 

64 

1 

-6  dB 

10 

P3  1  7  64  5  s 

64 

5 

Signal  Only 

11 

P3  1  7  64  5  0 

64 

5 

OdB 

12 

P3  1  7  64  5  -6 

64 

5 

-6  dB 
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P4  POLYPHASE  CODE 


P4 

Number  of 
code  phases 
N2 

Number  of  cycles 
per  phase 
cpp 

SNR 

1 

P4  1  7  16  1  s 

16 

1 

Signal  Only 

2 

P4  1  7  16  1  0 

16 

1 

OdB 

3 

P4  1  7  16  1  -6 

16 

1 

-6  dB 

4 

P4  1  7  16  5  s 

16 

5 

Signal  Only 

5 

P4  1  7  16  5  0 

16 

5 

OdB 

6 

P4  1  7  16  5  -6 

16 

5 

-6  dB 

7 

P4  1  7  64  1  s 

64 

1 

Signal  Only 

8 

P4  1  7  64  1  0 

64 

1 

OdB 

9 

P4  1  7  64  1  -6 

64 

1 

-6  dB 

10 

P4  1  7  64  5  s 

64 

5 

Signal  Only 

11 

P4  1  7  64  5  0 

64 

5 

OdB 

12 

P4  1  7  64  5  -6 

64 

5 

-6  dB 

COSTAS  CODE 


COSTAS 

Sequence 

Number  of  cycles 
per  phase 
cpp 

SNR 

1 

C  1  15  10  s 

4716523 

10 

Signal  Only 

2 

C  1  15  10  0 

4716523 

10 

OdB 

3 

C  1  15  10  -6 

4716523 

10 

-6  dB 

4 

C  1  15  20  s 

4716523 

20 

Signal  Only 

5 

C  1  15  20  0 

4716523 

20 

OdB 

6 

C  1  15  20  -6 

4716523 

20 

-6  dB 

7 

C  2  17  10  s 

2638751 

10 

Signal  Only 

8 

C  2  17  10  0 

2638751 

10 

OdB 

9 

C  2  17  10  -6 

2638751 

10 

-6  dB 

10 

C  2  17  20  s 

2638751 

20 

Signal  Only 

11 

C  2  17  20  0 

2638751 

20 

OdB 

12 

C  2  17  20  -6 

2638751 

20 

-6  dB 
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FREQUENCY  SHIFT  KEYING/PHASE  SHIFT  KEYING  COMBINED  WITH 
COSTAS  CODE  (FSK/PSK  COSTAS) 


FSK/PSK  Costas 

Sequence 

Barker 

Bits 

Number  of 
cycles  per  bit 
cpp 

SNR 

1 

FSK  PSK  C  1  15  5  1  s 

4716523 

5 

1 

Signal  Only 

2 

FSK  PSK  C  1  15  5  1  0 

4716523 

5 

1 

OdB 

3 

FSK  PSK  C  1  15  5  5  s 

4716523 

5 

5 

Signal  Only 

4 

FSK  PSK  C  1  15  5  5  0 

4716523 

5 

5 

OdB 

5 

FSK  PSK  C  1  15  11  1  s 

4716523 

11 

1 

Signal  Only 

6 

FSK  PSK  C  1  15  11  1  0 

4716523 

11 

1 

OdB 

7 

FSK  PSK  C  1  15  11  5  s 

4716523 

11 

5 

Signal  Only 

8 

FSK  PSK  C  1  15  11  5  0 

4716523 

11 

5 

OdB 

FSK/PSK  COMBINED  WITH  TARGET-MATCHED  FREQUENCY  HOPPING 


FSK/PSK  Target 

Number  of 
random  hops 
in  frequency 

Number  of 
cycles  and 
phases 

SNR 

1 

FSK  PSK  T  15  128  5  s 

128 

5 

Signal  Only 

2 

FSK  PSK  T  15  128  5  0 

128 

5 

OdB 

3 

FSK  PSK  T  15  128  10  s 

128 

10 

Signal  Only 

4 

FSK  PSK  T  15  128  10  0 

128 

10 

OdB 

5 

FSK  PSK  T  15  256  5  s 

256 

5 

Signal  Only 

6 

FSK  PSK  T  15  256  5  0 

256 

5 

OdB 

7 

FSK  PSK  T  15  256  10  s 

256 

10 

Signal  Only 

8 

FSK  PSK  T  15  256  10  0 

256 

10 

OdB 
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